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20.  ABSTRACT  (concluded) 

without  proper  scrutiny  of  the  basis  for  the  conclusions  reached.  Some 
of  these  Inconsistencies  are  pointed  out  and  the  various  analytical 
approaches  used  in  describing  solid  particle  erosion  phenomena  are 
suaaurlzed. 


The  need  for  an  adequately  characterized  solid  particle  erosion 
testing  capability  Is  Indicated  due  to  the  range  of  test  parameters 
which  have  been  found  recently  to  significantly  influence  the  measured 
erosion  rates.  Meaningful  advancement  in  the  field  is  strongly  depen- 
dent on  the  development  of  a reliable  experimental  data  base  and  experi- 
mentation which  Isolates  and  clarifies  a particular  feature  of  the 
material  removal  process. 
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1.0  INTRODUCTION 


Within  the  last  decade  a number  of  Investigators  have  contributed 
to  a clearer  understanding  of  the  solid  particle  erosion  mechanisms  in 
both  ductile  and  brittle  materials.  The  erosion  process  can  be  qualita- 
tively described  for  several  materials,  however  most  theories  of  solid 
particle  erosion  are  based  on  material  removal  models  which  do  not  repre- 
sent the  actual  physical  processes  of  material  removal  as  they  are  pres- 
ently known.  It  became  evident  to  the  author  that  a variety  of  published 
results  have  been  generally  accepted  by  the  erosion  community  without 
proper  scrutiny  of  the  basis  for  the  conclusions  reached.  The  objectives 
of  this  study  are  to  re-examine  both  productive  and  non-productive  in- 
vestigations which  have  been  undertaken  In  the  field  of  solid  particle 
erosion  In  order  to  assess  their  acttial  contribution  with  respect  to  the 
current  state  of  the  art  and  to  summarize  the  analyses  of  the  particle 
Impact  process  and  target  response  which  have  been  pursued  as  a guide 
for  future  activities  of  this  type.  This  report  concentrates  on  the 
analysis  and  modeling  of  solid  particle  erosion  phenomena,  since  this 
aspect  of  the  subject  does  not  appear  to  have  been  adequately  addressed 
In  the  several  reviews  which  have  been  published  (Wahl  and  Hartsteln, 
1946;  Welllnger  and  Uetz,  1955;  Heymann,  1968;  Engel,  1976;  Preece  and 
Macmillan,  1977;  Evans,  1979;  Ruff  and  Wlederhorn,  1979). 

The  solid  particle  erosion  literature  spans  many  areas  of  applica- 
tion and  has  received  the  attention  of  Investigators  with  diverse  back- 
grounds. Motivated  for  the  most  part  by  the  need  for  evaluating  the  re- 
sponse of  materials  and  coatings  to  solid  particle  Impingement  rather 
than  In  understanding  the  material  removal  process  Itself,  a diversity 
of  test  configurations  have  been  utilized  to  generate  the  screening 
assessments  required.  Erosion  testing  has  generally  been  undertaken  to 
satisfy  a very  specific  need  for  a particular  application.  Typically  a 
test  method  Is  devised  representative  of  the  operational  environment. 

A material  selection  procedure  is  pursued  and  then  the  testing  Is 
terminated.  It  Is  difficult  to  construct  a coherent  series  of  test 


results  on  this  basis.  On  "le  other  hand  there  have  been  a few  programs 
which  have  been  carried  out  to  Investigate  the  erosion  process;  these 
have  been  referred  to  fairly  extensively  In  the  erosion  literature. 
Unfortunately,  the  available  solid  particle  erosion  test  results  have 
generally  been  reported  with  Insufficient  Information  pertaining  to  the 
test  conditions  to  provide  a valid  comparison  between  the  data  obtained 
In  one  laboratory  with  that  from  another  laboratory.  This  situation  Is 
complicated  further  by  the  fact  that  changes  In  the  construction  of  the 
same  test  apparatus  can  significantly  Influence  the  delivery  of  particles 
to  the  target  material  as  shown  by  Uuemols  and  Klels  (1975),  Wolak,  et 
al.  (1977),  and  Majl  and  Sheldon  (1979). 


An  effort  was  made  to  organize  and  compile  the  various  types  of 
test  apparatus  used,  the  test  procedures,  and  the  trends  In  the  test 
results,  however  It  was  found  that  too  many  voids  existed  In  most  publi- 
cations pertaining  to  both  the  description  of  the  apparatus  used  and 
Information  on  conducting  the  tests  to  make  this  a meaningful  undertaking 
Therefore  only  a limited  number  of  the  data  trends  In  the  available  ero- 
sion data  are  summarized  In  Section  2 which  are  relevant  to  the  detailed 
evaluations  of  the  analytical  treatments  of  solid  particle  erosion  In 
Section  3. 


The  scope  of  this  work  Includes  all  classes  of  materials  with  some 
emphasis  on  metals  and  ceramics  exposed  to  multiple  particle  Impacts 
ranging  from  less  than  a meter  per  second  (ms  ^)  to  an  upper  limit  around 
500  ms  ^ which  encompasses  most  of  the  work  on  solid  particle  Impinge- 
ment of  turbine  engines,  helicopter  rotor  blades,  aircraft  window  mate- 
rials, pneumatic  transport  of  solid  materials,  coal  conversion  processes, 
and  a variety  of  additional  commonly  encountered  conditions  for  which 
solid  particle  erosion  Is  an  Important  factor  In  operational  performance. 
The  historical  development  of  the  subject  Is  of  some  Interest,  but  It  Is 
not  the  purpose  of  this  report  to  provide  a complete  review  of  the 
available  literature  or  to  follow  an  historical  perspective. 


2.0  GENERAL  OBSERVATIONS 


Early  test  devices  have  been  descr.  ‘^ed  by  Wahl  and  Hartsteln  (1946). 
They  divided  the  test  devices  used  by  various  investigators  into  five 
categories:  sand  blast  equipment,  wear  chambers,  wear  pumps,  wear  noz- 
zles, and  other  jet  devices  which  pertained  primarily  to  water  jets. 

The  use  and  representative  examples  of  each  of  these  methods  for  pro- 
ducing wear  are  described  by  Wahl  and  Hartsteln. 

The  most  widely  used  test  configuration  for  solid  particle  erosion 
is  the  sand  blast  apparatus.  Although  some  version  of  this  device  has 
been  used  for  laboratory-scale  erosion  evaluations  for  over  seventy 
years,  it  has  only  been  quite  recently  that  the  test  parameters  associated 
with  particular  types  of  equipment  have  been  investlgaged.  The  earliest 
reference  provided  by  Wahl  and  Hartsteln  (1946)  for  laboratory  sand  blast 
equipment  is  to  the  work  of  Gary  in  1904  (Gary,  1904).  Thi  laboratory 
equipment  used  for  providing  a multiple  particle  field  generally  worked 
on  the  basis  of  ingestion  of  particulate  matter  into  a flowing  gas  stream 
which  subsequently  impinges  on  a test  piece.  The  Impact  velocity  and 
angle  of  attack  of  the  particles  in  the  gas  stream  are  governed  in  part 
by  their  size,  shape,  and  density. 

Solid  particles  can  also  be  injected  into  a flowing  gas  stream  as 
in  a wind  tunnel.  The  specimen  is  then  completely  engulfed  by  the  par- 
ticle-laden gas  which  offers  the  advantage  of  eliminating  the  transition 
region  around  the  eroded  spot  on  the  test  piece  when  the  particles  are 
propelled  from  a nozzle  or  tube  in  the  sand  blast  arrangement.  A well- 
characterized  system  of  this  type  has  been  described  by  Grant  and 
Tabakoff  (1975)  and  Tabakoff  and  Wakeman  (1979), 

Rotating  arm  devices  in  which  the  specimen  is  mounted  at  the  top 
of  a rotor  blade  and  then  passes  through  a stream  of  falling  particles 
are  also  used  for  solid  particle  erosion  evaluations  (Klels,  1969; 
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Tilly  and  Sage,  1970;  Adler,  Morris,  and  Wahl,  1972).  The  rotating  arm 
configuration  has  the  advantage  of  having  the  specimen  Impact  a graded 
particle  distribution  '•'uch  as  specified  by  U.S.  Military  Specification, 
Mil-E-5007C,  which  Includes  particle  sizes  over  the  range  from  0 to 
1000  ym)  at  a uniform  velocity.  Such  conditions  are  characteristic  of 
high  speed  structural  components  In  a dust  field  where  the  particles  are 
traveling  at  considerably  lower  velocities  as  for  a helicopter  rotor. 

The  rotating  arm  configurations  also  have  potential  for  being  operated 
at  higher  Impact  velocities  compared  with  sand  blast  and  wind  tunnel 
test  facilities. 

The  variations  which  exist  In  solid  particle  erosion  test  apparatus, 
even  for  sand  blast  devices,  make  it  difficult  to  discuss  typical,  or 
standard,  test  conditions.  Generally  the  test  results  are  strongly 
related  to  the  test  apparatus  used.  In  addition  the  testing  procedures 
for  solid  particle  erosion  vary  widely  from  laboratory  to  laboratory,  so 
it  is  difficult  to  make  an  absolute  comparison  between  the  results  ob- 
tained for  one  set  of  test  conditions  to  another.  Also,  as  more  atten- 
tion is  paid  to  the  test  procedures  and  characterization  of  the  impact 
conditions,  additional  factors  which  influence  the  impacting  particle/ 
target  interactions  are  being  discovered  (Tilly,  1969;  Tilly  and  Sage, 
1970;  Uuemols  and  Kleis,  1975;  Grant  and  Tabakoff,  1975;  Wolak,  et  al. , 
1977;  Young  and  Ruff,  1977;  Tabakoff,  Kotwal,  and  Hamed,  1979;  Maji  and 
Sheldon,  1979;  Lapldes  and  Levy,  1979),  Therefore  at  the  present  time 
the  test  facility  and  the  test  procedure  must  be  considered  In  conjunc- 
tion with  the  test  results. 

A considerable  amount  of  testing  has  been  devoted  to  the  deter- 
mination of  the  so-called  "erosion  resistance"  of  a material.  Since, 
as  will  be  demonstrated  in  Sections  3 and  4,  no  satisfactory  analysis 
of  solid  particle  erosion  exists.  It  is  not  possible  to  provide  an  ex- 
pression representing  the  erosion  resistance  of  a particular  material 
which  can  be  used  In  a general  context  involving  different  particulates. 
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Impact  velocities,  angles  of  attack,  and  so  on.  Unfortunately,  without 
an  adequate  analysis  of  solid  particle  erosion,  the  erosion  resistance 
of  a material  according  to  a specific  definition  can  only  be  compared 
under  the  same  test  conditions  with  another  material.  Thus,  one  finds 
that  most  of  the  erosion  testing  reported  In  the  literature  Is  a ranking 
of  materials  for  a specific  environmental  condition.  Furthermore,  since 
the  majority  of  solid  particle  erosion  tests  have  been  for  the  purpose 
of  ranking  and  comparing  the  erosion  behavior  of  materials,  there  has 
been  little  concern  for  quantitatively  specifying  the  test  parameters  or 
providing  an  adequate  description  of  the  test  apparatus.  Faced  with  this 
situation  several  Investigators  have  commented  on  the  general  lack  of 
useable  test  results  In  the  field  of  solid  particle  erosion. 

The  general  erosion  behavior  of  a generic  material  eroded  under  a 
specific  test  condition  has  the  characteristic  response  Illustrated  In 
Figure  2.1.  The  specimen  Is  weighed  before  and  after  exposure  to  the 
eroding  environment.  The  mass  loss  as  a function  of  the  exposure  to  the 
eroding  environments  Is  the  raw  data  obtained  from  a test.  The  mass  loss/ 
exposure  time  curve  has  a general  form  which  consists  of  an  Incubation 
period,  acceleration  period,  and  the  maximum  rate  period.  Although 
stages  beyond  the  maximum  rate,  or  steady  state  rate,  of  material  removal 
have  been  defined,  the  material  at  this  time  Is  so  badly  eroded  that  Its 
usefulness  for  most  applications  Is  no  longer  of  Interest.  Most  of  the 
evaluations  made  for  solid  particle  erosion  continue  the  test  until  the 
maximum  rate  of  erosion  has  been  established.  This  Is  accomplished  by 
running  the  test  for  successive  equal  time  Increments  and  obtaining  the 
same  value  of  the  mass  loss  for  these  Increments.  Most  of  the  results 
reported  for  solid  particle  erosion  evaluations  are  expressed  as  the 
ratio  of  the  mass  of  material  removed  to  the  mass  of  the  Impacting  par- 
ticles for  a known  exposure  time  Increment.  The  ratio  of  the  volume  of 
material  removed  to  the  mass  of  the  Impacting  particles  or  the  volume 
removed  per  particle  Impact  are  also  used  as  measures  of  the  erosion  rate. 
Two  types  of  response  to  multiple  solid  particle  Impacts  have  been  noted. 


MASS  LOSS 


EXPOSURE  TIME 


Figure  2.1.  Characteristic  Maas  Removal  Curve  in  Terns  of  the  Tine 
of  Exposure  to  the  Erosive  Environment. 


The  first  Is  a ductile  response  generally  associated  with  metals,  while 
the  second  Is  a brittle  response  associated  with  glasses  and  ceramics. 

The  general  form  of  the  mass  removal  rate  as  a function  of  the  angle  of 
attack  for  each  response  Is  shown  In  Figure  2.2. 

It  Is  a general  observation,  however,  that  normally  brittle 
materials  exhibit  ductile  behavior  when  the  particle  diameters  are  less 
than  100  ym.  Sheldon  and  Finnic  (1966a)  state  that  this  phenomenon  occurs 
when  the  zone  of  material  affected  by  an  Impacting  particle  is  smaller 
than  the  spacing  between  material  flaws  or  Inhomogeneltles.  Experimen- 
tally this  transition  from  brittle  to  ductile  behavior  is  shown  by  plot- 
ting the  erosion  rate  versus  impact  angle  relationship  for  different 
sizes  of  Impacting  particles.  For  larger  particles  the  typical  brittle 
curve  is  found  with  its  maximum  at  normal  incidence.  For  very  small  par- 
ticles the  same  material  may  show  a dependence  more  like  the  typical 
ductile  curve  having  a maximum  at  other  than  normal  incidence.  Another 
Indication  of  ductile  behavior  Is  the  formation  of  surface  ripples  which 
Is  characteristic  of  ductile  erosion  around  the  impact  angles  at  which 
the  maximum  rate  of  erosion  occurs  on  brittle  materials  eroded  by  very 
small  particles. 

The  formation  of  ripples,  often  very  distinct  and  regular,  has 
been  reported  by  several  investigators.  Ripples  occur  only  for  ductile 
response  and  only  at  impact  angles  which  are  oblique.  Finnie  and  Kabil 
(1965)  have  proposed  a theoretical  explanation  in  terms  of  the  interac- 
tion between  surface  irregularities  and  the  equations  of  motion  of  im- 
pacting particles.  They  show  that  peaks  will  be  attacked  less  severely 
than  valleys  so  that  irregularities  tend  to  grow.  Initially  the  irregu- 
larities will  be  random,  but  a certain  wavelength  is  most  favored  for 
growth  and  ripples  at  that  wavelength  will  eventually  be  established. 

The  ripple  pattern  is  predicted  to  move  downstream  as  erosion  progresses. 
Moore  (1968)  showed  that  conical  brass  targets  with  a 90”  cone  angle  de- 
velop almost  perfect  right-angled  steps,  whose  width  and  pitch  diminish 
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with  Increasing  distance  from  the  cone  apex  when  Impacted  by  glass  beads. 
He  concluded  from  photomicrographs  that  the  step  formation  Is  a metal 
flow  phenomenon. 

There  are  many  factors  which  Influence  the  material  removal  process; 
however  the  particle  size  exponent  and  the  velocity  exponent  will  be  of 
primary  concern,  since  these  parameters  are  common  to  the  majority  of  the 
solid  particle  erosion  analyses  which  can  be  found  In  the  literature. 

The  lack  of  agreement  between  the  results  of  the  theories  based  on  kinetic 
energy  exchange  between  the  Impacting  particle  and  the  target  has  been  a 
perplexing  Issue  In  the  field  of  solid  particle  erosion.  The  general  form 
of  the  expression  for  volume  removal  (volume  removed  per  particle  Impact) 
can  be  written, 

U - g V®  R^’fCa)  (2.1) 

o 

where  g Is  a functional  representing  the  dependence  on  all  of  the  param- 
eters not  explicitly  stated 


Is  the  particle  Impact  velocity 
R Is  the  radius  of  the  Impacting  particle 

f(a)  Is  a function  of  the  attack  angle  a 

a Is  the  velocity  exponent 

b Is  the  particle  size  exponent 

The  velocity  exponent  has  been  generally  found  to  range  from  two  to 
three  for  materials  displaying  ductile  behavior  for  a wide  range  of  par- 
ticle Impact  conditions,  however  higher  exponents  (up  to  4)  have  been  ob- 
tained by  Grant  and  Tabakoff  (1975)  and  Tabakoff,  et  al.  (1979).  Taba- 
koff's  results  are  based  on  detailed  observations  of  the  aerodynamics  of 
the  particle  flow  and  the  reboimd  characteristics  of  the  particles  Im- 
pacting the  eroding  surface  of  the  target.  The  velocity  exponents  for 
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brittle  materials  range  from  1 to  6.5  (cp.  Gulden,  1979a;  Sargent,  et 

al.,  1979).  The  velocity  exponent  depends  on  the  size,  material,  shape,  . j 

and  angle  of  attack  of  the  Impacting  particle,  however  It  does  not  appear 

to  be  significantly  affected  by  the  microstructure  of  the  target  material.  _ i 

The  particle  size  exponent  has  not  been  adequately  evaluated.  It  is  h 

generally  stated  that  the  volume  removal  In  metals  Is  Independent  of  par- 
ticle size  for  particle  sizes  In  excess  of  100  pm  (Tilly,  1969).  The 
variations  In  particle  size  effects  has  already  been  noted  for  smaller 
particle  sizes  in  brittle  materials  (Sheldon  and  Flnnle,  1966a).  Gulden's 
test  data  (Gulden,  1979a,b)  shows  bS4  for  magnesium  fluoride  Impacted  by 
quartz  particles,  however  hot-pressed  silicon  nitride  Impacted  by  silicon  ] 

carbide  particles  displayed  a b~4  dependence  while  Impacts  by  quartz 
produced  a b>3  dependence  as  well  as  a change  In  a*4  to  a-1  for  par- 
ticle sizes  from  8 to  940  pm.  A definite  threshold  for  achieving  an 

4 4 < 

R relation  was  found  for  reaction-bonded  silicon  nitride  Impacted  by 

quartz  particles.  The  above  results  are  for  normal  Impacts  on  the  target 

surface.  The  dependence  of  the  particle  size  exponent  on  the  particle 

material  Is  evident,  but  It  Is  not  known  how  It  would  vary  with  the  angle 

of  attack. 

The  complete  characterization  of  a particular  Impact  condition 
(particle  type  and  target  material)  Involving  the  particle  size.  Impact 
velocity,  and  Impingement  angle  requires  numerous  erosion  tests.  Thus 
a reliable  experimental  data  base  for  the  development  as  well  as  compari- 
son of  analytical  results  Is  rarely  found  In  the  solid  particle  litera- 
ture which  leaves  many  aspects  of  the  subject  unsettled.  The  analytical  I 

approaches  which  have  been  undertaken  to  describe  material  removal  due  to  I 

solid  particle  erosion  will  be  described.  I 
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3.0  EROSION  CORRELATIONS  AND  ANALYSES;  DUCTILF-  MATFBTAT.s 


A number  of  concepts  have  been  used  as  the  basis  for  the 
development  of  analytical  models  of  the  solid  particle  erosion  process. 

In  general  these  models  do  not  Incorporate  an  accurate  experimentally 
based  description  of  the  material  removal  process  Into  the  mathematical 
formulation.  The  mode  of  material  removal  Is  simply  a conjecture  or  Is 
not  specified  at  all  In  several  of  these  seml-emplrlcal  models.  How- 
ever within  the  last  few  years  experimental  studies  pertaining  to  ero- 
sion mechanisms  have  been  pursued  and  so  a better  physically-based  con- 
cept of  what  Is  taking  place  can  be  proposed.  These  studies  are  either 
clever  experiments  or  microscopic  observations  of  metallographlcally 
treated  surfaces  exposed  to  a small  number  of  particle  Impacts.  There- 
fore It  should  be  possible  within  the  next  few  years  to  develop  more 
physically-based  concepts  of  how  material  Is  removed  from  various  mate- 
rial surfaces  and  to  begin  to  develop  a new  generation  of  erosion  models. 

The  mathematical  development  for  many  of  the  existing  models  will 
be  given  In  a reasonably  unified  form.  The  historical  progression  of 
the  well-known  models  of  Flnnle,  Bitter,  Nellson  and  Gilchrist  Is  cited 
religiously  In  the  Introduction  to  numerous  solid  particle  erosion  papers. 
These  correlations  and  others  will  be  considered  here  also  but  In  a some- 
what more  critical  way.  However  a strictly  historical  approach  will  not 
be  followed.  Recent  results  on  the  mechanics  of  the  Interactions  which 
can  take  place  between  an  Impacting  particle  and  a solid  surface  will  be 
used  to  provide  perspective  for  the  real  significance  of  the  erosion 
models  or  theories. 

It  Is  Important  to  note  that  the  physical  basis  or  Initial  concept 
for  a particular  erosion  model  Is  based  on  a moderately  detailed  analysis 
of  a single  particle  Impact  and  the  material  removal  process  due  to  multi- 
ple Impacts  Is  estimated  by  simply  multiplying  the  material  removed  by 
a single  particle  by  the  total  nundier  of  particles  or  the  number  of 
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particles  multiplied  by  an  empirical  correction  (adjustment)  factor. 

The  primary  defect  in  this  approach  is  that  the  surface  being  eroded  by 
multiple  particle  collisions  is  quite  different  from  the  initial  surface 
for  which  the  single  particle  mass  loss  was  originally  derived.  The 
detailed  modeling  effort  thereby  degenerates  to  a semlemplrical  correla- 
tion. The  transition  from  Isolated  Impacts  to  a few  overlapping  Impacts 
has  been  Investigated  to  some  extent  but  certainly  requires  considerable 
attention  (Adler,  1974b,  1976b). 

The  general  categories  for  the  basic  elements  in  the  solid  particle 
erosion  analyses  are  the  dynamical  equations  of  motion,  the  Hertzian 
theory  of  Impact,  energy  concepts,  plasticity,  and  regression  analysis. 
Due  to  the  prominence  of  the  Hertzian  theory  in  these  modeling  efforts, 
the  Hertzian  relations  for  a spherical  body  Impacting  an  elastic  half- 
space are  provided  in  Appendix  A for  reference. 

3.1  Hutchings  and  Co-workers 

Hutchings  and  co-workers  (Hutchings  and  Winter,  1974,  1975;  Winter 
and  Hutchings,  1974;  Hutchings,  et  al.,  1976;  Hutchings,  1977,  1979) 
have  made  significant  advances  in  providing  Insights  into  the  material 
removal  process  due  to  both  single  spherical  and  angular  particle  im- 
pacts by  means  of  a series  of  clever  experiments.  The  impact  conditions 
are  idealized  and  the  particles  are  scaled  to  much  larger  dimensions 
than  those  characteristic  of  multiple  particle  erosion  for  high  speed 
photographic  analysis  of  the  dynamics  of  the  impacting  particle.  Each 
case — spherical  particles  and  angular  particles — will  be  considered  in 
turn.  In  both  cases  a compressed  gas  gun  with  a 16  mm  bore  is  used  to 
propel  the  particle.  The  spherical  balls  are  a hardened  steel  (with  a 
Vickers  hardness  of  9.0  GPa)  3 and  9.5  mm  in  diameter.  A few  different 
angular  particle  configurations  have  been  used  in  these  experiments 
(Winter  and  Hutchings,  1974;  Hutchings,  1975,  1979).  The  most  recent 
is  to  propel  square  plates,  8x8x1. 5 mm,  at  a stationary  Inclined  target. 
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Three  major  categories  of  oblique  particle  Impacts  on  ductile 
materials  have  been  experimentally  Identified  by  Hutchings  and  co- 
workers In  their  experimental  Investigations.  The  first  Is  a "ploughing" 
action  associated  with  a spherical  particle,  a rounded  surface  on  an  Ir- 
regular shaped  particle,  or  the  properly  oriented  flat  face  of  an  angular 
particle  as  shown  In  Figure  3.1  a and  b.  The  remaining  two  categories 
pertain  to  angular  particles  producing  cutting  deformations  (Figure  3.1c). 
Hutchings  and  co-workers  make  a clear  distinction  between  the  angle  of 
Impact,  a , and  the  orientation  of  the  particle  at  Impact  described  In 
terms  of  the  rake  angle  6 . The  rake  angle  Is  used  In  metal  cutting 
technology  and  Is  defined  In  Figure  3. Id.  The  angle  6 Is  negative  If 
It  Involves  a clockwise  rotation  between  the  normal  to  the  surface  and 
the  front  face  of  the  particle  and  positive  for  a counterclockwise  rota- 
tion. The  remaining  categories  In  Hutchings'  classification  are  "Type  I 
cutting"  when  a forward  rotation  Is  Imparted  to  the  angular  particle  and 
"Type  II  cutting"  when  a backward  rotation  occurs  after  the  particle 
strikes  the  target.  No  criterion  has  been  given  In  terms  of  impact 
velocity.  Impact  angle,  and  rake  angle  for  determining  iidilch  of  the  three 
categories  will  prevail.  A limited  number  of  transitions  can  be  found  In 
Hutchings  experimental  evaluations  and  computer  simulations  of  the  Impact 
process. 

Characteristic  crater  profiles  are  shown  In  Figure  3.2.  Only 
Type  II  cutting  results  In  material  removal  for  a single  Impact  due  to 
a mlcromachlnlng  action.  For  an  Impact  angle  of  30®,  Type  II  cutting 
only  occurs  at  rake  angles  ranging  from  0 to  -17";  Type  I cutting  takes 
place  for  rake  angles  ranging  from  -17®  to  -90®.  These  observations  are 
for  a square  plate.  As  Illustrated  In  Figure  3.1b  the  geometry  of  the 
Impacting  particle  will  also  contribute  to  the  form  of  the  deformation 
which  may  occur. 

Hutchings  (1979)  estimates  that  Type  I cutting  due  to  the  restricted 
range  of  rake  angles  over  which  It  occurs  would  only  be  observed  In  about 
one-sixth  of  the  Impacts  on  randomly  oriented  plates  at  a 30®  Impact 
angle.  For  Type  II  cutting  and  ploughing  deformations  additional  Impacts 
are  required  to  remove  the  material  displaced  from  the  crater. 
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Figure  3«1.  Deformation  Modes  Identified  by  Hutchings  and  Co-workers. 

a)  Flouring  Deformation  with  a Sphere,  b)  Ploughing 
Deformation  with  an  Angular  Particle,  c)  Cutting  Defor- 
mation with  an  Angular  Particle,  and  d)  Definition  of 
the  Rake  Angle,  6:  The  Rake  Angle  showr.  Is  Negative  by 
Convention. 
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3.1.1  Spherical  Body  Impacts  on  Ductile  Targets 

The  single  spherical  particle  Impact  experiments  were  for  the 
purpose  of  understanding  the  general  character  of  the  Indentation  In  a 
ductile  target  and  to  Identify  the  conditions  under  which  material  Is 
removed  from  the  surface.  A quantitative  study  was  made  of  the  depen- 
dence of  the  crater  volume  on  the  particle  Impact  angle  and  velocity. 

Since  the  Impact  craters  produced  by  particles  less  than  500  pm  In  diam- 
eter are  difficult  to  examine  In  detail  or  to  make  quantitative  measure- 
ments on  them,  Hutchings,  et  al.,  (1976)  elected  to  use  considerably 
larger  diameter  balls  In  their  experiments.  Hutchings,  et  al.,  found 
that  glass  spheres  larger  than  1 mm  in  diameter  tend  to  fragment  upon 
Impacting  mild  steel  targets  at  velocities  from  100  to  300  ms  Steel 
spheres  were  therefore  used. 

The  hypothesis  that  the  material  removal  mechanisms  for  large 
particles  are  representative  of  the  particulates  of  practical  Interest 
Introduced  by  Hutchings,  et  al.,  (1976)  Is  better  substantiated  than  In 
the  earlier  work  of  Sheldon  and  Kanhere  (1972).  Refer  to  Section  3.2. 

The  basis  for  the  comparison  made  by  the  authors  Is  Tabor's  analysis  of 
a rigid  sphere  Impacting  normally  on  a rigid  plastic  plane  (Tabor,  1951). 
Tabor  predicts  that  the  crater  volume  will  remain  constant  If  the  relation 

2 

density  x (velocity)  * constant  (3.1) 

Is  satisfied  for  spheres  of  the  same  diameter.  Hutchings,  et  al.,  make 
two  unsupported  extensions  of  this  result:  (1)  that  It  Is  applicable  to 
oblique  Impacts;  and  (2)  that  for  the  oblique  Impact  of  different  diam- 
eter spheres,  the  dimensions  of  the  crater  normalized  with  respect  to 
the  sphere  diameter  will  remain  constant  as  long  as  Eq.(3.1)  Is  satis- 
fied. These  conjectures  are  demonstrated  to  be  valid  based  on  a limited 
amount  of  experimental  evidence.  According  to  Hutchings,  et  al.  (1976), 
this  comparison  provides  justification  for  the  experimental  evaluations 
on  an  enlarged  scale  and  for  comparisons  with  the  solid  particle  erosion 
literature. 
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A series  of  Impact  craters  In  mild  steel  due  to  a 9.5  mm  steel 
ball  Impacting  over  a range  of  Impact  angles  and  velocities  was  evaluated. 
The  mass  of  metal  removed  from  the  target  was  measured  by  we  alng  the 
target  before  and  after  Impact.  The  volume  of  the  Impact  crater  was  de- 
termined by  machining  away  all  of  the  material  above  the  level  of  the 
undisturbed  metal  and  filling  the  depression  with  Plasticine  which  was 
then  removed  and  weighed.  Thus,  direct  measurements  are  made  of  the 
mass  loss  and  crater  volume  as  a function  of  Impact  angle  and  velocity. 
High  speed  photographic  records  were  used  to  evaluate  the  energy  lost 
during  the  Impact  process. 

A numerical  analysis  of  the  collision  process  was  also  developed 
In  conjunction  with  the  experimental  program.  The  computational  model 
will  be  described  and  then  the  results  of  Hutchings'  work  (Hutchings, 
et  al.,  1976)  will  be  discussed. 

The  temporal  progression  of  a sphere  entering  and  finally  leaving 
the  target  Is  determined  numerically  by  solving  the  equations  of  motion 
for  Incremental  time  steps  with  a computer.  Hutchings,  et  al.,  (1976) 
write  the  equations  of  motion  for  a rigid  sphere  Impacting  a ductile 
target  obliquely  as  shown  In  Figure  3.3.  They  assume  the  contact  sur- 
face between  the  sphere  and  the  deformed  material  Is  a circle  with  diam- 
eter AB.  The  resisting  force  acting  on  the  sphere  during  the  penetra- 
tion phase  Is  taken  to  be  proportional  to  the  coli^ct  area  determined 
by  AB  and  directed  through  the  center  of  the  sphere.  When  a frictional 
force  Is  also  allowed  to  act  on  the  contact  area,  the  equations  of  motion 


mx  ■ -N8ln(0+<j»)  - uNcos(B'H^) 
my  - -Ncos(0+<j))  + viNsln(&f4) 


(3.2) 


where  m Is  the  mass  of  the  Impacting  sphere,  I Is  the  moment  of 
Inertia  of  the  sphere,  and  the  remaining  quantities  are  defined  in 
Figure  3.3. 
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Figure  3.3.  Diagram  Showing  the  Region  of  Contact  and  the 

Forces  Assumed  to  Act  on  the  Sphere  In  the  Con- 
stant Indentation  Pressure  Model  (Hutchings, 
et  aL,  1976). 
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In  Hutchings'  equations,  Eq.(3.2),  the  mass,  radius.  Initial 
velocity,  and  Initial  angle  of  approach  of  the  sphere  are  prescribed 
quantities.  The  value  of  the  normal  force,  N , and  the  coefficient  of 
friction,  y , are  the  only  additional  parameters  which  have  to  be  speci- 
fied In  order  to  Integrate  the  equations  of  motion.  Based  on  experimen- 
tal data,  y was  assigned  a value  of  0.05.  The  force  N Is  determined 
from  the  relation 

N - p*A  (3.3) 

where  p Is  the  mean  Indentation  pressure  and  A Is  the  area  of  the 
contact  surface  (determination  of  the  contact  area  has  been  described, 
but  will  be  discussed  at  greater  length  shortly) . The  pressure  p Is 
assumed  to  be  constant  throughout  the  cratering  phase  and  Is  set  equal 
to  4.0  GPa  In  the  computations  reported  by  Hutchings,  et  al.  (1976). 

The  target  material  was  a work-hardened,  low-carbon  steel  with  a Vickers 
hardness  of  2.35  GPa.  The  magnitude  of  the  quaslstatlc  Indentation  pres- 
sure was  adjusted  to  account  for  the  high  strain  rates,  estimated  to  be  i 

on  the  order  of  10^  to  10^  sec  encountered  during  the  Impact  process 
(Hutchings,  1977). 

The  determination  of  the  contact  area,  A , In  Eq.(3.3),  requires  | 

further  consideration.  The  actual  contact  area  Is  proportional  to  the  ] 

area  of  a spherical  segment  Instead  of  the  approximation  of  a circle 
with  diameter  AB  Introduced  Into  the  analysis  by  the  authors.  In  pro- 
gressing through  the  penetration  process,  one  finds  the  difference  be- 
tween these  two  views  of  the  contact  area  can  approach  more  than  a factor 
of  two:  the  area  determined  by  AB  being  the  smaller  of  the  two.  The 
calculation  also  neglects  the  additional  contact  area  due  to  the  bulld-up 
of  deformed  material  around  the  Impact  crater:  especially  In  the 
direction  of  particle  motion. 
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The  significant  difference  In  the  magnitude  of  the  contact  area 
can  affect  the  magnitude  of  the  indentation  pressure  required  to  utch 
the  experimental  data.  It  should  of  course  be  realized  that  the  magni- 
tude of  the  pressure  p Is  not  going  to  be  uniform  over  the  contact 
area:  Its  value  will  be  high  In  the  central  portion  of  the  contact  zone 
and  It  will  decrease  toward  the  periphery.  The  contact  area  Introduced 
by  Hutchings,  et  al.,  would  probably  encompass  the  primary  zone  of  maxi- 
mum pressure.  On  the  other  hand,  the  model  assumes  a rigid-plastic  mate- 
rial so  that  It  is  necessary  for  any  material  In  contact  with  the  sphere 
to  be  displaced  from  Its  path  by  plastic  flow.  Hence,  on  the  basis  of 
the  material  response  Introduced  for  the  target  material,  the  pressure 
p should  be  reached  over  the  complete  contact  zone  in  order  to  maintain 
continuity  with  the  spherical  body.  When  the  kinetic  energy  of  the 
sphere  Is  no  longer  able  to  overcome  the  resisting  force  In  the  direc- 
tion normal  to  the  target's  surface,  the  sphere  begins  to  move  outward 
from  the  target  by  means  of  the  energy  In  the  tangential  direction. 

The  assumption  of  a rlgld-plastlc  target  and  rigid  projectile  Implies 
that  for  an  impact  normal  to  the  surface  of  the  target,  the  projectile 
would  penetrate  to  Its  maximum  depth  and  remain  embedded  at  this  depth 
In  view  of  the  absence  of  any  restoring  forces  acting  on  the  body.  This 
limiting  case  Illustrates  the  Importance  of  the  tangential  velocity  com- 
ponent In  the  model  and  the  necessity  for  elastic-plastic  response  If 
normal  Impacts  on  real  materials  are  to  be  described  as  considered  by 
Mamoun  (Section  3.7).  The  difference  between  the  calculated  and  experi- 
mental rebound  velocities  should  provide  a reasonable  estimate  of  when 
the  elastic  response  Is  significant.  For  this  reason,  Hutchings,  et 
al.,  did  not  carry  their  computations  beyond  Impact  angles  of  75*. 

The  experimental  data  for  the  mass  loss  and  crater  volume  as  a 
function  of  the  impact  velocity  and  angle  Is  shown  In  Figure  3.4  to 
3.7.  Excellent  agreement  Is  found  between  the  experimental  determina- 
tion of  the  energy  lost  In  the  collision  process  and  the  model  predic- 
tions. The  energy  lost  computation  Is  relatively  Insensitive  to  changes 


Figure  3.4.  Maas  Losses  for  9.5  am  Figure  3.5.  Crater  Volumes  for  9.5  mn 
Steel  Balls  Striking  Mild  Steel  Steel  Balls  Striking  Mild  Steel  Specl- 

Speclmens  at  30*  Impact  Angle.  The  mens  at  30*  Impact  Angle.  The  Solid 
Solid  Line  Represents  the  Relation  Line  Is  Theoretically  Derived  from 
H * 5 82x10**^^  V Hutching* s Model  and  has  a Slope  of 

o 2.4  (Hutchings,  et  al.,  1976). 

(Hutchings,  et  al.,  1976). 
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IMPACT  ANQLE/DEGREES 


Figure  3.6 


The  Variation  of  Mass  Loss  with  Impact  Angle  for  9.5  nm 
Steel  Balls  Striking  Mild  Steel  at  270  ms-l.  The  Broken 
Line  Represents  the  Results  of  Klels  (1969)  for  Multiple 
Impacts  of  Round  Iron  Shot  on  Mild  Steel  Scaled  to  270  mi 
Impact  Velocity  (Hutchings,  et  al.,  1976) 


CRATER  VOLUM^^miS 


Figure  3 


IMPACT  ANQLE/DEQREES 


7,  The  Variation  of  Crater  Volume  with  Impact  Angle  for  9.5  mm 
Steel  Balls  Striking  Mild  Steel  at  270  ms-l.  The  Solid 
Curves  are  Theoretically  Derived  from  the  Model,  Assuming 
the  Plastic  Indentation  Pressure  to  be  2.35  GPa  and  3.50  GPa 
Note  the  difference  in  form  compared  with  the  experimental 
data:  dashed  curve  (from  Hutchings,  et  al.,  1976). 


in  Che  value  of  Che  plasClc  IndenCaClon  pressure  p : che  experlmenCal 
daCa  Is  flC  qulCe  well  when  p > 4.0  GPa  . On  Che  oCher  hand,  Che  model 
predlcClons  in  Figure  3.7  do  noc  reflecC  Che  general  Crend  In  Che  daCa. 
The  macch  appears  Co  be  reasonable  for  a^SO”  , however  for  a>30'* 
sizable  devlaClons  In  Che  predlcCed  and  experlmenCal  craCer  volumes 
occur.  The  craCer  volume  as  a funcClon  of  ImpacC  veloclCy  shown  In 
Figure  3.5  when  a*20’‘  happens  Co  be  a condlClon  which  shows  excellenC 
agreemenc  beCween  experlmenC  and  cheory  when  p»4.0  GPa.  HuCchlngs,  ec 
al.  (1976),  aCCempC  Co  flc  Che  experlmenCal  craCer  volume  measuremencs 
In  Figure  3.7  by  varying  Che  value  of  p , however  Che  general  Crend  In 
Che  daCa  does  noC  appear  Co  be  reproduced  as  long  as  p remains  con- 
sCanC.  They  IndlcaCe  ChaC  Che  craCer  volumes  are  inversely  proporClonal 
Co  p . They  sCaCe  ChaC 

”IC  appears  from  Chese  resulCs  ChaC  Che  value  of  p required 
Co  give  Che  besc  flC  wlch  Che  experlmenCal  daCa  decreases  ap- 
preciably wlch  Increasing  ImpacC  angle.  Ic  may  be  argued  ChaC 
Chls  Is  conslscenc  wich  Che  deformaClon  paCCem  for  dlfferenC 
angles  of  ImpacC;  aC  low  ImpacC  angles  Che  deformaClon  becomes 
more  localized  In  Che  surface  and  Che  overall  shear  s Crain  and 
hence  sCraln  race  become  higher." 

The  model  calculaClons  do  noC  apply  Co  Che  mass  loss,  since  chere 
Is  no  mechanism  In  Che  model  for  dlsclngulshlng  beCween  Che  maCerlal 
displaced  and  Che  maCerlal  which  Is  acCually  removed.  The  naCure  of 
Che  maCerlal  removal  process  IdenClfled  by  HuCchlngs  and  co-workers  will 
now  be  described  for  rounded  body  ImpacCs  based  on  chelr  Idealized 
experlmenCal  condlclons. 

The  major  ImpeCus  for  HuCchlngs  Co  use  large  parclcles  was  Co 
beCCer  underscand  Che  maCerlal  removal  process.  The  value  of  Chls  under- 
sCandlng  Is  Co  Idenclfy  Che  way  In  which  Che  properCles  of  a maCerlal 
conCrol  iCs  erosive  response.  For  essenClally  rigid  spheres  lapacClng 
mild  sCeel  CargeCs  obliquely,  meCal  is  displaced  from  Che  resulclng 
craCer  by:  (1)  acconmodaClon  of  che  displaced  maCerlal  Chrough  bulk 


deformation  of  the  material  surrounding  the  crater,  and  (2)  creation  of 
a shear  lip  at  the  end  of  the  crater  where  the  particle  leaves  the  sur- 
face. The  voluM  of  material  comprising  the  shear  lip  ranges  from  1/12 
to  1/4  of  the  total  displaced  volume.  Hutchings'  experiments  Indicate 
that  a critical  velocity  exists  for  each  Impact  angle  below  which  the 
lip  remains  attached  and  above  which  It  Is  removed.  A band  of  Intense 
shear  deformation  Is  observed  at  the  base  of  the  lip.  Separation  of 
the  lip,  when  It  takes  place,  occurs  along  this  shear  band. 

Sheldon  and  Kanhere  (1972)  made  an  essentially  similar  observation 
for  the  material  removal  mechanism  for  20  and  90  degree  Impacts  on  alumi- 
num targets  using  2.3  mm  steel  spheres  and  3 mm  glass  beads.  These  in- 
vestigators observe  for  a 20®  Impact  angle  and  velocity  of  232  ms  ^ 

(760  fps)  that  the  surface  material  flows  In  front  and  to  the  sides  of 
the  advancing  particle  until  It  Is  strained  sufficiently  to  fracture 
resulting  In  material  loss.  They  note  that  material  builds  up  on  the 
forward  side  of  the  crater,  but  the  volume  of  this  material  is  much 
less  than  the  total  displaced  volume.  The  forward  lip  appears  to  frac- 
ture sooner  for  the  work-hardened  6061  aluminum  compared  with  specimens 
In  the  annealed  state  and  most  of  the  displaced  material  fractures 
leaving  a cavity  with  no  extruded  material  at  the  periphery.  The  greater 
susceptibility  to  fracture  of  the  work-hardened  surface  appears  In  the 
higher  mass  loss  recorded  for  this  condition  over  tb  annealed  state. 

The  experimental  data  shows  that  there  Is  an  Impact  velocity  below  which 
the  work-hardened  surface  may  reduce  the  mass  loss.  It  may  be  that  the 
mass  loss  differences  are  related  to  the  strain  rates  which  increase 
with  Impact  velocity. 

Sheldon  and  Kanhere  find  the  same  material  removal  mechanism 
applies  to  90®  Impacts  at  232  ms  They  speculate  that  material  now 
flows  evenly  around  the  periphery  of  crater  due  to  the  Impact  by  glass 
and  steel  spheres  and  observe  that  a relatively  small  ridge  exists  at 
the  periphery  of  the  crater.  The  details  of  the  material  removal  process 
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are  not  as  well  documented  as  In  the  work  of  Hutchings,  et  al.,  1976. 

From  Sheldon  and  Kanhere's  mass  loss  measurements  It  Is  seen  that  there 
is  relatively  little  variation  In  the  absolute  magnitude  of  the  mass 
loss  between  the  20*  and  90*  glass  bead  Impacts  over  the  velocity  range 
Investigated.  A wider  variation  between  the  mass  loss  for  these  two  Im- 
pact conditions  Is  found  for  steel  shot.  It  would  be  of  Interest  to 
differentiate  between  mass  loss  and  crater  volume  In  the  manner  developed 
by  Hutchings  to  try  to  better  understand  the  material  removal  process  In 
aluminum  and  to  examine  a broader  range  of  Impact  angles. 

It  Is  significant  to  note  that  according  to  Hutchings,  et  al., 
only  a portion  of  the  displaced  volume  Is  removed  during  a single  colli- 
sion. The  work  of  Sheldon  and  Kanhere  would  Indicate  a significantly 

larger  portion  of  the  displaced  volume  Is  removed  from  6061  AJl.  Most 
analyses  for  erosion  rates  develop  expressions  for  the  volume  of  mate- 
rial displaced,  but  then  employ  heuristic  arguments  to  make  the  material 
removal  rates  conform  to  the  experimentally  determined  values. 

3.1.2  Angular  Plate  Impacts  on  Ductile  Targets 

The  Investigation  of  the  Impact  of  angular  particles  on  ductile 
targets  undertaken  by  Hutchings  and  Winter  (Winter  and  Hutchings,  1974; 
Hutchings,  1977,  1979)  parallels  the  Investigation  of  spherical  body 
collisions.  The  first  point  to  be  considered  Is  the  scaling  of  the  Im- 
pact damage  due  to  small  particulates  on  the  order  of  100  urn  to  that 
produced  by  steel  plates  on  the  order  of  10  mm.  This  point  Is  not  ade- 
quately substantiated.  The  Initial  comparison  (Hutchings  and  Winter, 
1974)  Is  for  spherical  bodies  and  this  comparison  Is  simply  adopted  for 
angular  bodies  (Hutchings,  et  al.,  1976).  Winter  and  Hutchings  (1974) 
note  that  It  Is  found  experimentally  that  the  erosion  of  metals  by  an- 
gular particles  greater  than  100  Urn  Is  Independent  of  particle  size 
(Tilly,  1969).  However  this  leaves  open  the  question  as  to  whether  the 
mode  of  material  removal  for  a single  particle  Impact  Is  identical  to 


i i' 

; t 

i i 


I 

\ 


I 

t 


1 ' 
1 ; 


\ 


1 


that  due  to  multiple  particle  impacts  on  a completely  obscured  surface. 
These  Investigators  also  refer  to  the  experimental  Justification  for 
scaling  presented  by  Sheldon  and  Kanhere  (1972).  However,  once  again, 
the  work  of  Sheldon  and  Kanhere  pertains  to  spherical  bodies  and,  as  will 
be  shown  in  Section  3.2,  the  basis  for  their  verification  is  question- 
able. Thus,  adequate  justification  cannot  be  found  in  Hutching's  pub- 
lications for  the  exactness  of  the  representation  of  the  mode  of  defor- 
mation associated  with  angular-shaped  erosive  particulates  by  geometri- 
cally idealized  particles  with  dimensions  of  several  millimeters.  On 
the  other  hand,  the  qualitative  features  of  the  damage  produced  under 
the  idealized  experimental  conditions  Hutchings  employs  are  quite  in- 
structive in  understanding  the  events  which  can  take  place  when  an 
angular  particle  Impacts  a ductile  target. 

The  effect  of  particle  orientation  can  be  expressed  in  terms  of 
the  rake  angle  as  Indicated  in  Figure  3.1.  However,  in  addition  to  the 
Impact  angle,  a , and  the  rake  angle,  6 , the  location  of  the  center 
of  mass  of  the  particle  is  also  a factor  in  determining  which  of  the 
three  classifications  (ploughing.  Type  I cutting,  and  Type  II  cutting) 
identified  by  Hutchings  (1979)  will  occur  (Winter  and  Hutchings,  1974). 
The  effect  of  the  shape  of  the  particle  (distribution  of  mass)  on  the 
material  removal  process  is  Illustrated  in  Figure  3.1.  The  influence 
of  Che  frangiblllty  of  the  impacting  particle  was  also  described  by 
Winter  and  Hutchings  (1974). 

In  subsequent  work  (Hutchings,  1977),  the  particle  geometry  was 
restricted  Co  a square  plate  which  removed  the  influence  of  particle 
shape,  other  than  the  angularity  of  Che  particle,  on  the  material  re- 
moval process.  The  plate  was  propelled  from  a gas  gun  such  that  a 
comer  of  the  plate  Impacted  the  inclined  plate  at  a predetermined  rake 
angle.  The  relevant  two-dimensional  analysis  is  analogous  to  that  al- 
ready described  for  spherical  body  impacts.  The  plate  is  assumed  to  be 
rigid  impacting  upon  a rigid-plastic  target.  The  equations  of  motion 


for  this  condition  are  then  formulated  and  integrated  incrementally  with 
respect  to  time  to  calculate  the  position  and  orientation  of  the  plate 
while  in  contact  with  the  target.  The  time-dependent  contact  area  is 
computed  and  a uniform  indentation  pressure  is  assumed  to  act  on  this 
area.  Frictional  forces  are  neglected.  Several  possibilities  are 
likely  to  occur  for  the  square  plate  Impact  as  depicted  in  Figure  3.8. 
Determination  of  the  contact  areas  for  the  conditions  in  Figure  3.8  b, 
c,  and  d become  a part  of  the  problem.  The  equations  of  motion  cannot 
be  written  explicitly,  but  depend  on  the  temporal  evolution  of  both  the 
translation  and  rotation  of  the  plate  while  in  contact  with  the  target. 
Hutchings  (1977)  divides  the  general  problem  into  four  cases  represent- 
ing the  current  location  of  the  contact  zone  and  its  new  location  due 
to  an  incremental  rotation.  The  computational  scheme  reformulates  the 
governing  differentiation  equations  after  each  time  step.  The  volume 
displaced  during  the  complete  Impact  sequence  can  be  computed  as  a func- 
tion of  the  initial  Impact  angle,  impact  velocity,  and  rake  angle. 

Energy  partitioning  before  and  after  Impact  is  readily  evaluated. 

Examples  of  the  computational  results  can  be  found  in  Hutchings  (1977). 

The  solid  particle  erosion  studies  undertaken  by  Hutchings  and 
co-workers  at  the  University  of  Cambridge  over  the  last  five  years  have 
made  a significant  contribution  to  understanding  the  modes  of  material 
removal  for  angular  particles  impacting  ductile  targets  obliquely.  Prior 
analyses  have  been  quite  conjectural  as  to  the  nature  of  the  particle/ 
target  Interactions.  The  computer  studies  in  conjunction  with  micro- 
graphic  studies  of  single  and  multiple  particle  crater  configurations 
associated  with  carefully  graded  particulates  should  provide  much  addi- 
tional valuable  information  on  the  true  volume  removal  rates  for  a 
variety  of  impact  conditions. 
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MNECTION 
OF  MOTION 


FINAL  CONFIGURATION 


(a)  TRANSLATION  WITH 
LITTLE  ROTATION 
FRONT  CONTACT  ONLY  (AB) 


(b)  TRANSLATION  AND 
FORWARD  ROTATION 
FRONT  AND  REAR  CONTACT 
(AB  AND  CD) 


(c)  TRANSLATION  AND 
BACKWARD  ROTATION 
FRONT  AND  REAR  CONTACT 
(AB  AND  CD) 


(d)  TRANSLATION  AND 
FORWARD  ROTATION 
FRONT  CONTACT  ONLY  (AB) 


(a)  TRANSLATION  WITH 
NO  ROTATION 
FULL  FRONT  AND  REAR 
CONTACT  (AB  AND  CD) 
(RESULTANT  VERTICAL) 


Figure  3.8.  Dynamics  of  Square  Plate  Impacts  on  a 
Ductile  Target  (Hutchings,  1978). 
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3.2  Sheldon  and  Kanhere 


Sheldon  and  Kanhere  (1972)  preceded  Hutchings  and  his  associates 
In  elevating  the  scale  of  the  particle  Impact  and  removal  process  In 
order  to  understand  the  erosion  process  for  ductile  metals  by  Investi- 
gating the  effect  of  large  single  particle  Impacts  on  an  aluminum  alloy. 
In  the  same  paper  they  presented  an  analysis  based  on  the  volume  of  mate- 
rial displaced  In  a single  spherical  Indent  to  obtain  a velocity  expo- 
nent of  3 Instead  of  2 which  was  Inherent  In  the  available  analyses  but 
which  was  not  representative  of  the  experimentally-determined  values  as 
described  In  Section  2.  The  main  features  of  both  the  experimental  and 
analytical  parts  of  Sheldon  and  Kanhere's  work  will  be  reviewed.  The 
primary  results  cannot  be  substantiated. 

3.2.1  Single  Sphere  Impact  Experiments  for  Ductile  Targets 

Sheldon  and  Kanhere  (1972)  support  their  contention  that  steady- 
state  material  removal  can  be  determined  from  single  spherical  bead  Im- 
pacts on  ductile  target  materials  on  the  basis  of  two  observations. 

The  first  Is  that  the  normalized  weight  loss  measurements  for  a single 
3 mm  diameter  glass  bead  Impacting  a 6061-T0  aluminum  alloy  Is  "approxi- 
mately equivalent"  to  the  measured  erosion  rates  quoted  from  the  litera- 
ture for  multiple  glass  bead  Impacts  on  aluminum  alloys.  The  second 
observation  Is  that  the  erosion  characteristics  are  quite  similar  when 
either  spherical  or  very  angular  particles  are  used.  The  erosion  char- 
acteristics they  Identify  as  being  comparable  are  the  variation  of  mass 
loss  with  Impact  angle,  the  dependence  on  particle  velocity  and  size, 
and  the  rippled  appearance  of  the  eroded  surface.  These  latter  obser- 
vations are  the  authors'  personal  views  In  that  they  are  simply  stated 
without  any  supporting  evidence  to  convince  the  reader  that  they  are  In- 
deed valid  comparisons.  The  authors  conclude  on  the  basis  of  the  above 
that  the  same  material  removal  mechanisms  are  operative  when  either  an- 
gular or  rounded  particles  are  the  eroding  medium.  This  conclusion 
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considerably  broadens  the  relevance  of  their  single  particle  Impact 
studies.  The  basis  for  this  conclusion  will  be  considered  In  terms 
of  the  data  Sheldon  and  Kanhere  present  In  their  paper. 
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In  support  of  the  equivalence  between  the  material  removed  by  a j 

single  glass  bead  Impact  and  multiple  glass  bead  Impacts,  Sheldon  and  | 

Kanhere  provide  the  data  In  Table  3.1,  which  has  been  augmented  by  addl-  | 

tlonal  values  which  can  be  obtained  from  the  references  cited.  Our  In-  I 

terpretatlon  of  this  data  Is  provided  In  the  lower  half  of  Table  3.1. 

It  Is  not  stated  whether  pure  aluminum  or  an  aluminum  alloy  was  used  In 

Nielson  and  Gilchrist’s  erosion  tests  (Nielson  and  Gilchrist,  1968).  j 

It  Is  thought  that  Ail  2618  was  used  In  the  work  of  Tilly  and  Sage  (1970) 

based  on  the  composition  reported  In  their  paper.  The  experimental  data  | 

from  these  two  papers  Is  shown  In  Figure  3.9.  The  comparison  between  ^ 

the  rounded  and  angular  particle  data  Is  shown  In  Figure  3.9b  for  Impact  | 

angles  from  0 to  90*.  This  direct  comparison  from  Tilly  and  Sage  (1970)  i 

shows  a decisive  difference  between  the  two  conditions  which  In  Itself  I 

i 

weakens  the  relationships  Sheldon  and  Kanhere  are  trying  to  establish.  ^ 

One  also  notes  a significant  change  In  the  shape  of  the  experimental  I 

curves  reported  by  Nielson  and  Gilchrist  (Figure  3.9a)  and  Tilly  and 
Sage  (Figure  3.9b)  for  glass  beads.  Furthermore,  there  appears  to  be 
an  Inconsistency  In  the  data  reported  by  Tilly  and  Sage  In  that  the 
values  of  the  mass  loss  for  90°  Impact  data  for  glass  beads  does  not 
agree  with  the  values  read  from  plots  of  the  velocity  dependence  of  the 
mass  removal  rates  found  In  the  same  paper.  The  mass  loss  read  from  the 
plot  of  the  velcotly  dependence  of  125-130  ym  glass  beads  Impacting  nor- 
mally onto  aluminum  targets  Is  shown  In  parentheses  In  the  last  line  of 
Table  3.1.  The  relevant  multiple  glass  bead  erosion  data  does  not 
appear  to  be  very  consistent  In  toto. 

Sheldon  and  Kanhere  state  that  the  mass  loss  due  to  single  particle 
Impacts  was  determined  by  accurately  weighing  the  test  specimen  before 
and  after  Impact.  Three  particles  mounted  In  a styrofoam  sabot  were 
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I TABLE  3.1.  Comparison  of  Single  and  Multiple  Glass  Particle 

, Erosion  Tests  on  Aluminum. 


PARTICLE 

DATA 

EROSION 

(mg/g) 

VELOCITY 
fps (ms-l) 

IMPACT  ANGLE 
(degrees) 

MULTIPLE  IMPACTS 

SINGLE  IMPACTS 

260  ( 79) 

20 

0.36  * 

0.12 

800  (244) 

20 

** 

0.40 

1.8 

800  (244) 

90 

2.5 

1.8 

260  ( 79) 

20 

* 

0.35 

0.11 

800  (244) 

20 

** 

.36 

1.7 

260  ( 79) 

90 

.073  * 

0.14 

800  (244) 

90 

2.4  (5.9)* 

1.9 

* 

Nellson  and  Gilchrist  (1968). 
Tilly  and  Sage  (1970). 
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EROSION  (mg/g)  EROSION  (mg/g) 


I 


IMPACT  ANGLE  (DEGREES) 


Figure  3.9.  Impact  Angle  Dependence  for  Erosion  of  Aluminum 
Alloys,  a)  475  um  Glass  Spheres  Impacting  at 
79  ms~^  (Mellson  and  Gilchrist,  1968);  b)  Impact 
Velocity:  244  ms”!  (Tilly  and  Sage,  1970). 


fired  from  a gas  gun  for  each  test.  The  average  from  four  shots 
(totaling  12  Impacts)  constitute  the  experimental  value  for  each  test 
condition  they  record  In  their  paper.  Accepting  for  the  moment  this 
method  for  the  mass  loss  evaluation,  the  comparison  between  the  authors' 
single  particle  mass  loss  evaluations  and  the  erosion  data  from  other 
Investigations  Is  not  very  convincing.  Unfortunately  the  two  sets  of 
erosion  data  [Nielson  and  Gilchrist  (1968)  and  Tilly  and  Sage  (1970)] 
show  distinctly  different  locations  for  the  maximum  erosion  rates  with 
respect  to  the  Impact  angle  (refer  to  Figure  3.9),  so  little  can  be  said 
regarding  Che  trend  In  the  single  bead  Impact  data  except  that  It  ap- 
pears to  show  relatively  little  variation  with  Impact  angle: 
considerably  less  chan  either  multiple  Impact  evaluation. 


3.2.2  Analysis  of  Material  Removal  for  Normal  Impacts  of  Spherical 
Beads  on  Ductile  Targets 


Sheldon  and  Kanhere  provide  an  analysis,  based  on  the  indentation 
hardness  of  the  target  material,  to  establish  that  Che  material  removal 
per  gram  of  uniformly  sized  particles  Is  proportional  to  Che  cube  of 
Che  Impact  velocity  and  not  to  the  velocity  squared  as  typically  found 
from  energy  considerations.  However,  It  will  be  shown  the  cubic  velocity 
dependence  Is  not  supported  by  their  own  analysis.  Without  going  through 
Che  details  of  the  derivation,  Sheldon  and  Kanhere  show  that  [Sheldon 
and  Kanhere  (1972),  Eq.(6)] 


1/2  f 1 9 1 

,(V»v)  -l-23[i(3X-X^)J 


(3.4) 


where  Is  Che  particle  Impact  velocity 

Is  the  density  of  the  Impacting  particle 
Is  the  Vickers  hardness  of  the  target  material 
X - (X<1) 

q Is  Che  depth  of  the  Impact  crater 
D Is  Che  diameter  of  the  Impacting  sphere. 
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This  relation  Is  obtained  from  an  energy  balance  between  the  kinetic 
energy  of  a sphere  Impacting  normal  to  the  surface  of  the  target  and 
the  work  done  In  creating  the  final  crater  dimensions  In  accordance 
with  Meyer's  hardness  relation.  For  X<p.3  , the  final  crater  depth  Is 
proportional  to  the  Impact  velocity.  For  0.3<X<1.0  , the  term  Invol- 

3 

vlng  X becomes  more  significant  as  X Increases  to  Introduce  some 
nonlinearity  into  the  functional  dependence  which  is  not  indi- 

cated on  the  authors'  plot  of  Eq.(3.4).  On  the  basis  of  Eq.(3.4)  and 
the  experimental  results  for  single  bead  Impacts  on  6061-T0  aluminum, 
Sheldon  and  Kanhere  conclude  that 

q “ (3.5) 

This  result  does  not  appear  unreasonable  although  one  wonders  why  the 
experimental  data  for  20°  Impacts  appears  to  cluster  around  the  analy- 
tical values  for  90°  impacts  (Curve  I in  Figure  3.10),  while  the  trend 
In  the  crater  depth  measurements  from  90°  Impact  experiments  (Curve  II 
in  Figure  3.10)  are  displaced  from  the  computed  depth-velocity  relation. 
It  Is  to  be  noted  that  the  experimental  data  Is  quite  sparse.  In 
Figure  3.10  a straight  line  Is  drawn  through  data  points  representing 
four  different  Impact  conditions.  It  Is  seen  that  If  taken  on  an  Indi- 
vidual basis,  each  impact  condition  shows  a distinctly  different  trend 
than  that  for  the  consolidated  data. 

Sheldon  and  Kanhere  make  the  perplexing  statement  that  "since 
the  dimensions  of  the  crater  formed  by  Impact  are  all  proportional  to 
q , and  since  the  amount  of  material  removed  Is  nearly  the  full  crater 
size,  then.  It  Is  to  be  expected  that  material  removal  per  particle," 

W , would  be  proportional  to  the  cube  of  the  maximum  penetration  depth, 

q > 

W « q^  (3.6) 
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Figure  3.10.  Conparlson  of  the  Hardness  Penetration  Theory 
(Curve  I)  with  Experimental  Erosion  Results. 

Note  discrepancy  between  theoretical  prediction 
and  experimental  data  for  normal  Impacts,  (from 
Sheldon  and  Kanhere,  1972). 
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So,  according  to  Eq.(3.S),  W “ V , The  volume  of  the  Indent  is 
actually  more  strongly  proportional  to  the  square  of  the  penetration 
depth  as  shown  in  Eq.(B4)  of  Appendix  B. 


Returning  to  Eq.(3.4),  the  term  in  brackets  is  identified  as 
being  proportional  to  the  volume  of  an  Incomplete  spherical  cap,  v)  , 
referring  to  Eq.(BA)  in  Appendix  B 


/D  2 1 3\ 

VJ  - q - - q j 


(3.7 


which  can  be  rewritten 


u - 2|-[i  (3x'  - X^] 


(3.8 


Substituting  into  Eq.(3.4), 

3 


(1.23)  8 'h  ' 


(3.9 


The  corresponding  expression  obtained  by  Sheldon  and  Kanhere  la 

^ 3/2 


U « q' 


(1.23)^  8 ° 


(3.11 


H 


[with  minor  corrections].  However,  this  result  is  based  on  a highly 
erroneous  measure  of  the  crater  volume.  A similar  observation  was  made 
by  Hutchings  in  his  Ph.D.  thesis  (Hutchings,  1978). 

Eq.(3.9)  is  a consistent  result  for  the  model  used  by  Sheldon  and 
Kanhere  as  the  basis  for  their  analysis:  the  material  removal  rate  re- 
mains proportional  to  the  square  of  the  impact  velocity.  The  incorrect 
assumption  made  in  Eq.(3.6)  is  the  reason  for  the  velocity  cubed 


dependency  which  the  authors  state  as  a principal  conclusion  from  their 
analysis.  According  to  Eq.(3.9)y  the  Vickers  hardness  Is  once  again 
inversely  proportional  to  the  rate  of  material  removal  consistent  with 
Flnnie's  experimental  data  (Flnnle,  Wolak,  and  Kabll,  1967)  and  not  to 
the  (-3/2)  power  which  the  authors  cannot  Justify  on  the  basis  of  the 
available  experiment  data. 


In  summary,  Sheldon  and  Kanhere's  work  does  not  develop  a valid 
experimental  basis  for  establishing  a correspondence  between  single 
spherical  particle  Impacts  and  the  material  removal  rates  for  multiple 
Irregular  shaped  particles.  In  addition,  the  paper  does  not  provide  a 
sound  analytical  derivation  for  demonstrating  that  material  removal 
rates  proportional  to  the  cube  of  the  particle  Impact  velocity  can  be 
achieved.  Therefore,  both  the  experimental  and  analytical  results 
contained  In  this  paper  are  of  questionable  value. 


3.3  Flnnle 


Flnnle  provided  a simple  analysis  of  the  erosion  of  ductile  metals 
due  to  angular  particle  Impacting  at  oblique  angles  of  attack  (Flnnle, 
1958;  1960) . The  equations  obtained  In  1958  have  remained  essentially 
unchanged  to  the  present  time.  The  general  approach  will  be  described. 
The  derivation  of  the  basic  equations  was  reformulated  In  1965  (Flnnle 
and  Kabll,  1965)  and  summarized  with  the  assumptions  necessary  for  their 
derivation  In  1972  (Flnnle,  1972).  A modification  of  the  original  equa- 
tions was  proposed  in  1977  which  provided  a variable  velocity  exponent 
which  was  greater  than  2 (Flnnle  and  McFadden,  1978).  These  equations 
have  been  reviewed  recently  by  Flnnle  (Flnnle,  Levy,  and  McFadden,  1979) 
Flnnie's  model  and  associated  analysis  have  been  widely  referenced  In 
the  solid  particle  erosion  literature  over  the  intervening  two  decades. 
Relatively  few  new  conceptual  ondels  or  detailed  particle  Impact  experi- 
ments have  provided  better  Insights  Into  the  material  removal  process 
during  this  period,  except  for  the  Innovative  studies  of  Hutchings  and 
his  colleagues. 


The  erosion  of  ductile  metals  was  described  by  Flnnle  In  terms  of 
the  volume  removed  by  a single  angular  particle  Impacting  obliquely. 

The  volume  removed  Is  determined  by  writing  the  equations  of  motion  for 
a single  particle  Interacting  with  the  surface  and  then  estimating  the 
volume  removed  from  the  particle's  trajectory  as  It  cuts  through  the 
surface  layer  of  the  target  material.  The  analysis  assumes  that  the 
target  material  Is  rlgld-plastlc  and  that  material  Is  removed  due  en- 
tirely to  the  displacing  or  cutting  action  of  the  Impacting  particle. 
Flnnle's  analysis  Is  a special  case  of  Hutchings’  analysis  for  angular 
particles.  The  conceptual  model  upon  which  Flnnle's  analysis  Is  based 
Is  Type  I cutting  as  defined  by  Hutchings  (1979).  Hutchings'  results 
would  Indicate  that  Flnnle's  cutting  mechanism  Is  a single  particle 
material  removal  mechanism,  but  the  conditions  necessary  for  Its  opera- 
tion In  a random  distribution  of  Impacting  particles  may  occur  for  only 
a small  number  of  particle  impacts. 

Flnnle's  concept  of  an  angular-shaped  particle  Impacting  a 
ductile  surface  at  an  acute  angle  Is  shown  In  Figure  3.11.  The  dynamical 
equations  governing  the  motion  of  the  particle  after  It  strikes  the 
surface  are: 

mx  + F * 0 (3.11) 

C X 

my^  + KF^  - 0 (3.12) 

1$  + F^r  - 0 (3.13) 

The  equations  are  written  for  motion  relative  to  the  center  of  mass  of 
the  particle.  The  particle  has  a uniform  width,  b . The  horizontal 
resisting  force  of  the  ductile  target  material,  F^  , Is  equal  to  the 
flow  stress,  0^  , for  the  ductile  target  times  the  area  of  the  cutting 
face  of  the  particle  In  contact  with  the  material  as  the  particle  pene- 
trates the  target.  The  area  of  the  contact  zone  normal  to  the  horizontal 
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force  component  Is,  in  accordance  with  Figure  3.11, 


( 

I 


i 

I 


I 


1 


A - f,b 


(3.14) 


so 


p'I'by. 


(3.15) 


Flnnle  assumes  that  during  the  cutting  process  the  ratio  Is 

a constant  and  that  (Flnnle,  1958).  The  latter  relatJ^on  re- 

quires that  the  particle  be  rigid  and  Its  rotation  during  the  cutting 
process  Is  negligible.  Flnnle  further  assumed  that  the  ratio  of  the 
vertical  force  component,  F , to  the  horizontal  force  component,  F , 

y ^ 

Is  a constant  K . Arguments  based  on  metal  cutting  are  presented  to 
provide  numerical  estimates  for  the  constants  and  K (Flnnle,  1958, 
1960) . 


It  Is  seen  that  the  model  Is  essentially  for  a rigid  cutting  tool 
removing  a chip  of  material.  This  correspondence  Is  further  emphasized 
in  that  the  particle  rotation  during  the  period  of  contact  Is  taken  to 
be  small.  It  Is  now  known  from  the  work  of  Hutchings  (Hutchings,  1977, 
1979)  that  this  is  a somewhat  restrictive  assumption.  In  Flnnle's  model 
the  entire  forward  face  of  the  particle  which  penetrates  the  target  de- 
fines the  contact  area.  This  condition  will  not  prevail  If  particle 
rotations  are  permitted  as  shown  In  Figure  3.8.  The  assumption  of  a 
constant  ratio,  K , between  the  vertical  to  horizontal  force  fixes  the 
direction  of  the  force  vector  acting  on  the  contact  area.  Thus,  there 
Is  no  provision  In  Flnnle's  analysis  for  variations  In  the  orientation 
of  the  particle  prior  to  Impact,  These  two  primary  assumptions  enable 
Flnnle  to  solve  the  governing  differential  equations  analytically. 


The  volume  removed  by  the  cutting  action  of  the  particle  Is  simply 


^ ■ b j y^dx^ 


(3.16) 


as  Illustrated  In  Figure  3.11.  There  are  two  conditions  which  can  be 
realized:  the  particle's  kinetic  energy  Is  dissipated  before  the  chip 
Is  removed  so  the  particle  motion  terminates  while  It  Is  still  In  con- 
tact with  the  surface,  or  the  particle  completes  the  cut  and  leaves  the 
surface.  These  conditions  will  be  used  In  conjunction  with  Eq.(3.16) 
upon  completion  of  the  Integration. 

In  Flnnle  (1958)  the  equations  of  motion  were  Integrated  and 
explicit  expressions  obtained  for  x^.  and  y^  which  were  then  used  to 
evaluate  Eq.(3.16).  A slightly  more  direct  approach  was  used  In  Flnnle 
and  Kabll  (1965)  and  utilized  In  the  continued  discussion  of  the  model 
presented  In  Flnnle  (1972).  The  latter  derivation  will  be  given  In 
moderate  detail  to  Indicate  the  basic  elements  of  Flnnle's  theory. 


For  a rigid  particle  In  Flnnle's  analysis 


x^  “ X + r$ 
t c 


Substituting  Eq.(3.11)  and  (3.13)  Into  Eq.(3.17)  yields 


(3.17) 


-Ob  * 


(3.18) 


where  F^  Is  given  explicit  form  from  Eq.(3.15).  Substituting  Eq.(3.18) 
Into  Eq.(3.16), 


b r ^ . 

O - - -p  J Xj.dXj. 


(3.19) 


I 


.iwipiwiqpiljjff  lyipn 


where  is  the  value  of  for  which  the  cutting  action  terminates 


at  a time  t . Note  that 
c 


/"  xdx  -/  X X.  t - i / ‘=^(x) 


1 r c d ,..2.^  1 . 2 

2 J “ 2 * 

o 


The  initial  condition 


V cos  a 
o 


(3.20) 


can  be  Introduced,  and  Eq.(3.19)  becomes 


„ b 2 2 .2  \ 

U .^(Vo  cos  a-Xj^  ) 


(3.21) 


It  is  now  necessary  to  evaluate  the  critical  conditions  as  indicated  in 
the  discussion  following  Eq.(3.16).  The  first  condition  is  simply 


I 


(3.22) 


the  second  condition  is  the  value  of  x^  when  y^"0  . The  expression 
for  y^  is  determined  from  Eq.(3.12)  which  can  be  written 


y’t  + - 0 


(3.23) 


where  g ■ O-Kil^b/m.  The  solution  to  this  differential  equation  Is 

* r 


. I 

: i 


* Cj^sin  t + C 2^08  /?t 


(3.24) 


subject  to  the  initial  conditions 


I - 0 and  y^ 
t-0  t-0 


V slna 
o 


(3.25) 


Hence 


V sina 


. -0_ 


sin  t 


(3.26) 


and  it  is  seen  that  y -0  when  t - where  n»0,+l,+2...  Flnnle 

selects  the  value,  explicit  expression  for  x can  be 

/g  t 

obtained  from  Eq.(3.l8),  so 


- X, 


t«t  t-0 

c 


— V slna 
g o 


|cos«^) 


- cos/^t 


t-t 


t-0 


- V cosa 
o 


2r 

-r-  V sina 
g o 


(3.27) 


t-t 
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The  two  conditions  represented  by  Eq.(3.22)  and  (3.27)  can  now  be 
substituted  Into  Eq.(3.21)  to  give  the  final  result. 


h „ 2 2 

U - V cos  o 
o 


bV 


u> 


(sln2a  - sln^a) 


Introducing  the  expressions  for  F and  g , 


mV  - 

o 2 

u*  ;r  cos  a 


(3.28) 


(3.29) 


(3.30) 


(3.31a) 


(3.31b) 


using  the  notation  from  Finnle  (1972)  where  P**K/  ^1+  — j . In  his 
work  Finnle  has  provided  some  rationale  from  an  analogy  with  metal  cut- 
ting for  letting  i/^2  and  K-2  . The  value  ip^2  is  implicitly  incor- 
porated into  the  results  presented  in  Finnle  (1972).  In  1958  Finnle 
2 2 

selected  I “ , while  I - is  used  in  Finnle  (1972);  otherwise, 

the  expressions  are  identical.  These  mass  moments  of  Inertia  do  not 
relate  to  a particular  particle  geometry  but  appear  to  be  determined 
on  the  basis  of  Finnle 's  Intuition  as  to  the  forms  which  would  be 
representative  of  the  angular  particles  in  his  model. 
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The  expressions  In  Eq.(3.30)  and  (3.31),  or  equivalently  Eq.(3.28) 
and(3.29),  provide  the  functional  relations  for  the  angular  dependence 
of  the  steady-state  erosion  rates.  Indicate  the  volume  removed  is  pro- 
portional to  the  kinetic  energy  of  the  Impacting  particle,  and  contain 
only  one  material  parameter,  the  flow  stress  0^  , representing  the  re- 
sponse of  the  target.  The  general  features  of  the  curves  corresponding 
to  Eq.(3.30)  and  (3.31)  will  be  examined. 

The  curves  corresponding  to  Eq.(3.30)  and  (3.31)  are  plotted  In 
Figure  3.12  for  0<a£90  . The  solid  curve  is  the  one  used  by  Flnnle 

to  correlate  erosion  data;  the  conditions  for  the  point  at  which  the 
two  curves  cross  and  the  maximum  in  the  erosion  curve  can  be  obtained 
from  the  explicit  relations  for  the  angular  dependence.  The  maximum 
value  Is  found  by  differentiating  Eq.(3.31)  with  respect  to  a and 
setting  the  result  to  zero,  then  the  value  of  a corresponding  to 
maximum  erosion  Is  determined  from 


tan2o  ■ P 
m 


(3.32) 


The  expressions  In  Eq.(3.30)  and  (3.31b)  are  equal  for  the  value  of 


a determined  from 
o 


tana  ■ P/2 
o 


(3.33) 


The  relative  locations  of  the  angles  represented  by  Eq.(3.32)  and 

(3.33)  are  shown  In  Figure  3.12.  The  variation  in  the  magnitude  of 

these  angles  with  Flnnle 's  specifications  for  I can  be  evaluated. 

For  K«2  the  maximum  erosion  occurs  when  a -16*- 51'  and  the  erosion 

m 

curves  are  equal  when  a -18*-26'  with  I - . The  corresponding 

® 2 2 

values  of  a when  I ■ -y  are  13*-17'  and  14*-2'  . 
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FINNIFS  EROSION  CURVE 


a (d«gr«M) 


Figure  3.12 


Fom  of  the  Erosion  Curve  for  Oblique  Impacts 
on  Ductile  Targets  based  on  Finnic 's  Analysis 

inr^ 

when  K-2  and  I - 1/3  ^ . 


The  general  range  of  applicability  of  Flnnle's  model  is  taken  to 
be  for  rigid  particles  greater  than  100  pm  and  for  values  of  a from 
0 to  45*  (Finnie  and  McFadden,  1978;  Flnnle,  et  al.,  1979).  The 
governing  equations  determine  two  conditions:  the  completion  of  a cut 
in  the  target  material  at  low  angles  of  attack,  Eq.(3.31),  and  embed- 
ding of  the  particle  at  larger  angles  of  attack,  Eq.(3.30).  The  tra- 
jectory followed  by  the  particle  in  passing  through  the  material  deter- 
mines the  quantity  of  material  removed.  The  width  b of  the  contacting 
face  of  the  particle  does  not  appear  in  the  final  form  of  the  equations. 

Finnie  and  McFadden  (1978)  modified  Finnie' s original  derivation 
in  an  effort  to  obtain  an  analytical  expression  for  the  rate  of  mate- 
rial removal  which  agrees  more  closely  with  the  experimentally-deter- 
mined velocity  exponent  which  for  ductile  metals  falls  into  the  range 
from  2 to  3.  In  view  of  Hutchings  (1977,  1979)  more  accurate  represen- 
tations (Figure  3.8)  of  the  conditions  which  are  likely  to  arise  for 
angular  particle  Impacts,  the  modification  Introduced  by  Flnnle  and 
McFadden  is  essentially  an  artifice  to  achieve  their  objective  formu- 
lated so  that  the  governing  equations  can  be  solved  analytically.  The 
only  change  introduced  into  the  initial  derivation  is  that  the  force 
acting  on  the  contact  area  is  applied  to  the  centroid  of  the  contact 
area  instead  of  at  the  cutting  tip.  The  resulting  analytical  expres- 
sion is  then  equated  to  the  power  law  relation  for  volume  removal  which 
is  essence  is  equivalent  to  assuming  the  result  they  want  to  derive. 

The  direction  of  the  applied  force  and  orientation  of  the  contact 
face  are  held  constant  throughout  the  time  the  particle  is  in  contact 
with  the  target.  This  change  does  not  influence  Eq.(3.11)  or  (3.12) 
but  Eq.(3.13)  is  modified  as  follows 

• • 

!♦  + F^(r-y)  - 2KF^y  - 0 (3.34) 
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where  y Is  the  distance  to  the  centroid  of  the  contact  area.  The 
same  procedure  Is  followed  for  solving  the  revised  equations  as  has 
been  outlined  for  the  Initial  formulation.  The  new  equation  correS' 
ponding  to  Eq.(3.31b)  Is 


sln2a-4sln  a I -t- 


where  as  before  3 ~ 
(compare  Eq. (3.31b)) 


and  a value  of  P 


has  been  assumed 


This  expression  for  the  volume  removed,  U 
the  general  form 


u - (A+BV  )V 


On  the  basis  of  the  experimental  data,  Flnnle  and  McFadden  assume 


The  velocity  exponent  Is  determined  by  equating  the  analytical  result 
to  the  power  law  relation  representative  of  the  experimental  mass-loss 
data  Indicated  In  Eq.(3.37),  so 


The  explicit  form  of  Eq.(3.38)  expressed  as  the  ratio  of  the  volume 
removed  at  two  different  Impact  velocities,  and  , Is  then  used 
to  evaluate  the  exponent  n , 


sin2a-4sln  ort-(15/4)Ti(U,/U,  )X,sln  a 


sln2a-4sln  oH-(15/4)TrX,8ln  a 


UjSina 

where  X,  ■ which  is  the  maximum  depth  of  cut  relative  to  the 

^ ^r 
particle  radius. 

According  to  Flnnle  and  HcFadden  the  value  predicted  for  the 
effective  exponent  n for  a given  a and  X^  is  quite  insensitive 
to  the  ratio  of  the  impact  velocities,  so  the  volume  removed  plotted 
against  Impact  velocity  should  be  essentially  a straight  line  on  a 
log- log  scale. 

This  is  the  result  when  the  particle  leaves  the  surface  and 
applies  for  angles  of  attack  less  than  approximately  20® . The  eval- 
uation of  the  expression  for  the  volume  removed  corresponding  to  the 
condition  when  the  particle's  trajectory  terminates  in  the  target  is 
somewhat  more  involved.  Representative  values  of  n were  not 
evaluated  by  Flnnle  and  McFadden  for  this  latter  condition. 

Flnnle's  analysis  (Flnnle,  1958)  for  ductile  materials  was  the 
first  of  its  kind  in  the  solid  particle  erosion  literature  and  there- 
fore deserves  recognition  for  the  contribution  made  in  providing  a 
basis  for  correlating  erosion  data.  However  the  restriction  of  being 
able  to  integrate  the  equations  of  motion  analytically  removes  many 
features  of  the  dynamic  penetration  process  occurring  for  oblique  im- 
pacts of  angular  particles  and  may  no  longer  be  a realistic  objective. 
A velocity  exponent  greater  than  2 comes  out  of  Hutchings  numerical 
analyses  which  describe  more  representative  Impact  conditions  fairly 
directly  as  opposed  to  the  imposition  of  the  highly  idealized  impact 
conditions  assumed  by  Flnnle  and  McFadden  (1978)  in  order  to  obtain 
an  analytic  solution. 
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3.4  Bitter 


Flnnle's  theory  of  erosion  of  ductile  materials  Is  based  on  a 
model  of  the  cutting  action  by  the  Impacting  particle  and  therefore 
did  not  account  for  material  removal  for  particle  Impacts  at  normal 
Incidence.  The  experimental  data  shows  that  material  removal  does 
occur  for  angular  particle  Impacts  on  ductile  materials  at  normal 
Incidence  which  Is  a limitation  on  the  application  of  Flnnle's  cor- 
relation to  values  of  a greater  than  45  * as  previously  noted. 

Bitter  (1963a)  provided  a fairly  detailed  analysis  for  the  normal 
collision  of  an  elastic  sphere  on  an  elastic-plastic  half-space  as  a 
model  for  material  removal  In  a solid  particle  erosive  environment. 

The  main  elements  of  Bitter's  approach  will  be  dfP'.rlbed,  but  refer- 
ence to  the  original  publication  is  advised  in  order  to  follow  the 
details  of  Bitter's  discussions  of  the  material  response.  The  general 
features  of  Bitter's  concept  of  the  material  removal  mechanisms  are 
presented  as  Interpreted  by  the  author,  but  the  physical  process  by 
which  material  Is  removed  Is  not  reflected  In  the  final  seml-emplrlcal 
correlation. 

Bitter  divides  the  general  erosion  process  Into  two  categories: 
deformation  wear  and  cutting  wear.  Flnnle's  analysis  of  cutting  wear 
was  already  published,  however  Bitter  developed  an  analytical  approach 
for  deformation  wear  as  well  as  an  alternative  derivation  of  the  cut- 
ting wear  analysis.  Bitter's  analyses  are  for  a single  spherical  par- 
ticle Impacting  an  elastic-plastic  target  material.  The  deformation 
wear  analysis  is  based  on  an  energy  balance  whereby  the  kinetic  energy 
of  a single  Impacting  particle  Is  divided  Into  recoverable  elastic  and 
elastic-plastic  deformatlonal  energy  and  non-recoverable  plastic  de- 
formation. The  cutting  wear  relations  are  obtained  from  the  equations 
of  motion  for  a single  spherical  particle  striking  an  elastic-plastic 
target  obliquely.  Bitter  (1963a, b)  notes  that  both  de  Haller  (1939) 
and  Wei  linger  (1949)  recognized  that  the  material  removal  process  due 


to  multiple  particle  Impacts  was  Influenced  by  the  hardness  of  the  target 
material  (the  flow  stress  In  Flnnle's  analysis  for  ductile  targets)  and 
the  Impingement  angle. 

3.4,1  Deformation  Wear  Model 

Bitter  (1963a)  described  the  material  removal  mechanism  due  to 
deformation  wear  as  a process  of  strain  hardening  and  embrittlement  of 
the  ductile  target.  He  argues  that  no  material  removal  will  occur  when 
the  states  of  stress  In  the  target  material  do  not  exceed  the  elastic 
limit  except  for  fatigue  damage  due  to  the  repetitive  Impact  cycles. 

When  the  particle  Impacts  are  severe  enough  to  produce  plastic  deforma- 
tion, the  repetitive  Impacts  produce  a plastically  deformed  surface  layer 
which  becomes  work  hardened  and  brittle  and  loses  Its  capacity  for  fur- 
ther plastic  deformation.  During  the  embrittlement  process  the  elastic 
limit  Increases  until  It  represents  the  strength  of  the  target  material. 
Stresses  In  excess  of  the  elastic  limit  for  the  strain-hardened  material 
can  no  longer  be  tolerated  through  further  plastic  deformation  and  there- 
fore fracture  Is  Initiated  apparently  In  a brittle  mode.  Bitter's  analy- 
sis was  strongly  Influenced  by  Davies'  analysis  of  a solid  sphere  Im- 
pacting a ductile  half-space  using  the  Hertzian  theory  of  Impact  (Davies, 
1949).  Using  either  the  Tresca  or  von  Mlses  yield  condition,  Davies  found 
that  the  maximum  principal  stress  difference  (onset  of  yielding)  occurred 
along  the  axis  of  symmetry  at  a point  approximately  one-half  the  contact 
radius  for  a steel  ball  Impacting  a steel  plate.  Thus  locating  the  maximum 
stress  condition  within  the  target  for  the  onset  of  yielding.  Bitter 
equates  this  to  the  point  at  which  the  strength  of  the  material  Is  ex- 
ceeded and  envisions  cracks  developing  parallel  to  the  contact  area  In 
the  plastically-deformed  surface  layer.  Material  fragments  are  assumed 
to  be  removed  In  the  form  of  flakes.  The  mathematical  development  Is 
guided  by  the  conceptual  model  describe'!'  above. 


r 


) 

I 

! 


The  general  single  particle  Impact  at  normal  Incidence  Is  shown 
In  Figure  3.13  as  Idealized  by  Bitter.  For  low  Impact  velocities  the 
Impact  Is  purely  elastic  and  the  Hertzian  theory  of  Impact  Is  assumed 
to  be  applicable  (Appendix  A) . The  applied  pressure  distribution  Is 
paraboloidal.  VRien  the  elastic  load  limit  Is  exceeded,  plastic  stresses 
develop  as  shown  In  Figure  3.13b  for  an  elastlc-perfectly  plastic 
material. 

Bitter's  expression  for  material  removal  due  to  deformation  wear 
Is  based  on  a partitioning  of  the  Initial  kinetic  energy  of  the  Impacting 
particle  Into  plastic,  elastic-plastic,  and  plastic  deformations  of  the 
target  material.  The  energy  balance  Is 

Q . i .Vj  . (3.40) 

The  term  Is  the  fully  recoverable  elastic  deformation  energy  asso- 

ciated with  the  loading  function  up  to  the  load  level  at  which  yielding 
In  the  target  Initiates.  Up  to  this  load  level  the  Impact  Is  described 
In  terms  of  the  Hertzian  theory  of  Impact  which  has  a major  role  In 
Bitter's  analysis.  The  elastic-plastic  response  of  the  target  Is  con- 
sidered to  comprise  a further  elastically  recoverable  component,  Q , 

pe 

and  a pure  plastic  (permanent)  deformation  of  the  target  absorbing  the 
energy  . The  rebound  velocity  of  the  spherical  projectile  from  the 
elastic-plastic  target,  V , Is  governed  by  the  elastic  deformatlonal 
energy  stored  In  the  target,  so 

i “V^  - Qg  + Qpg  (3.41) 

Using  the  explicit  representations  for  Q^,  and  Qp  , Bitter 

Is  able  to  show  that 

On,  ■ 1.936VQ.Q„  ■ 


(3.42) 


Eliminating  Qp^  from  Eq.(3.40)  , 


Q - Qe  + + Qp  - o.oeAVQ^Qp 


(3.43) 


which  Bitter  approximates  by 


(3.44) 


Bitter  defines  the  measure  of  deformation  wear,  , as 


(3.45) 


where  C Is  an  empirical  parameter  representing  the  amount  of  energy 
required  to  remove  one  unit  volume  of  material,  is  due  solely  to 
the  energy  absorbed  by  plastic  deformation  of  the  target  material. 


With  the  Introduction  of  the  empirical  parameter  5 > Eq.(3.45)  Is 
essentially  no  more  than  an  assumed  representation  for  the  deformation 
wear  response.  Substituting  Eq.(3.44)  into  (3.45),  using  the  definition 
of  Q In  Eq.(3.40)  and 


Q - y mV  ^ 
e 2 e 


(3.46) 


which  defines  the  portion  of  the  Initial  kinetic  energy  contributing  to 
the  purely  elastic  deformation  of  the  target  material,  Eq.(3.45)  becomes 


%- 


M(V  sina-V  ) 
2 o e 


(3.47) 


It  Is  assumed  that  only  the  normal  component  of  the  Initial  particle  Im- 
pact velocity,  , will  be  effective  in  producing  deformation  wear  and 
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that  the  steady-state  erosion  rate  Is  simply  a linear  accumulation  of 
the  damage  produced  by  a single  particle  with  the  Introduction  of  the 
total  mass,  M , for  all  the  particles  striking  the  target.  No  erosion 
takes  place  when  V^slna£V^  . The  original  definition  of  deformation 
wear,  Eq.(3.45),  has  now  been  generalized  to  oblique  Impacts  and  to 
multiple  spherical  body  Impacts. 


Although  the  Hertzian  theory  Is  used  extensively  In  Bitter's 
analysis,  only  the  velocity  Is  related  explicitly  to  the  material 

properties  of  the  Impinging  spherical  bodies  and  the  target  In  accord- 
ance with  the  following  relation 


V 

e 


1 


Pi 


2 

(kj^+k2) 


(3.48) 


This  expression  Is  derived  from  the  straightforward  application  of  the 
Hertzian  theory  of  Impact  which  Is  taken  to  be  valid  up  to  the  value  of 
the  average  applied  pressure,  , at  which  yielding  within  the  target 

material  Is  Initiated.  (Bitter  refers  to  as  the  elastic  load  limit 

and  denotes  It  by  the  symbol  y .)  The  remaining  parameters  are  the  den- 
sity of  the  Impinging  sphere,  , and  the  elastic  moduli  of  the  sphere 
and  the  target,  k^  and  k2  , as  defined  In  Eq.(A3)  of  Appendix  A. 

Despite  the  lack  of  a sound  basic  derivation  of  the  deformation 
wear  expression  In  Eq.(3.47),  It  has  been  one  of  the  major  correlating 
functions  for  solid  particle  erosion  test  results. 

Rebound  tests  for  hardened  steel  balls  Impacting  plates  of  various 
materials  are  used  to  Indicate  the  Independence  of  on  the  particle 

size  (mass).  Hardened  steel  balls  of  different  mass  were  dropped  from 
a constant  height  of  85.7  cm  and  the  rebound  height  was  measured.  The 
rebound  height  Is  taken  to  be  dependent  on  the  recoverable  elastic  energy 
In  Eq.(3.40)  which  Is  the  sum  of  and  . The  energy  balance  can 
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now  be  written 


1 

2 


(3.49) 

e pe  p 

2 12  1 ^ 

« ^mV''  + i m(V  -V^) 
o z z o e 

(3.50) 

where  the  definitions  In  Eq.(3.40)  and  (3.41)  and  Bitter's  result  In 
Eq.(3.44)  have  been  used.  Bitter  obtains  from  Eq.(3.50) 


2 2 
V - 2V  V - V 
o e e 

which  can  be  solved  for  V , 

e * 


(3.51) 


(3.52) 


where  e » •—  is  the  coefficient  of  restitution.  According  to  this 
*o 

expression  the  functional  dependencies  of  are  Identical  to  those 

of  e (Goldstein,  1960?  Adler,  1961).  When  e=l  , the  collision  is 
perfectly  elastic  and  , so  according  to  Bitter's  hypothesis  no 

material  loss  will  occur.  All  of  the  Initial  kinetic  energy  will  be 
dissipated  In  plastic  deformation  when  e=0  . In  general  e varies 
with  both  the  mass  of  the. Impacting  particle  and  its  velocity,  as  well 
as  numerous  other  properties  of  the  particle  and  target  material:  e 
is  a dynamics  parameter  which  represents  all  of  the  modes  of  energy 
dissipation  which  occur  for  a particular  impact  condition  in  real 
materials . 


The  experimental  determinations  of  V are  used  to  evaluate  o . 

e eS, 

from  Eq.(3.48).  The  empirical  parameter  ^ Is  determined  from  the  ero- 
sion test  data.  Bitter  elected  to  use  the  experimental  erosion  rate 
(mass  removal /mass  of  Impinging  particles)  at  a-bO”  to  obtain  the 
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value  of  C . Eq.(3.47)  then  provides  the  explicit  variation  in  angular 

dependence  for  the  erosion  rate  associated  with  solid  particulates  im- 


pacting at  a specified  velocity  on  a specific  material. 


W. 


"2% 


M 


^2  2 
-rp  (V  sina-V  ) 
It,  o e 


(3.53) 


where  P2  is  the  density  of  the  target  material.  W-  is  expressed  in 


terms  of  the  mass  of  material  removed  per  unit  mass  of  impacting  parti- 
cles. Some  material  properties  are  included  explicitly  in  , Eq.(3.48), 
however  the  material  properties  influencing  the  target's  response  in  an  ero- 
sive environment  are  primarily  represented  in  the  parameter  ^ . Bitter 
provides  a brief  discussion  of  this  factor. 


3.4.2  Cutting  Wear  Model 


Bitter  derived  the  equations  governing  cutting  wear  on  the  basis 
of  an  elastic-plastic  analysis.  Bitter  considered  the  two  particle  im- 
pact conditions  which  Flnnle  described;  namely,  the  particle  cuts  through 
the  surface  with  a nonzero  final  velocity,  or  it  becomes  embedded  in  the 
surface  before  completion  of  the  cutting  action.  In  the  first  case 


u . ^2.,,  2 2 .2. 

"Cl  “"*f> 


(3.54) 


is  the  tangential  velocity  which  the  particle  has  as  it  leaves 


where 

the  surface  and  (J)  is  the  cutting  wear  factor  which  is  the  energy  re- 
quired to  "scratch"  a unit  volume  of  material  from  the  surface.  Again, 


'^Cl  expressed  in  terms  of  the  mass  of  material  removed  per  unit  mass 


of  impacting  particles.  Bitter  does  not  provide  any  additional  informa- 
tion on  the  cutting  wear  factor  except  that  it  depends  on  the  mechanical 


properties  of  the  target.  For  the  second  Impact  condition  when  x^-0 


before  completion  of  the  cut. 


w * 

C2 


i 

COS  a 


2(J) 


(3.55) 
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Bitter  finds  that  a correction  factor  is  required  in  Eq.(3.55)  if  the 
condition  x»0,yK)  is  achieved  when  the  maximum  depth  of  cut  occurs. 


Bitter’s  cutting  wear  analysis  involves  the  determination  of  the 
final  horizontal  component  of  velocity,  , in  Eq.(3.54)  and  the  indi- 
cated correction  factor  for  Eq.(3.55).  The  approach  used  is  based  on 
the  oblique  impact  of  a spherical  particle  which  is  equivalent  to  an 
angular-shaped  abrasive  particle.  The  size  of  the  irregular-shaped  body 
is  approximated  by  a spherical  radius  R . The  particle  is  assumed  to 
have  rounded  edges  and  to  be  fairly  equlaxed.  A smaller  rounded  portion 
of  the  particle  with  a radius  of  curvature  r penetrates  the  ductile 
target  and  the  mass  of  the  particle  is  expressed  in  terms  of  this 
apparent  radius,  so 


4 , 3 

m = ir  p'  r 

p'  = (R/r)^  p 


(3.56) 


where  p’  is  the  fictitious  density  of  the  equivalent  particle  and  p 
is  the  density  of  the  actual  particle.  The  velocity  x^  is  then  deter- 
mined from  the  equations  of  motion  for  the  oblique  impact  of  the  equiva- 
lent spherical  body. 


mx  = -A(J)  = 


F 

X 


(3.57) 


my 


2^22 

TTr  a 0+  - a .= 

p eJl  e 3 eil 


F 

y 


(3.58) 


where  F^  is  the  horizontal  shearing  force  resisting  motion  tangential 
to  the  surface  of  the  target  by  the  spherical  body  and  F^  represents 
the  elastic-plastic  penetration  force  after  the  maximum  Hertzian  pressure 
exceeds  the  elastic  stress  limit,  , of  the  target  material.  The  area 
A is  for  the  vertical  cross  section  of  that  part  of  the  particle  which 
has  penetrated  the  target  (roughly  equivalent  to  the  region  shown  in 


Figure  Bl,  Appendix  B)  , and  4)  Is  the  cutting  wear  factor.  The 
vertical  resisting  forces  correspond  to  the  condition  In  Figure  3.13 
already  described  for  the  deformation  wear  analysis.  Bitter's  equations 
of  motion  exclude  frictional  effects  and  rotation  of  the  spherical  body 
as  It  cuts  through  the  surface  of  the  material.  The  depth  of  the  cut  Is 
determined  from  the  normal  velocity  component. 

It  Is  noted  In  passing  that  Surette  (1971)  Identified  the  primary 
mode  of  volume  removal  to  be  failure  of  the  Impacted  surface  In  shear 
due  to  sliding  friction  between  the  Impinging  particle  and  the  surface. 

He  Indicates  that  his  model  Is  based  on  first  principles  and  that  the 
parameters  required  are  combinations  of  readily  obtainable  engineering 
data.  A somewhat  confusing  derivation  of  the  governing  relations  Is 
presented. 

Surette  writes  the  equations  of  motion  for  a rigid  spherical  body 
Impacting  a ductile  target  obliquely,  however  these  equations  require 
some  Interpretation  due  to  the  lack  of  precise  meanings  assigned  to  the 
notation  used.  Referring  to  Figure  3.14, 

mx  + ® (3.59) 

my  + N = 0 (3.60) 

1*6  - f (a  A.)R  - 0 (3.61) 

s s a 


Bitter  states  A>4/3r^H,  where  H Is  the  plastic  penetration  depth  In 

the  target,  however  the  author  was  not  able  to  derive  this  expression 
for  the  area  A.  It  Is  not  obvious  what  approximations  were  used  In 
Bitter's  derivation. 
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where 


N Is  the  normal  force  developed  at  the  Interface  and  Is 
stated  explicitly  in  terms  of  the  Hertzian  relation 
for  static  indentations  in  elastic  materials 

f^  is  the  coefficient  of  sliding  friction 

A^~2irRy  is  the  area  of  the  contact  surface  for  a rigid 

sphere  impressing  an  elastic  half-space  (Eq.(B3) 
in  Appendix  B 

2 2 

I ■ mR  is  the  moment  of  inertia  for  a sphere 
m is  the  mass  of  the  sphere 

R is  the  effective  radius  of  the  particle  allowing 
for  irregular-shaped  particles 

The  parameter  is  simply  defined  as  the  shear  stress,  so 

represent  a force.  It  is  not  clear  why  this  poorly  defined  shearing 
force  is  then  multiplied  by  the  coefficient  of  sliding  friction.  This 
definition  of  the  force  generated  during  the  collision  is  inconsistent 
with  the  available  analyses  for  elastic  bodies  in  contact  and  subjected 
to  tangential  loadings  (Mlndlln,  1949;  Johnson,  1955).  The  problem  of 
a sphere  impacting  an  elastic  surface  obliquely  has  been  more  accurately 
treated  by  Mamoun  (1975b)  in  relation  to  development  of  an  erosion  model. 

We  note  that  the  equations  of  motion  suggested  by  Burette  (1971) 
can  be  formulated  as  a special  case  of  Hutchings'  equations  when 
3+^“0  in  Eq.(3.2),  then 

mx  - -pN 

my  - -N  (3.62) 

•• 

10  - RpN 

The  imposed  condition  corresponds  to  the  force  acting  along  the  y axis 
which  is  representative  of  a normal  collision  rather  than  representing 
the  contact  area  which  would  result  as  a rigid  bead  ploughs  through  a 
ductile  target  during  an  oblique  impact.  In  Burette's  analysis  the 
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equations  of  motion  are  very  Idealized  for  oblique  Impacts  and  do  not 
conform  to  a realistic  physical  model  of  the  collision  process.  However 
this  simplification  permits  the  equations  of  motion  to  be  solved  anal- 
ytically. The  lack  of  clarity  In  Burette's  work  eventually  frustrated 
this  author  and  his  work  Is  left  to  others  to  decipher. 

Comparing  Eq.(3.58)  and  (3.60)  It  Is  seen  that  both  Bitter  (1963b) 
and  Burette  define  the  normal  force  In  terms  of  the  Hertzian  Impact  pres- 
sure but  according  to  Eq.(3.57)  and  (3.59)  define  Independent  forces  for 
the  resistance  which  develops  parallel  to  the  surface  of  the  target. 

This  Is  contrary  to  Burette's  model  but  Is  consistent  with  Bitter's  ap- 
proach where  the  particle  actually  penetrates  the  surface.  Bitter  as- 
sumes that  no  material  Is  removed  If  the  oblique  impact  Is  purely  elastic. 

Heuristic  arguments  are  used  with  the  general  equations  for  the 
assumed  model  to  obtain  the  final  expressions  for  target  volume  removal 

0 / 

2MC(V  slna-v  ) V coso  2MC  (V  slna-V  ) 

^C1  - — • ° — (p  a<a  (3.63) 

/V  slna'  V slna  ° 
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u_„  “ V ^cos^a-K, (V  slna-v  ) 

C2  2(j)  L o 1 o e . 


a>a  (3.64) 
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where 


0.288 


ei  ''eZ 


(3.65) 


2/0  / 1-v/  1-v  ^ V 

'l  - “et  (fr) 


(3.66) 


and  Vg  Is  defined  In  Eq.(3.48).  When  and  K,  the  angle  a^ 

can  be  found  from 
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DetermlnaClon  of  the  range  of  application  of  Eq.(3.63)  and  (3.64) 

is  not  as  simple  as  for  Flnnle's  analysis.  Bitter's  expression  In 

Eq.(3.63)  shows  the  dependence  of  the  erosion  rate  on  Impact  velocity 

5/2 

to  be  essentially  since  the  second  term  on  the  right-hand  side 

will  generally  be  small  compared  to  the  first  term.  Thus  a velocity 
exponent  different  than  2 Is  a result  of  Bitter’s  analysis.  Removal  of 
this  velocity  squared  restriction  was  the  objective  for  the  models  and 
corresponding  analyses  undertaken  by  Sheldon  and  Kanhere  (1972)  and 
Flnnle  and  McFadden  (1978).  It  Is  also  Interesting  to  note  that  the 
first  and  dominant  term  in  Eq.(3.63)  can  be  evaluated  explicitly  without 
knowing  the  empirical  factor  ((i  . 


The  experimental  erosion  rate  data  obtained  as  a function  of 
Impingement  angle  Is  now  represented  by  the  sum  of  the  contributions 
due  to  deformation  wear  and  cutting  wear.  Thus,  from  Eq.(3.53),  (3.63) 
and  (3.64), 


when  a < a 


when  a > a 


In  order  to  fit  the  erosion  data  the  empirical  parameters  ^ and  4* 
have  to  be  determined.  Two  experimental  data  points  are  required.  From 
Eq.(3.S3)  It  Is  seen  that  ^ can  be  evaluated  directly  when  a*90”  . 

The  curve  can  then  be  constructed.  A second  data  point  for 

a <a«90*  provides  the  equation 


can  then  be  constructed  if  0 


Is 


for  the  evaluation  of  (^  . 
known  for  the  target  material.  The  value  of 
from  the  Intersection  and  U_.  . 


el 

Is  found  graphically 


Bitter  (1963b)  shows  that  Eq.(3.64)  Is  equivalent  to  Flnnle's 
Eq.(3.30)  when  . However  he  notes  that  Flnnle's  Eq.(3.31)  for 

a<a^  Includes  the  Implicit  assumption  that  the  time  required  for  pene- 
tration of  the  particle  Is  equal  to  the  time  required  for  Its  egress 
which  can  be  seen  upon  differentiating  Eq.(3.26).  Bitter's  observation 
Is  that  If  the  normal  velocity  of  the  parti -''•e  becomes  zero  before  the 
tangential  velocity  component,  the  particle  Is  pressed  out  of  the  target 
by  the  elastic  reaction  forces  while  Its  cutting  action  continues  (due 
to  the  non-zero  tangential  velocity).  This  Is  consistent  with  Flnnle's 
analysis  but  Is  not  consistent  with  the  assumption  that  the  target  mate- 
rial Is  perfectly  plastic.  For  this  case  there  Is  no  elastic  restoring 
force;  It  can  be  presumed  the  cutting  action  of  the  particle  would  con- 
tinue tangential  to  the  target's  surface  until  the  tangential  velocity 
(and  corresponding  kinetic  energy)  was  dissipated.  Bitter's  formulation 
of  the  problem,  Eq.(3.57)  and  (3.58)  accounts  for  this  condition.  Bitter 
thus  comments  that  Flnnle  has  utilized  a relation  describing  the  contact 
duration  for  the  collision  of  a spherical  body  with  a plane  surface, 
although  his  model  Is  based  on  oblique  Impacts  by  angular  particles. 


3.5  Nellson  and  Gilchrist 


Nellson  and  Gilchrist  (1968)  devised  a purely  empirical  scheme  for 
fitting  the  experimental  erosion  data  (erosion  rate  as  a function  of  Im- 
pact angle)  based  on  their  own  test  results  and  the  analyses  of  Flnnle 
and  Bitter.  Their  correlation  procedure  will  be  described. 

Nellson  and  Gilchrist  note  that  the  following  observations  should 
be  Included  In  any  correlation  of  the  experimental  data. 


The  kinetic  energy  of  the  impacting  particles  normal  to 
the  target  is  imparted  to  the  specimen  and  accounts 
for  deformation  wear. 


There  is  a threshold  velocity  component  normal  to  the 
surface,  k , above  which  erosion  occurs:  principally 
for  deformation  wear.  Hie  kinetic  energy  of  the  impact- 
ing particles  parallel  to  the  target  is  responsible  for 
cutting  wear. 

There  is  a critical  Impingement  angle,  , above  which 
the  particles  come  to  rest  in  the  target  and  below  which 
the  particles  leave  the  surface  with  a velocity 
parallel  to  the  surface. 

Nellson  and  Gilchrist  further  observe  that  tests  on  ductile  materials 

using  a constant  Impact  velocity  Indicate  that  as  the  impingement  angle, 

a , increases  from  zero,  the  erosion  rate  increases  initially  at  a rapid 

rate  which  decreases  as  the  impingement  angle  Increases.  They  assume  ' 

that  the  cutting  action  of  the  particles  when  0(<a  can  be  represented  ’ 

by  a function  of  the  form 

2 ' ' 
) - 1-sln  n a (3.71) 

where  n is  a constant.  Then, 

Vp  - when  o - 0 (3.72) 

V • 0 and  ot  ■ ~ when  a ■ a (3.73) 

P O M O 

The  above  observations  are  contained  in  the  following  relations. 
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V 

E_ 

Vcosa 


(3.74) 


W. 


■ -^(v  slno-ic)^ 
Z::  o 

or,  introducing  Eq.(3.71), 

W_  - ^(V  sina-K)^ 
T O 

and 


12  2 2 
+ -^(V  cos  a-V  ) 

2*  0 p 


12  2 

+ cos  asln  n a 

2v  o 


12  2 
+ V cos  a 
24)  o 


when  a<a 
— o 


a>a 
— o 


(3.75) 

(3.76) 


where  <^  denotes  the  amount  of  kinetic  energy  which  has  to  be  absorbed 
to  remove  one  unit  mass  of  material  from  the  target  by  the  cutting  wear 
mechanism  and  H is  the  corresponding  parameter  for  deformation  wear. 
Neilson  and  Gilchrist  described  the  general  features  of  Eq.(3.75)  and 
(3.76)  and  provided  a procedure  for  evaluating  the  parameters  n,<(  and 
E when  <>0.  A self-contained  approach  for  carrying  out  this  procedure 
will  be  described  for  use  by  the  Interested  reader. 


First  consider  erosion  rate  curves  representative  of  ductile 
response  when  the  cutting  mechanism  dominates.  In  terms  of  Neilson  and 
Gilchrist's  approach,  . The  maximum  erosion  rate  in  the  ductile 

response  curve  will  be  determined.  Differentiating  Eq.(3.75)  with 
respect  to  a and  setting  the  result  equal  to  zero,  yields 


sin  na 

m 


ncos  na 
2tan  a 

m 


m 


(3.77) 


where  a < a is  the  value  of  a at  which  the  maximum  erosion  rate 
m — o 

occurs.  The  values  of  are  plotted  in  Figure  3.13  as  a function  of 

n and  ^/E  . When 
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♦/E  - 1 , 

a “ a - 

m o 2n 


(3.78) 


The  procedure  for  evaluating  the  parameters  n,<I>  and  E when 

4/H  ^ 1 requires  experimental  data  to  establish  the  erosion  rate  when 

Cl  * tt/Z  , the  value  of  a , and  the  maximum  In  W . The  steps  In  con- 

m c 

strutting  the  erosion  rate  curve  are  Illustrated  In  Figure  3.16. 


1. 


Evaluate  E from  Eq.(3.76)  when  a 


o 


O 


2 

(exp) 


it/2  , 


2. 


Knowing  E , draw  the  curve  for  the  erosion  rate  associated 
with  the  deformation  mechanism  of  removal. 


W, 


TV  ' 
2 o 


, 2 
sin  a 


3.  Subtract  from  W^(exp)  data  points  to  approximate  , 

the  erosion  rate  due  to  the  cutting  mechanism  of  removal,  so 

A.  Determine  a for  . 

m C 

5.  Using  the  value  of  a I determine  the  value  of  n from 

"'“c 

Figure  3.15  when  4/E  ~ 0 which  corresponds  to  the  condi- 
tion when  the  erosion  Is  due  entirely  to  the  cutting  mechan- 
ism. The  value  of  a (4/E  « 1)  can  now  be  found  from 

o 

Figure  3.15  for  this  value  of  n . 
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6.  Knowing  n and  having  the  experimental  value  of  a for 

n 

, the  value  of  4>/H  can  be  obtained  from  Figure  3.15. 
The  value  of  ^ can  now  be  evaluated  since  H Is  known. 


The  curve  for  based  on  the  parameters  n,4>,H,  and  can  be  con- 

structed as  an  analytical  representation  for  the  experimental  data. 


When  ♦/“  > 1 , — > 0 for  all  a , the  deformation  mechanism 
uu 

dominates  and  the  erosion  rate  curve  displays  behavior  characteristic 
of  brittle  materials.  For  this  case  Nellson  and  Gilchrist  replace  the 
condition  for  determining  a corresponding  to  a maximum  erosion  rate 
with  the  arbitrary  selection  of  the  value  of  a (denoted  by  0)  cor- 
responding to  the  erosion  rate  which  Is  one-half  the  erosion  rate  for 
(X"Tf/2  . Substituting  this  condition  Into  Eq.(3.75)  and  (3.76)  yields. 


cos  0sln  nB 

1/2  - sln^B 


1/2  - sin  B 


when  B < 01 
— o 


when  B > ot 
— o 


(3.79) 


(3.80) 


The  values  of  B are  plotted  In  Figure  3.17  as  a function  of  n and 

along  with  Eq.(3.78).  To  construct  the  empirical  curve  for  experi- 
mental data  displaying  brittle  response,  steps  1 through  5 are  used, 
however  the  experimental  value  of  corresponding  to  B Is  used  with 

the  known  value  of  n to  determine  the  value  of  ♦/£  from  Figure  3.17. 

Nellson  and  Gilchrist  observe  that  $ and  E are  dependent  on 
the  particle  Impact  velocity  and  particle  shape  but  do  not  appear  to 
depend  on  the  particle  size.  It  Is  obvious  that  many  additional  factors 
would  Influence  4,H,  and  n . 
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If  the  parameter  K is  equated  to  Bitter's  determination  of  the 
maximum  velocity,  , for  which  the  target  material's  response  In  elas- 
tic, then  an  explicit  determination  of  K can  be  made  using  Eq.(3.48). 
Typically  the  velocity  for  ductile  materials  Is  quite  small,  so 

neglecting  this  effect  Is  not  very  significant  unless  long-term  erosion 
at  very  low  velocities  Is  of  Interest. 

Prior  to  the  publication  of  Neilson  and  Gilchrist,  Wood  (1966) 
simplified  Bitter's  Eq.(3.69)  to  the  form 

W - -i(V  sina)^  + ~(V  cosa)^  for  a > 15“  (3.81) 

it,  o i(p  o — 

which  he  used  to  correlate  his  erosion  data  for  Impacts  on  metals  when 
15“  ^ OL  ^ 90“ . Wood  evaluated  the  value  of  e and  <|)  for  all  of  his 
test  (lata  which  included  six  quartz  dusts  with  different  size  distribu- 
tions ranging  from  an  average  size  of  3 to  50  ym,  and  nine  metals,  some 
with  different  heat  treatments.  Both  Impact  velocity  (122  to  244  ms  ^) 
and  temperature  (from  ambient  to  648“C  (1200“F))  were  varied.  Wood 
found  that  the  deformation  wear  factor  (fi  and  the  cutting  wear  factor 
5 do  not  show  any  correlation  with  the  Knoop  hardness  or  strain  energy. 

A thirty-fold  difference  was  represented  In  the  hardness  values  for  the 
target  materials  selected  and  a fifty-fold  difference  in  strain  energy. 
The  factors  <t>  and  C did  however  display  a dependence  on  particle 
size  and  temperature:  the  erosion  resistance  first  Increased  with  In- 
creasing temperature  and  then  decreased  as  the  temperature  Increased 
further. 
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Tabakoff  and  Co-workers 


Tabakoff  and  his  associates  have  studied  and  modeled  solid  particle 
flow  trajectories  through  fairly  complex  turbomachinery  (Clevenger  and 
Tabakoff,  1974;  Grant  and  Tabakoff,  1975;  Tabakoff,  et  al.,  1976;  Taba- 
koff and  Hamed,  1977;  Wakeman  and  Tabakoff,  1979).  Their  more  recent 
work  examined  the  effects  of  solid  particle  erosion  on  gas  turbine 
and  compressor  components.  In  order  to  model  the  multiple  Interactions 
a particle  may  have  with  the  stationary  and  moving  surfaces  It  encounters 
as  It  passes  through  a gas  turbine  for  example,  Tabakoff  finds  that  It 
Is  necessary  to  Investigate  the  rebound  characteristics  of  the  particle 
collisions  represented  by  the  coefficient  of  restitution  (Ball  and  Taba- 
koff, 1974),  The  coefficient  of  restitution  is  taken  to  be  a measure  of 
the  kinetic  energy  exchanged  between  the  two  colliding  bodies  which  char- 
acterizes the  particle/target  Interaction.  The  Impact  and  rebound  pro- 
cess Is  statistical  In  nature,  so  statistical  distributions  were  obtained 
experimentally  and  empirical  equations  were  fitted  to  the  data.  In  order 
to  complete  the  general  model  of  turbomachinery  erosion,  it  is  necessary 
to  prescribe  the  amount  of  material  which  will  be  removed  by  particle 
collisions  (Grant  and  Tabakoff,  1973). 

Tabakoff  and  his  co-workers  found  that  it  was  the  coefficient  of 
restitution  defined  In  terms  of  the  tangential  components  of  the  impact 
^ I and  rebound  velocity  which  was  related  to  the  magnitude  of  the  material 

; ' removed:  maximum  erosion  occurred  In  ductile  metals  when  the  coeffl- 

I 

I dent  of  restitution  for  the  tangential  velocity  components  was  minimum, 

f I They  conjectured  that  the  coefficient  of  restitution  based  on  the  normal 

' components  of  the  impact  and  rebound  velocities  represented  the  kinetic 

energy  dissipated  by  plastic  flow  of  the  target  material  without 
^ significant  material  removal. 

I 

I Let  be  the  magnitude  of  Impact  velocity  and  a the  angle  of 

inclination  of  the  velocity  vector  to  the  surface  as  previously  defined. 
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and  be  the  magnitude  of  the  rebound  velocity  and  3 Its 

Inclination  angle  with  the  target's  surface,  then  the  tangential  restl' 
tutlon  ratio,  R_  , Is 


The  experimental  data  for  410  stainless  steel,  T1-6A£-4V,  and  2024 
aluminum  (Ball  and  Tabakoff,  1973;  Grant  and  Tabakoff,  1973)  Indicated 
that 


for  quartz  and  alumina  particle  collisions 


The  erosion  data  for  2024  aluminum  (Grant  and  Tabakoff,  1973)  was 
found  to  be  described  by  the  following  equation 


where 


The  angle  is  the  impingement  angle  for  maximum  erosion,  Kj^,K2,  and 

are  empirical  constants,  and  Is  determined  from  Eq.(3.83). 

Eq.(3.84)  was  formulated  on  the  assumption  that  the  erosion  process  Is 
characterized  by  two  erosion  mechanisms:  one  dominant  for  small  Impinge 
ment  angles  and  the  other  applicable  for  large  Impingement  angles.  It 
can  be  seen  from  Eq.(3.84)  that  the  velocity  dependence  will  be  greater 
than  the  velocity  squared  for  low  angles  of  attack  and  that  the  velocity 


exponent  will  be  4 when  a-90*.  The  velocity  exponent  of  4 was 
found  to  correlate  the  erosion  test  data  for  a 2024  aluminum  alloy  Im- 
pacted by  quartz  particles  (20  to  200  ym)  for  particle  Impacts  ranging 
from  300  to  590  ms  ^ (Grant  and  Tabakoff,  1973,  1975).  For  ductile 
metals  this  Is  a high  exponent  which  Illustrates  the  strong  dependence 
on  the  test  procedure  and  detailed  characterization  of  the  collision 
process,  since  Wood  (1966)  was  able  to  show  that  his  data  for  several 
metals  Including  a 2017  aluminum  alloy  also  using  quartz  particles  with 
a maximum  size  of  50  ym  correlated  with  a velocity  exponent  of  2 when 
a-90®  as  predicted  by  Bitter's  correlation  (Bitter,  1963a).  Wood  used 
the  same  values  of  the  empirical  constant  ^ in  Eq.(3.53)  over  a veloc- 
ity range  from  122  to  244  ms  With  this  wide  variation  in  the  form 
of  the  correlation  for  normal  collisions,  the  magnitude  of  the  discrep- 
ancies possible  between  the  test  procedure  used  and  the  analysis  of  the 
erosive  environment  are  evident.  The  consistency  of  the  data  reported 
by  both  investigators  and  that  of  others  essentially  rules  out  the  pos- 
sibility that  these  differences  are  the  result  of  the  material's  actual 
response  to  the  erosive  environment. 

The  empirical  constants  were  found  for  a number  of  cases  (Grant 

and  Tabakoff,  1975;  Tabakoff,  et  al.,  1979).  For  example,  for  quartz 

particles  impacting  2024  aluminum  these  constants  are  K *3.67x10 

-12 

K2*0.585,  and  K^*6.0xl0  when  the  measure  of  erosion  is  in  mg/gm  of 
impacting  particles  (Grant  and  Tabakoff,  1973,  1975).  A slightly  modi- 
fied form  of  Eq.(3.85)  was  used  for  correlating  the  experimental  data 
obtained  by  Tabakoff,  et  al.,  1979.  There  is  no  explicit  dependence  on 
the  particle  size  in  Eq.(3.84),  however  Grant  and  Tabakoff  (1973)  indi- 
cate there  is  a significant  dependence  on  the  specimen  size  used  in  the 
Tabakoff 's  wind  tunnel  facility.  The  specimen  size  influences  the  aero- 
dynamics of  the  flow  field  in  the  test  section.  Increasing  the  specimen 
width  decreased  the  erosion  rate.  The  specimen  size  effect  appears  to 
be  dependent  on  the  impact  parameters  as  reflected  in  the  particle  inter- 
actions at  the  surface  and  the  nature  of  the  flow  field.  Tabakoff  has 
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not  defined  an  absolute  erosion  rate  for  his  experimental  data  but 
reports  the  specimen  size  with  the  erosion  data.  There  does  not  appear 
to  be  a simple  correlation  between  the  test  results  obtained  from 
specimens  of  different  sizes. 

Tabakoff's  experimental  facility,  test  procedures,  data  analysis, 
and  general  correlation  all  require  further  consideration,  since  they 
represent  a significant  departure  from  the  more  typical  sand  blast  test 
procedures  and  data  analysis.  As  stated  In  Section  2 a greater  aware- 
ness Is  developing  in  solid  particle  erosion  testing  concerning  the  in- 
fluence of  the  test  parameters  on  the  measured  material  removal  rates. 

In  this  regard  Tabakoff's  work  deserves  more  attention. 

Menguturk  and  Sverdrup  (1979)  also  were  concerned  with  the  effects 
of  solid  particle  erosion  in  gas  turbines.  They  Introduced  a generaliza- 
tion of  Bitter's  equations  (Bitter,  1963a, b),  Eq.(3.53)  and  (3.54),  pro- 
ceeding along  the  lines  of  Nellson  and  Gilchrist  (1968),  using  the  final 
forms  given  In  Eq.(3.75)  and  (3.76),  by  raising  the  cutting  wear  and  de- 
formation wear  terms  to  arbitrary  exponents,  m^^  and  m2  , respectively. 
Neglecting  the  maximum  normal  velocity  for  elastic  response  and  setting 
mj^“m2“m  , Menguturi&and  Sverdrup  find 


u ■ K,  (V  c^a)”  sinna  + K,(V  sina)"  for  a < a 
1 o \ 2 o — o 


u - K.(V  cos  (V  + K,(V  sina)™ 

AO  Z O 


(3.87) 


for  a > a 
— o 


where  K^,K2,  and  m are  em^rlcal  constants  which  have  to  be  determined 
As  in  Eq.(3.73)  n ■ 2^  • Thd  procedure  for  determining  the  required 
constants  Is  analogous  ?o  that  outlined  f«^r  evaluating  Eq.(3.75)  and 


3. 7 Mamoun 


Mamoun  (1975a, b,c,d;  1976)  developed  estimates  of  material  removal 
due  to  spherical  particle  Impacts  for  both  ductile  and  brittle  materials 
using  the  Hertzian  theory  of  Impact  and  seml-emplrlcal  relations  for  the 
fatigue  life  as  a function  of  strain  amplitude.  The  ultimate  objective 
of  Mamoun* s modeling  study  was  to  develop  analytical  predictions  for  the 
material  removal  and  mechanical  property  degradation  In  components  sub- 
jected to  the  eroslon/corroslon  environment  associated  with  coal  gasifi- 
cation. The  material  response  was  categorized  as  follows:  ductile  mate- 
rials at  load  levels  below  the  elastic  limit  (case  1),  fully  plastic  de- 
formations In  nons train-hardening  materials  (case  2)  and  strain-hardening 
materials  (case  3),  elastic-plastic  deformations  (case  A),  brittle  mate- 
rials at  load  levels  below  the  fracture  threshold  (case  5),  and  brittle 
materials  above  the  fracture  threshold  (case  6).  Only  the  erosion  due 
to  normal  Impacts  by  spherical  particles  Is  considered  In  the  published 
reports,  however  analyses  for  oblique  Impacts  of  angular  particles  have 
apparently  been  formulated  but  remain  unpublished  (Mamoun,  1979).  Thus 
only  the  general  relations  for  normal  Impacts  will  be  outlined  here. 

The  approach  used  by  Mamoun  Is  based  on  the  general  assumptions 
that  the  Impact  event  Is  represented  by  a quaslstatlc  loading  condition 
and  that  the  deformations  during  the  collision  are  small.  The  quasl- 
statlc condition  Is  hypothesized  to  be  valid  If  the  Impact  duration  Is 
at  least  four  times  longer  than  the  natural  period  of  the  Impacting 
sphere.  For  the  case  of  a char  or  ash  particle  Impacting  a steel  plate 
the  Impact  velocity  for  which  quaslstatlc  conditions  pertain  Is  on  the 
order  of  600  ms  The  magnitude  of  the  deformations  Is  estimated  from 
the  condition  that  the  distance  q the  two  bodies  approach  each  other 
be  less  than  one-fourth  the  radius  of  the  Impacting  sphere.  The  range 
of  conditions  evaluated  for  the  fully  plastic  case  when  a char  sphere 
Impacts  a steel  plate  shows  this  condition  pertains  for  Impact  veloci- 
ties ranging  from  90  to  245  ms  On  the  basis  of  this  comparison  Mamoun 
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concludes  that  his  approach  Is  valid  for  typical  coal  gasification 
conditions  where  the  flow  velocities  would  not  be  expected  to  exceed 
90  ms 


The  limitation  on  the  Impact  conditions  encompassed  by  the  results 
of  the  model  Is  arbitrarily  based  on  a factor  of  four  and  the  bounds  are 
established  In  terms  of  the  analysis  which  Is  being  evaluated.  The  hy- 
pothesized response  of  the  particle  collision  may  be  somewhat  distinct 
from  the  purely  elastic  properties  assigned  to  the  particle  without  frac- 
ture or  crushing  being  a possibility.  The  Imposed  Impact  velocity 
constraint  requires  more  adequate  substantiation. 

Continuing  with  Mamoun's  straightforward  mechanics  analysis  of  the 
erosion  problem,  the  volume  of  material  removed  for  the  six  material  re- 
sponse regimes  previously  enumerated  will  be  described.  The  multiple 
particle  effects  are  taken  Into  account  in  terms  of  a fatigue  process 
for  all  of  the  cases  considered  except  for  case  6:  single  particle  frac- 
ture of  brittle  materials.  The  volume  removed  is  expressed  In  the  general 


volume  loss/impact  ■ v/N^ 


(3.88) 


where  v denotes  the  elemental  volume  removed  and  is  the  number  of 

impacts  required  to  remove  this  volume.  The  volume  removed  Is  estimated 
on  the  basis  of  assumptions  made  with  regard  to  the  material's  response 
in  a multiple  Impact  environment  and  the  number  is  taken  to  be  the 

fatigue  life  for  the  material  determined  for  uniaxial  stress  loadings. 

Is  thus  prescribed  In  terms  of  seml-empirlcal  fatigue  relationships. 
Some  examples  follow.  The  details  of  Mamoun's  calculations  are  quite 
lengthy  and  parallel  the  work  of  Bitter  In  many  aspects  of  the  Hertzian 
theory;  only  the  general  form  of  the  results  will  be  reported  here. 
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For  case  1,  when  Is  less  than  the  threshold  velocity  for 

yielding  to  occur,  Mamoun  suggests  that  mlcrosllp  takes  place  producing 

fatigue  cracks  which  eventually  reach  a depth  of  0.5  a , where  a 

max  max 

is  the  maximum  contact  radius  for  a sphere  Impacting  an  elastic  half- 
space (Appendix  A).  Mamoun  assumes  the  elemental  volume  removed  is  then 
equal  to  a semlelllpsold,  which  is  equivalent  to  Eq.(B4)  if  a^^^«  1 , 


V ” 


1, 

2 


(3.89) 


and  for  purely  elastic  deformations 


"*(211  ^max  ) 


1/q 


(3.90) 


where  is  the  maximum  pressure  in  the  Hertzian  pressure  profile 

(obtained  from  Eq.(A4)  and  Eq.(A35)  for  spherical  body  impacts),  and 
H and  q are  experimental  fatigue  constants  which  depend  on  the  target 
material.  The  final  expression  for  the  volume  removed  per  Impact  is 
completely  deterministic  when  the  parameters  H and  q are  prescribed 
(refer  to  Table  3.2). 


(^) 

- f(H,v^,V2,E^,E2P^)  ' (3.91) 

where  Is  the  density  of  the  sphere  of  radius  R impacting  at  a 

velocity  , and  Young's  moduli  and  Poisson’s  ratios 

for  the  sphere  and  the  half-space,  respectively. 

The  volume  removed  per  particle  Impact  for  a fully  plastic 
collision  Is  determined  from  the  total  kinetic  energy  imparted  to  the 
target  by  the  Impacting  sphere.  The  sphere  does  not  rebound  from  the 
surface  for  this  Impact  condition.  The  volume  removed  per  Impact,  u^. 

Is  estimated  to  be 


TABLE  3.2.  Constants  Employed  In  Determining 
the  Critical  Number  of  Impacts . 


I 

! 


! I 


Type  of  Material 

CONSTANTS 

q 

H 

s 

I 

AISI  4130  (soft) 

-0.0577 

0.00884 

-0.627 

1.1397 

AISI  304  (hard) 

-0.132 

0.02715 

-0.693 

0.9765 

AISI  4340  (hard) 

-0.0996 

0.02041 

-0.596 

0.7482 

AISI  4340  (annealed) 

-0.0626 

0.00824 

-0.538 

0.6245 

AISI  52100 

-0.0692 

0.02376 

-0.563 

0.2373 

AISI  4130  (hard) 

-0.0529 

0.01258 

-0.693 

1.1155 

AISI  310  (annealed) 

-0.127 

0.01124 

-0.567 

0.8128 

Inconel  X 

-0.138 

0.0239 

-0.708 

1.1679 

Titanium  6A1-4V 

-0.130 

0.0405 

-0.738 

1.5453 

Beryllium 

-0.0655 

0.00314 

-0.353 

0.02541 

1100  Aluminum 

-0.0506 

0.00278 

-0.685 

2.258 

2014-T6  Aluminum 

-0.12 

0.02794 

-0.649 

0.5296 
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(3.92) 


where  is  the  fatigue  life  of  the  target  material  In  terms  of  the 
plastic  strain  amplitude  and  is  the  plastic  flow  stress.  Hamoun 
uses  the  fatigue  relations 
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(3.94) 


where  I,s  are  constants  for  a given  material.  Representative  values 
are  listed  In  Table  3.2. 


According  to  Mamoun's  definition  of  the  plastic  strain  amplitude, 


(3.95) 


and  in  view  of  the  magnitude  of  the  values  for  I listed  in  Table  3.2, 
the  condition  that  Is  only  possible  when  a^^^>R  which  Is  physi- 

cally unrealistic.  Thus,  In  contrast  to  the  expressions  derived  by 
Mamoun  using  Eq,(3.92)  to  (3.94),  it  is  only  necessary  to  consider 
Eq.(3.92)  and  (3.94)  which  would  only  exclude  the  erosion  of  beryllium 
except  for  very  small  deformations.  In  order  for  Eq.(3.95)  to  be  con- 
sistent with  Mamoun's  criterion  for  small  deformations,  ^ 0.14  . 

For  strain-hardening  materials,  Mamoun  finds 

48-1  2(48-1) 

/ , > 8(4+n)  _ 8(4+n) 


where  n is  the  strain-hardening  exponent.  The  expression  In  Eq.  (3.96) 
reduces  to  the  nons train-hardening  condition  (case  2)  when  n~0  . 

A simple  but  Inadequate  modification  of  Eq.  (3.96)  Is  used  for 
treating  elastic-plastic  deformations  (case  4).  This  Is  the  condition 
which  Bitter  Investigated.  Mamoun  also  Investigated  material  removal 
due  to  normal  collisions  by  spherical  bodies  on  brittle  materials  (cases 
5 and  6)  which  will  be  described  In  Section  4. 

The  final  expressions  for  the  volume  removed  per  Impact  In  Eq.  (3.91) 
and  (3.96)  show  that  the  volume  removed  Is  directly  proportional  to  the 
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mass  of  the  Impacting  sphere  (R  dependency)  but  the  velocity  exponent 
Is  a function  of  the  fatigue  characteristics  of  the  target  material. 

The  range  of  velocity  exponents  for  elastic  and  fully-plastic,  non- 
straln-hardenlng  (n»0)  materials  will  be  considered  using  the  data 
provided  In  Table  3.2. 

For  the  elastic  case,  -0.138^q<  - 0.05  . The  velocity  exponent  in 
Eq.  (3.91)  therefore  ranges  from  4.2  to  9.2:  the  lower  value  corresponds 
to  Inconel  X and  the  higher  value  to  1100  aluminum.  These  values  would 
appear  to  be  fairly  high,  although  the  magnitude  of  the  coefficient  ob- 
tained by  evaluating  the  function  f would  have  to  be  determined.  In 
addition,  the  range  of  application  for  these  equations  for  ductile  metals 
is  only  on  the  order  of  10  ms  ^ at  which  point  plasticity  effects  become 
dominant.  For  the  nons train-hardening  fully  plastic  case  the  fatigue 
parameter  s has  a smaller  range,  -0.738  £ s £ -0.538,  excluding  beryl- 
lium. The  velocity  exponent  according  to  Eq.  (3.96)  when  n~0  ranges 
then  from  2.68  to  2.93.  The  Introduction  of  the  work-hardening  parameter 
n would  bring  the  computed  velocity  exponents  within  the  range  of 
values  typically  quoted  for  ductile  metals.  The  comparison  Is  interesting 
but  will  require  more  definitive  confirmation. 


The  analytical  predictions  for  spherical  particles  impacting  normal 
to  the  target's  surface  are  compared  with  hard  angular  particle  Impacts 
(silicon  carbide  and  alumina  particles)  Impacting  normal  but  in  some  cases 
oblique  Impacts  are  Included  in  the  limited  data  base  Mamoun  claimed  was 
available  at  the  time  (Mamoun,  1975d,  1976).  A more  extensive  literature 
survey  would  have  revealed  a somewhat  expanded  amount  of  experimental  re- 
sults. The  materials  for  which  erosion  data  was  reported  are  304  and  310 
stainless  steel  (annealed)  and  1100  aluminum.  The  previously  determined 
volume  losses  are  converted  to  mass  loss  per  mass  of  Impacting  particles. 
However,  Mamoun  notes  that  only  a percentage  of  the  stated  mass  of  im- 
pacting particles  actually  impact  Che  specimen:  this  percentage  ranges 
from  33  to  50%.  Mamoun's  reasoning  is  that  particle-particle  collisions 
occur  as  well  as  a particle  shielding  effect  which  takes  place  at  the 
face  of  the  target  as  Illustrated  by  Moore  (1968).  The  incoming  particles 
collide  with  the  particles  rebounding  from  the  surface  for  near  normal 
Impacts  and  may  therefore  Impact  the  surface  along  a different  trajectory 
and  Impact  velocity  than  assumed  in  the  experiment  or  not  at  all.  How- 
ever Che  specific  determination  of  the  magnitude  of  these  effects  would 
have  to  be  based  on  the  manner  in  which  Che  experiment  was  conducted  and 
the  particle  concentration  in  the  flow  field.  The  comparison  between 
experiment  and  theory  is  quite  tenuous  and  certainly  requires  further 
support  before  any  conclusions  can  be  drawn  regarding  the  merits  of 
Mamoun's  approach. 

3,8  Head  and  Co-workers 


After  completion  of  an  extensive  testing  program  (Head,  et  al., 
1967),  Head  and  Harr  (1970)  attempted  to  develop  a model  applicable  to 
natural  soils  which  Included  consideration  of  the  energy  transferred 
from  the  impinging  particles  to  the  target,  Che  nature  of  the  response 
of  the  target,  and  Che  nature  of  the  erosive  agent  including  pertinent 
descriptors  of  composition,  angularity,  hardness,  and  size  distribution. 
The  result  of  this  effort  was  the  development  of  a weighted  regression 
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analysis  program  (WRAP  model)  In  which 


/ H B \ 

, R,  sina,  cosa,  ^ . g / 


(3.97) 


where  A is  the  erosion  rate  measured  In  volume  removed  per  unit  mass 
of  dust,  V is  the  effective  velocity  of  the  mix,  R is  the  effective 
roundness  of  the  dust  particles,  a is  the  particle  impingement  angle, 

B is  the  hardness  of  the  target,  and  E is  the  erosion  resistance  per 
unit  volume  of  target.  Head  and  Harr  use  the  effective  parameters  which 
they  define  as  the  sum  of  the  parameters  associated  with  the  individual 
size  components  of  the  dust  weighted  on  the  basis  of  the  grain-size  dis- 
tribution of  the  dust.  A number  of  different  materials  were  eroded  with 
four  test  dusts  with  a particle  size  distribution  ranging  from  10  to 
120  Vim  to  supply  the  data  base  for  the  statistical  analysis.  A different 
form  of  interaction  between  the  dependent  variables  was  found  for  brittle 
and  ductile  materials. 

The  erosion  resistance  of  the  target  material  was  measured  by 
simple  expressions  of  the  form 


E - f(a,B,MT) 


where 


a +0 


MT  “ e 


(3.98) 


(3.99) 


is  the  modulus  of  toughness  determined  from  a uniaxial  tension  test. 

The  stress  and  strain  at  rupture  are  denoted  by  and  , respec- 
tively, and  o is  the  yield  point  stress  determined  at  2%  offset. 

The  uniaxial  data  was  obtained  for  the  metallic  specimens  at  a strain 
rate  of  0.05  in/mln  up  to  the  yield  point  and  then  the  strain  rate  was 
increased  to  0.1  in/mln  until  rupture  occurred.  For  brittle  elastic 
materials,  such  as  glass,  the  value  of  MT  is  determined  by  assuming  a 


85 


For  this  case,  It  appears  that  Head  and  Harr 


The  basic  concept  behind  this  approach  which  was  considered  In 
earlier  work  (Flnnle,  et  al.,  1967)  Is  that  the  resistance  of  the  target 
material  to  erosion  was  proportional  to  the  hardness  of  the  material  for 
small  Impingement  angles  and  to  Its  ability  to  absorb  energy  as  repre- 
sented by  the  modulus  of  toughness  at  normal  Impingements.  This  Idea  Is 
most  simply  represented  by  a straight  line  variation 


E 


B - (B-MT) 


(3.101) 


Four  versions  of  the  model  were  formulated  based  on  the  parameters 
appearing  In  Eq,  (3.97):  two  for  brittle  materials  and  two  for  ductile 
materials.  The  model  Is  simply  a statistical  correlation  of  the  test 
data  generated  by  Head  and  his  colleagues  without  any  consideration  for 
the  mechanism  for  material  removal.  The  data  base  for  brittle  materials 
was  not  too  extensive;  the  simpler  step-wise  regression  analysis  Indicated 
the  particle  shape  factor,  R , was  not  significant  and  that  the  erosion 
rate  Increased  as  the  erosion  resistance  of  the  target,  E , Increased. 


The  stepwise  regression  analysis  fitted  the  available  data  for  ductile 
materials  quite  well;  however,  although  the  parameters  appearing  In  the  j 

final  expression  followed  physically  reasonable  trends,  the  erosion  rate  j 

showed  an  even  stronger  correlation  for  Increasing  as  the  target's  ero-  i 


slon  resistance  Increased.  The  resulting  expressions  are 


. v3.06  ^2.69  jj2.08  ^0.03  ^-3.64 


(3.102) 
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for  the  stepwise  regression  analysis  for  brittle  materials,  and 

(3.103) 

for  the  stepwise  regression  analysis  of  ductile  materials.  The  velocity 
exponent  Is  3.06  for  brittle  materials  and  4.34  for  ductile  materials 
which  Is  highly  Inconsistent  with  the  most  of  th«*  erosion  data  from 
other  sources  where  the  velocity  exponent  for  brittle  materials  Is  higher 
than  that  for  ductile  materials  for  which  the  velocity  exponent  typically 
falls  In  the  range  from  2 to  3. 

A weighted  regression  analysis  (WRAP  model)  was  then  carried  out 
for  both  ductile  and  brittle  materials.  The  WRAP  model  Indicated  that 
the  volume  removed  would  decrease  as  the  erosion  resistance  Increased. 

For  brittle  materials  the  volume  loss  Is  weakly  dependent  on  the  Impact 
velocity  with  a velocity  exponent  of  two.  For  ductile  materials  the 
volume  removal  Is  fairly  strongly  dependent  on  the  Impact  velocity  to 
the  fourth  power. 

The  erosion  rates  obtained  previously  for  stainless  steel,  17-4PH 
(AMS5643)  heat  treated  to  a hardness  of  RC38,  using  three  natural  dusts 
and  silica  flour  (Head,  et  al.,  1967)  were  used  as  a basis  for  evaluating 
the  WRAP  model  predictions.  The  particle  shape  factor,  R values  were 
not  available  for  the  three  test  soils:  only  for  silica  flour.  The 
roundness  parameters  for  the  test  soils  were  calculated  from  Eq. (3.103) 
for  each  value  of  a and  the  values  of  R averaged  over  the  number  of 
Impact  angles.  The  effective  particle  velocities  and  effective  particle 
hardnesses  were  computed  and  the  erosion  rates  calculated  from  the  WRAP 
expression  for  ductile  materials.  The  predicted  and  observed  erosion 
rates  agreed  within  a factor  of  two  and  showed  the  same  general  trend 
for  Impingement  angles  between  15  and  75“  for  three  natural  soils,  but 
the  values  at  90“  were  significantly  In  error.  The  predictions  of  the 
WRAP  model  for  silica  flour  using  a properly  evaluated  shape  factor 
showed  the  closest  agreement  with  the  test  data  for  Impingement  angles 
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between  15  and  75”,  however  the  trend  in  the  predicted  and  observed  re- 
sults Is  definitely  reversed  In  terms  of  the  curvature  of  the  two  curves. 
Head  and  Harr  consider  the  agreement  to  be  quite  good  noting  that  the 
predicted  curve  yields  a maximum  erosion  rate  for  an  Impingement  angle 
of  30”  while  the  experimental  data,  according  to  them.  Indicates  the 
maximum  erosion  rate  occurs  at  15”.  This  Is  a surprising  comparison 
since  no  maximum  erosion  rate  has  been  established  for  the  experimental 
data:  the  curve  Is  still  Increasing  monotonlcally  at  15”.  The  data  re- 
ported by  Head  and  Harr  (1970)  for  the  erosion  of  17-4PH  stainless  steel 
by  silica  flour  displays  a similar  trend  as  the  predicted  curve  but  the 
erosion  rates  are  approximately  twice  as  large.  The  erosion  tests  carried 
out  by  Head  and  Harr  (1970)  appear  to  have  been  carried  out  using  a dif- 
ferent sand  blaster  than  used  by  Head,  et  al.  (1967),  although  both  ero- 
sion apparatus  had  an  8 mm  ID  blast  tube.  This  comparison  does  not  eval- 
uate how  well  the  velocity  dependence  Is  predicted  or  If  the  target  mate- 
rial properties  selected,  hardness  and  toughness,  can  adequately  differ- 
entiate the  erosion  response  of  one  material  from  another. 

The  predictive  capability  of  the  WRAP  model  was  subsequently 
evaluated  for  two  ductile  metals  (type  302  stainless  steel  and  6061-T6 
aluminum)  eroded  by  220  mesh  (68  ym)  alumina  particles  and  220  mesh  fluo- 
rite particles  (Head,  et  al.,  1973).  When  the  model  was  adjusted  for  a 
higher  maximum  particle  velocity,  300  ms  ^ as  opposed  to  185  ms  ^ in 
the  earlier  work  Imparted  by  the  carrier  gas  and  changes  In  the  dimen- 
sions of  the  blast  tube,  agreement  between  the  model  predictions  and 
the  experimental  data  for  alumina  particles  Is  claimed  for  Impingement 
angles  between  20”  and  60”  with  good  estimates  of  the  maximum  erosion 
rate;  however,  the  general  trend  of  the  model  predictions  and  the  ex- 
perimental data  as  a function  of  impingement  angle  show  little  resem- 
blance to  each  other.  The  comparison  Is  less  encouraging  for  fluo- 
rite particles  where  the  predicted  maximum  erosion  rates  underestimate 
the  actual  values  by  a factor  of  two  for  the  stainless  steel  and  a factor 
of  four  for  the  aluminum  alloy.  It  Is  Interesting  to  note  that  the  fluo- 
rite particles  with  a Moh  hardness  of  4.0  and  higher  particle  roundness 


than  the  alumina  were  significantly  more  erosive  than  the  alumina  par- 
ticles with  a Moh  hardness  of  9.0.  In  the  case  of  6061-T6  aluminum, 
there  is  nearly  an  order  of  magnitude  difference  in  the  measured  erosion 
rates  for  these  two  particles.  The  WRAP  model  predicts  that  the  erosion 
rate  Increases  with  particle  hardness  and  roundness.  There  does  not  ap- 
pear to  be  a sound  technical  basis  for  the  modified  WRAP  model;  the  re- 
sults are  simply  adjusted  to  match  the  magnitude  of  the  WRAP  predictions 
with  the  experimental  data.  Evidently,  there  are  other  factors  which 
must  be  taken  into  account  with  regard  to  particle  characterization  and 
testing  procedures  which  have  not  been  Identified  by  Head  and  his 
colleagues . 

The  statistical  analyses  developed  by  Head  and  his  co-workers  are 
significant  in  that  this  is  the  only  erosion  modeling  effort  which  ad- 
dresses erosion  due  to  a distribution  of  particle  sizes  rather  than  care- 
fully graded  particulates  generally  used  In  erosion  investigations. 

3.9  Correlations  with  the  Thermal  Properties  of  the  Target  Material 

Several  general  correlations  relating  to  the  thermal  properties  of 
the  target  material  have  been  proposed.  The  validity  of  these  correla- 
tions has  typically  been  determined  by  comparison  to  the  erosion  resis- 
tance for  a number  of  pure  metals  determined  by  Finnie,  Wolak,  and  Kabll 
(1967).  The  original  correlation  with  Vickers  hardness  will  be  described 
and  then  the  thermally-based  correlations  will  be  summarized. 

Finnie 's  theory  of  ductile  erosion  Indicates  that  the  erosion 
resistance  Is  inversely  proportional  to  the  flow  stress,  o , of  the 
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material  [refer  to  Eq.(3.30)  and  (3.31)]  which  is  directly  related  to  a 
measure  of  hardness.  The  experimental  results  of  Finnie,  Wolak,  and 
Kabll  (1967)  for  various  pure  metals  (and  some  steels)  are  reproduced  in 
Figure  3.18.  Note  that  on  this  log-log  plot  tbtire  are  two  distinct  and 
parallel  "lines,"  both  having  a slope  which  represents  a direct  propor- 
tionality between  erosion  resistance  and  Vickers  hardness.  The  one  line. 


Volume  Removed  (inrff  /g  of  abrasive) 


Vickers  Hardness 


Volume  Removed  as  a Function  of  Vickers  Hardness. 
Metals  Eroded  by  250  ym  (60-mesh)  Silicon  Carbide 
Particles  at  a>20*  . All  Metals  except  Cadmium 
were  In  an  Annealed  Condition.  (Data  from  Flnnle 
Wolak,  and  Kabll,  1967). 


which  the  authors  regard  as  the  "main  trend,"  Is  established  by  tin, 
allmlnum,  silver,  copper,  tantalum  and  (with  slight  deviation)  by  cad- 
mium and  nickel.  The  other  "line"  can  be  Imagined  to  pass  through  or 
close  to  the  points  for  bismuth,  magnesium.  Iron  and  molybdenum,  and 
represents  a resistance  about  one-third  that  of  the  "main  trend,"  for 
the  same  hardness.  Tungsten  exhibits  an  even  worse  resistance  relative 
to  hardness,  although  It  is  the  most  resistant  on  an  absolute  level. 
Lead,  on  the  other  hand.  Is  the  most  resistant  relative  to  hardness. 

The  authors  propose  explanations  for  these  results  In  terms  of  the  ml- 
crostructure,  work-hardening  properties  and  dynamic  properties  of  the 
various  materials,  and  make  comparisons  with  corresponding  results  from 
abrasive  wear  tests.  Goodman,  Tilly  and  Sage  (1969)  have  shown  that 
hardness  In  combination  with  other  properties  may  provide  some  correla- 
tion for  the  erosion  resistance  of  alloys. 

Flnnle,  Wolak,  and  Kabll  (1967)  also  Investigated  the  effect  of 
prior  work-hardening  on  the  erosion  resistance  of  several  pure  metals 
and  found  no  effect.  Furthermore,  the  hardening  by  heat-treatment  of 
several  steels  not  only  showed  no  beneficial  effect,  but  actually  re- 
duced the  erosion  resistance  slightly.  The  lack  of  Improved  erosion  re- 
sistance through  coldworking  and  heat  treatment  of  alloys  was  also  found 
by  Smeltzer,  et  al.  (1970),  and  to  a lesser  extent  by  Sheldon  and  Kan- 
here  (1972).  Eq.(3.30)  and  (3.31)  predicts  that  the  annealed  specimens 
with  lower  flow  stresses  should  erode  three  to  four  times  more  readily 
than  the  heat  treated  specimens.  Instead  the  difference  In  erosion 
rates  for  the  annealed  and  heat  treated  specimens  was  less  than  five 
percent  In  all  cases. 

Smeltzer,  et  al.  (1970),  purportedly  observed  evidence  for  target 
material  melting  based  on  microscopic  examinations  of  the  eroded  sur- 
faces of  several  alloys  [2024  Aluminum  (annealed),  410  stainless  steel, 
and  Tl-6Af^4V  (annealed)].  They  provide  a simple  calculation  to  demon- 
strate that  the  kinetic  energy  for  a representative  particle  and  Impact 


91 


condition  is  more  than  sufficient  to  melt  the  volume  of  material  (based 
on  pure  iron)  they  believe  is  removed  by  this  process.  Smeltzer  and  co- 
workers define  two  mechanisms  (Smeltzer,  et  al.,  1970,  p.  216): 

Mechanism  1 (Melting  and  splattering)  - Comer-oriented  particle 
Impacts  are  sufficient  to  cause  localized  target  melting  in  the 
crater  region  and  subsequent  spraying  or  splattering  of  micron- 
size  droplets  of  molten  and  semlmolten  target  metal  into  the  gas 
stream  under  the  high  pressures  generated  by  Impact.  Below  the 
threshold  energy  for  erosion,  there  is  insufficient  energy  avail- 
able for  significant  metal  melting  and/or  splattering.  Metal 
splattering  would  be  most  feasible  (reasonable)  for  the  comer- 
oriented  Impact,  where  the  particle  is  decelerated  over  the 
smallest  total  contact  area,  promoting  most  Intense  heat  buildup, 
maximum  pressurization,  and  maximum  opportunity  for  target 
liquating  and  spraying. 

Mechanism  2 (Melting,  bonding,  and  dlslodgement)  - Corner-oriented 
particle  Impacts  also  result  in  broken  particle  corners  being  em- 
bedded (mechanically  bonded)  in  the  impact  craters  they  formed. 
This  bonding  is  a transient  condition;  most  particles  subsequently 
are  extracted  from  the  craters  with  whatever  residual  energy  they 
possess  or  by  subsequent  particle  Impacts.  But  in  the  extraction 
process,  some  bonded  target  metal  may  adhere  to  the  particle 
surface  and  be  removed. 

Smeltzer  and  co-workers  could  not  establish  the  relative  contribu- 
tions from  each  of  these  mechanisms  to  the  total  erosion  rate  but  in- 
cluded both  mechanisms  in  an  empirical  correlation  based  on  the  idea 
that  the  material  removal  rate  was  proportional  to  the  particle  kinetic 
energy.  They  assume  that  a percentage  of  the  impacting  particles  will 
terminate  their  trajectories  while  in  contact  with  the  surface  of  the 
target  and  thereby  transfer  all  of  their  kinetic  energy  into  localized 
melting  of  the  target.  The  amount  of  material  melted,  W , expressed  in 
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where  Is  the  Impact  velocity,  Q is  the  energy  required  to  heat 

and  melt  one  gram  of  target,  and  <))(a)  is  the  probability  that  a par- 
ticle impacting  the  target  at  an  angle  a will  be  snagged  by  the  target. 
The  energy  Q required  for  melting  one  gram  of  target  material  may  be 
calculated  from  the  equation 


(3.105) 


where  is  the  specific  heat  of  the  target,  is  the  heat  of  fusion 

of  the  target,  and  AT  is  the  difference  between  the  actual  target  tem- 
perature and  its  melting  temperature.  Smeltzer,  et  al.,  assum<i 


♦ (a)  - <K90‘’)sina 


(3.106) 


Since  the  major  portion  of  the  melted  material  is  generated  in  the 
vicinity  of  point  P (Figure  3.19)  on  the  crater-particle  Interface, 
maximum  liquation  is  assumed  to  occur  in  the  vicinity  of  the  comer  of 
the  particle.  A portion  of  the  molten  material  is  squeezed  out  of  the 
vicinity  of  point  P . The  amount  of  material  removed  by  mechanism  1 
is  estimated  to  be 


where  n is  the  fraction  of  splattered  material  which  actually  escapes 
from  the  crater,  and 


is  the  probability  for  splattering  of  the  molten  material  at  point  P . 


mg/g  of  impacting  particles,  is  given  by  the  expression 
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Since  i))(a)  Is  the  probability  for  splattering  of  the  melted 
material,  then  !-<)>  is  the  probability  that  the  melted  material  will 
remain  in  the  vicinity  of  point  P , and  (l-iJ;)W  represents  the  amount 
of  target  material  remaining  at  the  interface.  Then 


W2  - C[l-tp(a)]  W 


(3.109) 


where  Q is  the  fraction  of  the  remaining  melted  material  which  is 
removed  in  accordance  with  mechanism  2. 


The  final  form  of  the  correlation  is  simply  the  sum 


+ W2 
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where 


_ n4»(90»)  , 
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It  is  seen  from  Eq.(3.110)  that  no  material  will  be  removed  when  a>90* 
unless  mechanism  2 is  Introduced. 


The  constants  K^,K2  are  evaluated  from  the  experimental  data  for 
as  a function  of  a . When  a»90” , 
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(3.111) 


Knowing  K2i  K^  is  determined  by  fitting  the  first  term  on  the  right-hand 
side  of  Eq.(3.110)  to  the  difference 
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W^(exp)  - -Y  sin^a 


(3.112) 


using  as  an  adjustable  constant.  The  coefficients  and  are 

evaluated  for  the  alloys  evaluated.  It  is  Interesting  to  note  that  the 
maximum  erosion  rates  in  the  experimental  data  occur  for  30*<a^<40”  in 
contrast  to  the  predictions  of  Finnic 's  model  and  experimental  data  of 
Finnic , Wolak,  and  Kabil,  1967. 


Smeltzer,  et  al.  (1970)  provide  a brief  discussion  of  the  relation 
of  and  to  more  fundamental  properties  of  the  target  material. 
They  Indicate  that,  when  the  erosion  loss  for  pure  metals  [the  erosion 
data  of  Finnic,  Wolak,  and  Kabil  (1967)]  is  plotted  against  such  target 
material  properties  as  Young's  modulus,  absolute  melting  temperature, 
heat  required  for  melting  a unit  mass  Q in  Eq. (3.105),  and  liquid 
metal  surface  tension,  they  all  yield  the  correct  relative  ranking  of 
the  erosion  resistance  but  without  a high  degree  of  absolute  predicta- 
bility. They  suggest  that  erosion  might  occur  simply  from  decohesion 
of  the  target  material. 

Ascarelll  (1971)  Improved  the  correlation  of  erosion  resistance 
and  the  melting  point  temperature  suggested  by  Smeltzer,  et  al.  (1970) 
by  defining  the  thermal  pressure  for  pure  metals.  Ascarelll  notes  that 
according  to  Finnic 's  analysis  the  volume  removed  per  particle  impact 
can  be  %fritten  in  the  form 


u » C 4 ® V ^ f(a) 
2 o 


(3.113) 


where  C is  a constant  dependent  on  the  material  properties  of  the 
target  material,  has  the  dimensions  of  the  inverse  of  an  energy  density 
or  pressure,  and  is  proportional  to  the  volume  removed  per  unit  energy 
of  the  impinging  particles.  In  order  to  determine  the  explicit  material 
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dependence  of  C , Ascarelll  Introduced  the  concept  of  target  melting 
suggested  by  Smeltzer  and  co-vorkers  but  did  not  prescribe  a specific 
mechanism  by  which  material  is  removed  except  that  the  melted  material 
will  be  carried  off  by  the  stream  of  Impinging  particles.  He  observes 
that  while  large  plastic  deformations  do  occur  as  part  of  the  erosion 
process,  they  may  contribute  little  to  the  actual  volume  removal  mech- 
anism. This  is  the  motivation  for  attempting  to  correlate  erosion  re- 
sistance with  the  physical  properties  of  the  target  material  which,  even 
without  a detailed  analysis  of  the  mode  of  material  removed,  may  be 
simply  related  to  such  properties  as  crystal  structure,  compressibility, 
melting  point,  thermal  expansion,  and  elastic  constants. 

These  observations  led  Ascarelll  to  define  the  thermal  pressure 
as  being  representative  of  the  pressure  corresponding  to  the  localized 
region  within  the  crater  being  produced  by  large  plastic  deformations 
which  reaches  the  melting  point  of  the  target  material.  Then 


(3.114) 


where 


a,.(T  -T^) 
_ _ t m o 

c a p_  = 

T X 


Is  the  linear  coefficient  of  expansion 
T^  Is  the  melting  temperature 
T^  Is  the  temperature  of  the  test 
X Is  the  Isothermal  compressibility 


By  Introducing  the  concept  of  thermal  pressure,  Ascarelll  states  that 
Isothermal  compressibility  Is  the  primary  parameter  characterizing  the 
erosion  resistance  of  both  brittle  and  ductile  metals. 

A comparison  Is  made  of  the  volume  removed  determined  by  combining 
Eq.(3.113)  and  (3.114)  with  the  experimental  data  of  Flnnle,  Wolak,  and 
Kabll  (1967)  where  0(~a  . The  experimental  data  falls  within  + 25Z  of 


the  predicted  values  for  volume  r»^val. 

Ascarelll  also  considered  the  form  of  the  function  f(a)  In 
Eq. (3.113)  which  he  relates  to  the  surface  roughness  which  evolves  as 
part  of  the  erosion  process.  The  derivation  of  Eq. (3.107)  Is  general- 
ized accounting  for  surface  roughness.  The  final  form  of  Ascarelll 's 
expression  for  volume  removal  Is 


where 


0.5m  V ‘ 

“ ■ “t'VV'x 


f(0)  . 3ln(2ksln^6)  (2cos^9-l)  ^ sln(k8ln^9)  _ _1 
4ksln^0  kslnO  ^ 


(3.115) 
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The  angle  6 replaces  a as  the  Impingement  angle  and  Is  defined  as 
the  actual  angle  of  Impingement  between  the  particle  and  the  eroded  sur- 
face where  '9-8*  ^ 9 ^ 0+9*  Is  taken  to  be  spread  In  the  values  of  the 
actual  Impingement  angle,  k Is  an  empirical  roughness  parameter  which 
has  to  be  determined  from  the  experimental  data  when  a-90**  which  As- 
carelll Indicates  Is  Inversely  proportional  to  the  ductility  of  the  mate- 
rial. It  Is  seen  that  using  Eq. (3.115)  eliminates  the  need  for  Intro- 
ducing Smeltzer's  mechanism  2,  Eq. (3.109),  which  yields  a non-zero  value 
of  the  material  removed  for  normal  collisions,  since  this  condition  Is 
accommodated  In  Eq. (3.116)  when  k)k)  . When  k>K)  , Eq. (3.109)  results. 

Ascarelll  (1971)  suggests  that  the  temperature  dependence  of 
erosion  results  from  the  competition  of  two  contrasting  effects.  The 
first  Is  that  the  ductility  of  metals  Increases  with  Increasing  tempera- 
ture which  leads  to  larger  deformations  (for  constant  Impinging  particle 
velocity)  and  reduces  the  probability  of  energy  localization.  However, 
this  Is  contrasted  by  the  reduced  energy  necessary  to  get  a region  of 
the  metal  to  the  melting  point.  He  Indicates  that  due  to  the  difficulty 
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In  predicting  or  describing  the  ductility  of  metals.  It  Is  Indeed  a 
difficult  problem  to  make  predictions  for  the  temperature  dependence  of 
erosion  for  temperature  larger  than  say  0.5  . According  to  Ascarelll 

temperature  dependent  erosion  behavior  should  be  similar  to  the  tempera- 
ture dependent  behavior  of  the  shear  modulus  or  of  Young's  modulus. 

Brauer  and  Krlegel  (1965)  and  Truitt  (1974)  found  some  correspondence 
between  erosion  resistance  and  Young's  modulus  for  certain  metals  and 
alloys. 

The  experimental  work  of  Jennings,  et  al.  (1976),  supported  the 
possibility  for  target  melting  previously  reported  by  Smeltzer,  et  al. 
(1970).  Following  a procedure  similar  to  that  used  by  Head  and  Harr 
(1970) , they  obtained  the  following  correlation  between  the  target  volume 
loss  per  unit  mass  of  Impacting  particles,  A , and  the  properties  of 
the  target  material. 
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where 


Is  the  kinetic  energy  transferred  from  the  impacting 
particles  to  the  target  per  unit  mass  of  particles 

R Is  the  particle  roundness  parameter 

G Is  the  gram  molecular  weight  of  the  target 

p Is  the  density  of  the  target 

k Is  the  thermal  conductivity  of  the  target 

T^  Is  the  melting  temperature  of  the  target 

AH^  Is  the  enthalpy  of  melting  for  the  target 

S Is  a dimensionless  scaling  factor  coupling  the 

calculated  values  with  the  erosion  facility. 


Comparisons  with  the  experimental  data  are  provided  but  will  not  be 
considered  here  . 

The  kinetic  energy  tr^-^sf erred  expressed  as 
V ^ 

- -f-  (1-e^)  , (3.118) 

where  e Is  the  coefficient  of  restitution,  was  evaluated  by  examining 
photographic  records  of  the  Incoming  and  rebounding  particles.  The 
average  conditions  were  estimated,  however  the  detailed  statistical 
analyses  undertaken  by  Tabakoff  and  his  co-workers  (Grant  and  Tabakoff, 
1973;  Ball  and  Tabakoff,  1974)  were  not  pursued.  A constant  velocity 
exponent  of  five  results  when  Eq. (3.118)  Is  substituted  Into  Eq.(3.117) 
which  Is  considerably  In  excess  of  any  of  the  experimental  values  re- 
ported for  ductile  metals.  This  observation  should  be  kept  In  mind  In 
relation  to  the  sensitivity  of  the  model  verification  provided  by 
Jennings,  et  al.  (1976). 

In  contrast  to  Smeltzer,  et  al.  (1970),  Jennings  et  al.  (1976) 
acknowledge  that  material  removal  is  occurring  by  both  mechanically  and 
thermally  activated  processes,  although  they  cannot  determine  the  portion 
of  the  volume  removed  which  should  be  contributed  to  each  of  these 
processes.  In  addition  to  the  evidence  they  claim  Is  observed  In  sup- 
port of  material  melting  based  on  microscopic  examination  of  eroded 
surfaces,  they  note  that  the  erosion  rates  for  ductile  materials  do  not 
change  appreciably  beyond  a material-dependent  level  of  the  kinetic 
energy.  After  this  condition  Is  reached  the  Impact  velocity  as  well  as 
the  particle  shape  are  Ineffectual  In  modifying  the  volume  loss. 

Hutchings  (1975)  finds  that  the  data  of  both  Finnic,  Wolak,  and 
Kabll  (1967)  and  that  of  Tadolder  (1966)  for  pure  metals  can  be  cor- 
related with  the  product 
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Cp(T  -T  ) 
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where  C Is  the  specific  heat  of  the  material  at  room  temperature 

p is  the  density  of  the  target  material. 

Hutchings  states  that  the  one  set  of  data  (Flnnle,  Wolak,  and  Kabil,  1967) 
can  be  fitted  within  the  limits  indicated  by  Ascarelll  in  terms  of  the 
thermal  pressure,  Eq. (3.114),  but  that  comparable  agreement  is  not  found 
for  Tadolder's  (1966)  results.  The  relation  in  Eq.( 3.119)  compares  more 
favorably  with  the  expanded  data  base.  For  pure  metals  Vijh  (1976) 
demonstrated  that  the  same  erosion  data  correlated  with  the  inverse  of 
the  atomic  bonding  energies  representative  of  the  cohesion  of  the  metals 
tested.  These  correlations  are  for  pure  metals;  modified  relations  would 
be  required  for  metal  alloys. 

Jennings,  et  al.  (1976),  among  others,  point  out  the  general 
observation  in  solid  particle  erosion  that  the  macro-mechanical  proper- 
ties of  the  target  do  not  significantly  affect  the  erosion  behavior  of 
ductile  materials.  However,  in  view  of  the  above,  there  does  not  seem 
to  be  a shortage  of  thermal  and  physical  material  properties  which  show 
a correspondence  with  the  erosion  data. 


4.0  EROSION  CORRELATIONS  AND  AMALYSES;  BRITTLE  MATERIALS 

! 

s 

In  recent  years  Che  fracture  morphologies  in  ceramics  associated 
with  Indentation  and  Impact  by  both  spherical  and  angular  particles  have 
been  Investigated  In  some  depth.  The  current  view  of  particle  Impinge- 
ment on  brittle  ceramics  divides  the  fractures  produced  Into  two  dis- 
tinct categories:  the  elastic  response  regime  and  the  elastic-plastic 
response  regime. 

The  elastic  response  regime  Is  characterized  by  segmented 
circumferential  cracks  surrounding  the  Impact  site  analogous  to  the 
Hertzian  cone  cracks  In  glass  when  Impacted  or  Indented  with  a spherical 
body.  Investigations  of  these  characteristic  cone-shaped  cracks  In  glass 
are  fairly  extensive  (for  example.  Tlllett,  1956;  Benbow.  I960;  Tsai  and 
Kolsky,  1967;  Oh  and  Flnnle,  1967;  Fisher,  1967;  Frank  and  Lawn,  1967; 

Hamilton  and  Rawson,  1970;  Evans,  1973;  Knight,  Swain,  and  Chandhrl, 

1977;  Chaudhrl  and  Walley,  1978;  Klrchner  and  Gruver,  1978). 

The  fractures  characteristic  of  the  elastic-plastic  response  regime 
are  typically  produced  by  pointed  Indenters  or  angular  particle  Impacts. 

Evans  and  Wllshaw  (1976)  Indicate  that  for  spherical  Indenters  of  small 
radii  the  elastic-plastic  regime  applies  and  the  hardness  of  the  target 
Is  the  controlling  material  parameter.  For  radii  above  a critical  value 
the  Hertzian  fracture  analysis  Is  appropriate  and  the  controlling  mate- 
rial properties  are  the  fracture  toughness  and  the  surface  flaw  size 
distribution.  The  fracture  sequences  are  material  dependent  and  have 
been  described  in  a number  of  papers  (Tsai  and  Kolsky,  1967;  Lawn  and 
Swain,  1975;  Lawn,  Swain,  and  Phillips,  1975;  Lawn  and  Fuller,  1975; 

Swain  and  Hagan,  1976;  Evans  and  Wllshaw,  1976,  1977;  Knight,  Swain  and 
Chaudhrl,  1977;  Perrott,  1977).  Fracture  development  due  to  a small 
spherical  Indenter  on  a seml-brlttle  ceramic  Is  shown  schemtlcally  In 
Figure  4.1.  Lawn  and  co-workers  (1975)  do  not  observe  the  radial  cracks 
as  the  Initial  fractures  In  the  brittle  materials  they  Investigated  when 
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loaded  with  a pointed  Indenter;  their  observations  were  made  on  silicon, 
quartz,  fused  silica,  and  soda  lime  glass.  As  the  sharp  Indenter  pene- 
trates the  target  It  Is  thought  to  produce  a localized  elastic-plastic 
zone  surrounding  the  Indentation  site.  It  Is  the  modem  view  of  Inden- 
tations In  ceramics  that  localized  plastic  deformations,  even  In  essen- 
tially brittle  materials,  play  a critical  role  In  the  crack  formation 
process.  During  loading  a penny-shaped  crack  forms  below  the  pointed 
Indenter,  assuming  a normal,  axlsymmetrlc  loading  condition.  During  the 
unloading  phase  the  penny-shaped  crack,  referred  to  as  a median  vent 
fracture  will  reseal,  and  lateral  cracks  will  propagate  outward  from  the 
Impact  site  due  to  residual  stress  effects.  These  cracks  Initiate  In 
the  vicinity  of  the  base  of  the  penetration  zone,  extend  outward  nearly 
parallel  to  the  face  of  the  specimen,  and  then  tend  toward  the  surface. 
These  are  referred  to  as  lateral  vent  fractures. 

The  fracture  sequence  shown  In  Figure  4.1  Is  characteristic  of 
Indentations  with  a small  spherical  Indenter  In  chemically  vapor  de- 
posited (CVD)  zinc  sulfide,  hot-pressed  (HP)  silicon  nitride,  single 
crystal  spinel,  and  single  crystal  sapphire  (Evans  and  Hllshaw,  1976). 

The  observations  by  Evans  and  Wllshaw  (1976)  for  glass  were  consistent 
with  those  of  Lawn  and  Swain  (1975)  In  that  radial  and  lateral  vent 
cracks  did  not  form  on  loading  and  on  unloading  the  median  vent  cracks 
propagated  to  the  surface  and  at  near  zero  load  lateral  cracks  developed. 
Klrchner  and  Gruver  (1978)  demonstrated  the  transitions  In  fracture  modes 
for  glass  plates  Impacted  by  3 mm  diameter  glass  beads  by  heating  the 
specimens  to  770*C.  Hertzian  cone  cracks  were  observed  at  temperatures 
up  to  about  700‘C.  Radial  cracks  were  also  found:  the  number  of  radial 
cracks  Increased  with  Increases  In  temperature  and  particle  Impact  veloc- 
ity. Returning  to  the  experiments  of  Evans  and  Wllshaw  (1976),  It  Is 
noted  that  radial  and  median  vent  cracks  In  CVD  zinc  selenlde  occurred 
In  accordance  with  the  sequence  shown  In  Figure  4.1,  but  lateral  vent 
fractures  would  only  result  through  spatial  Interaction  of  serveral 
adjacent  Indents. 


Analyses  of  the  fractures  described  here  have  been  formulated  by 
a number  of  authors  as  previously  cited.  For  the  elastic  response  regime 
the  Hertslan  cone  cracks  have  been  considered  as  the  basis  for  erosion 
theories  of  brittle  materials  by  Sheldon  and  Flnnle  (1966b),  Mamoun  (1975a), 
and  Adler  (1976a, b).  In  view  of  the  current  understanding  of  the  elastic- 
plastic  response  regime  for  ceramics,  the  median  vent  fractures  and 
radial  cracks  are  considered  to  contribute  to  general  strength  degrada- 
tion of  the  ceramic  component  when  subjected  to  quaslstatlc  loadings, 
while  the  lateral  vent  fractures  contribute  to  the  erosion  of  the  exposed 
surface. 

A.l  Elastic  Response  Regime 

Bitter  (1963a)  Indicated  his  deformation  wear  analysis  (Section  3.4) 
developed  for  ductile  materials  was  also  applicable  to  brittle  materials, 
such  as,  cement  and  glass.  According  to  Bitter  these  materials  have  no 
capacity  to  absorb  plastic  energy,  so  presumably  the  ultimate  embrittle- 
ment of  work  hardening  ductile  materials  Is  present  at  the  Initial  Impact. 

At  "not  too  high  velocities"  the  crack  pattern  In  glass  Is  taken  to  be 
the  same  as  In  embrittled  ductile  materials.  Bitter  concludes  that  the 
amount  of  material  removed  from  a glass  surface  Is  proportional  to  the 
diameter  of  the  Impacting  spherical  body  and  Inversely  proportional  to 
the  one-dlmenslonal  strain  energy  of  the  target  (as  obtained  from  a uni- 
axial tension  test).  The  maximum  principal  stress  difference  employed 
by  Bitter  as  a fracture  criterion  Is  a compressions!  stress  state  which 
Is  a secondary  failure  mode  In  glass  which  Is  more  susceptible  to  tensile 
^*^7'Ures.  The  maximum  tensile  stresses  for  an  elastic  sphere  Indenting 
an  elastic  half-space  occurs  at  the  periphery  of  the  contact  radius 
[Eq.(A33)].  On  the  basis  of  the  Hertzian  stress  distribution  summarized 
in  Appendix  A,  it  Is  seen  that  the  critical  stress  difference  defined  by 
®l^***>f  oot  proportional  to  the  radius  of  the  Impacting  sphere  but  Is 
a more  complex  function  of  the  sphere  radius.  Even  along  the  axis  of 
the  dependence  on  the  critical  stress  difference  would  not  be 
proportional  to  the  radius  of  the  sphere. 
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A theoretical  analysis  for  the  erosion  of  brittle  materials  was 
presented  by  Sheldon  ^uld  Flnnie  (1966b).  Whereas  Bitter's  deformation 
wear  analysis  is  based  on  an  elastic-plastic  model,  Sheldon  and  Flnnie 's 
analysis  is  based  on  an  elastic-brittle  model.  The  penetration  of  the 
particle  into  the  material  specimen  is  calculated  from  the  Hertzian 
theory  of  Impact  and  a region  of  cone  cracks  is  assumed  to  spread  around 
the  Impact  point. 

Sheldon  and  Flnnie  apply  the  Hertzian  theory  for  the  normal 
collision  of  a rigid  sphere  on  an  elastic  half-space  to  irregularly- 
shaped  particles  by  introducing  an  Idealized  sphere  of  radius  R*  as 
shown  in  Figure  4.2.  The  radius  R is  that  of  a sphere  of  equal  mass 
to  that  of  the  particle,  while  R'  is  the  radius  which  is  Involved  in 
causing  fracture.  As  the  indentation  process  continues  material  which 
is  already  cracked  will  transmit  load  to  as  yet  uncracked  material 
causing  further  fracture.  When  the  particle  finally  comes  to  rest,  the 
condition  depicted  in  Figure  4.2  is  the  basis  for  the  model.  The  new 
effective  radius  of  the  particle,  R , for  determining  Che  volume  re- 
moval is  then  estimated  to  be  the  radius  of  a sphere  which  just  produces 
cracking  when  indenting  to  a depth  . The  volume  removed  per  particle 
u , is  taken  proportional  to  the  volume  of  the  spherical  cap  defined  by 
R*,  a*  , and  q^ 


(4.1) 


According  to  Sheldon  and  Flnnie  (1966b)  the  value  of  the  coefficient 
should  be  on  the  order  of  unity  and  can  be  estimated  from  measurements 

of  volume  removal  per  particle  In  erosion  tests  and  calculations  of 

* 

a and  q . 


Flnnie  and  co-workers  (Oh,  et  al.,  1972)  subsequently  defined  the 
volume  removed  per  particle  as 


(4.2) 
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using  Eq.(BS)  which  includes  the  complete  cracked  zone,  however  they  then 
let  S be  proportional  to  , so  Eq.(4.1)  once  again  represented  the 
volume. 


The  erosion  of  brittle  materials  Is  based  on  the  Idea  that  the 
volume  removed  per  spherical  bead  Impact,  u , Is  proportional  to  the 
product  of  the  statistical  average  of  the  square  of  the  radius  of  the 
outermost  ring  crack  (a  ) and  the  depth  of  the  Inner  ring  fracture. 
Using  the  statistical  flaw  theory  of  Oh  and  Finnle  (1967)  and  the  usual 
Hertzian  relations,  the  primary  result  of  the  analysis  Is  that 

U - C R®  (4.3) 

where  for  spherical  particles 

a - (3m-2)/(m-2),  b » 0.8(3t»-2) /(m-2)  (4.4) 

and  for  angular  particles 

a - 1.2(3m-2)/(m-2),  b - 0.8(3m-2)/(m-2)  (4.5) 

The  parameter  m Is  the  exponent  for  a Welbull  distribution  which  Is 
assumed  to  be  constant  for  a given  material.  It  Is  seen  that  the  veloc- 
ity exponent  is  the  same  for  both  spherical  and  angular  particles.  The 
experimental  data  for  hardened  steel,  graphite,  alumina,  and  magnesia 
shows  the  velocity  exponent  to  range  from  2.53  to  2.74  for  SIC  grit  in 
agreement  with  the  theory  and  identical  values  for  graphite  and  magnesia 
^■pocted  by  steal  shot.  The  experimental  value  of  the  velocity  exponent 
for  steel  shot  Impacting  plate  glass  deviates  significantly  from  the 
calculated  value. 
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Ruff  and  Wiederhorn  (1979)  note  that 


"although  the  theory  by  Sheldon  and  Flnnle  provides  a reason- 
able description  in  brittle  materials,  its  physical  basis  must 
be  questioned  because  it  assumes  Hertzian  crack  formation, 
whereas  lateral  crack  formation  is  the  main  cause  of  material 
removal  during  erosion." 

This  is  a fairly  subjective  view.  The  Hertzian  crack  formations  are 
responsible  for  material  removal  from  brittle  material  surfaces  for 
spherical  particles  for  particular  ranges  of  bead  diameters,  materials, 
and  Impact  velocities.  Flnnle  (1970)  noted  that 

"in  our  treatment  of  erosion  due  to  perpendicular  Impingement 
by  solid  particles,  we  took  the  volume  removed  by  a single 
particle  as  proportional  to  the  cracked  volume  produced  by  a 
single  particle  striking  a smooth  surface.  Thus,  at  least 
one  erosion  test  must  be  made  on  the  chosen  material  before 
quantitative  predictions  can  be  made.  We  took  this  approach 
because  the  intersection  of  the  fracture  patterns  from  many 
particles  will  be  a complex  process,  and  more  work  needs  to 
be  done  on  this  topic  before  a complete  quantitative  treat- 
ment of  erosion  is  possible." 

Such  a program  was  pursued  by  Adler  (1974a, b;  1976b). 


Adler  (1976b)  investigated  the  progression  of  material  removal 
from  the  Initial  single  particle  Impacts,  to  a few  interacting  Impacts, 
and  then  to  significant  levels  of  mass  removal.  In  the  usual  solid  par- 
ticle erosion  tests  the  incubation  period  for  mass  loss  (described  in 
Section  2)  occurs  so  rapidly  that  the  interaction  of  the  impacts  on  the 
starting  surface  go  unnoticed.  It  is  only  when  the  particle  concentra- 
tions are  low  and  short  incremental  exposures  are  used  that  the  manner 
in  which  the  initial  surface  is  obscured  can  be  observed.  The  progression 
of  erosion  mechanisms  operative  on  flat  surfaces  of  borosllicata  glass 
specimens  (Coming  Glass  Works  designation  Pyrex  7740)  and  fused  silica 
specimens  (Coming  7940)  Impacted  by  70,  135,  and  290  micron  glass  beads 
were  established  through  detailed  microscopic  investigation  of  a number 
of  progressively  eroded  surfaces. 
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The  objective  for  this  research  was  the  development  of  an  analyti- 
cal model  for  the  overall  particle  erosion  process.  The  model  was  based  I 

on  a concept  of  preferential  erosion  pit  nucleatlon  and  growth  (Adler, 

1976a)  which  is  more  characteristic  of  liquid  drop  impingement,  however,  ^ 

the  initial  effort  simplified  the  erosive  medium  by  employing  glass  beads 
and  reduced  the  complexity  of  the  response  of  the  target  material  by  con- 
sidering elastic-brittle  materials.  The  functional  requirements  in  the 
general  theory  were  explicitly  evaluated  for  this  special  case.  The 
complexity  of  pursuing  a detailed,  mechanism-oriented  model  became  ap- 
parent (Adler,  1976b).  This  effort  represents  the  only  known  attempt 
in  the  solid  particle  literature  to  mechanistically  relate  the  initial 
single  particle  Impact  event  with  the  mode  of  removal  during  the  steady- 
state  period  for  material  removal. 

i 

Mamoun's  formulation  for  material  removal  due  to  normal  collisions 
of  spherical  bodies  on  brittle  materials  will  now  be  considered  (Mamoun, 

1975a).  As  noted  in  Section  3.7  Impacts  above  and  below  the  fracture 
threshold  are  treated  separately  (cases  5 and  6).  For  case  5,  Mamoun 
views  crack  nucleatlon  and  propagation  as  the  result  of  a mlcrosllp  proc- 
ess and  therefore  uses  the  same  formulation  for  the  material  removed  as 
for  case  1 (Section  3.7).  No  information  is  provided  on  how  the  elastic 
fatigue  constants  are  to  be  obtained  for  brittle  materials,  such  as,  in- 
organic glasses  or  alumina. 

A different  point  of  view  is  adopted  when  the  Impact  velocity  is 
sufficient  to  produce  brittle  fractures  which  are  taken  to  be  in  the  form 
of  Hertzian  cone  cracks.  Mamoun  comments  that  the  crack  depth  can  be  j 

determined  for  statically  applied  loads  using  fracture  criteria  from  frac-  | 

ture  mechanics  however  acceptable  fracture  criteria  are  not  available 
for  predicting  the  extent  of  the  fracture  due  to  an  Impacting  sphere.  j 

Mamoun  therefore  resorted  to  an  energy  relation  by  assuming  that  the  { 

energy  difference  between  the  initial  kinetic  energy  with  which  the  1 

sphere  Impacts  the  target  and  that  which  it  has  when  it  leaves  the  target 

i 

i 
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surface  is  equal  Co  the  energy  expended  in  creating  the  Hertzian  fracture 
surface.  The  Hertzian  fracture  is  in  Che  form  of  a truncated  cone  which 
increases  in  radius  with  depth.  These  fracture  surfaces  have  been  studied 
extensively  in  glass  (as  already  indicated  in  the  references  provided) . 

The  depth  of  penetration,  q , of  Che  cone  crack  is  determined  from  Che 
ratio  of  the  energy  supplied  to  that  required  to  form  a cylindrical  frac- 
ture surface  with  a radius  equal  to  Che  Hertzian  contact  radius. 


where  e is  Che  coefficient  of  restitution  and  y is  the  surface  energy 
per  unit  surface  area.  The  volume  removed  is  assumed  to  equal  Co  Che 
cylindrical  volume  described  above,  so  the  volume  loss  per  Impact  is 


Using  Eq.(4.6)  end  the  expression  for  a^ 
of  Impact 


based  on  the  Hertzian  theory 


1280it^pJ 


For  this  case  is  not  prescribed,  so  the  above  result  is  incomplete. 
The  form  of  the  mass  removal  expression  subsequently  presented  by  Mamoun 
(197Sd)  differs  from  the  form  in  Eq.(4.8),  however  only  a brief  descrip- 
tion is  provided  concerning  the  derivation  of  Che  revised  form. 


Maaoun  (1975d)  compared  his  analysis  with  what  he  claimed  was  the 
axlsclag  sspsrlmsntal  data  which  consists  of  a single  data  point:  127  ym 
silica  carbide  particles  impacting  hlgh-purlty  alumina  at  152  ms  ^ (Shel- 
don and  Finals,  1966a).  Although  Sheldon  and  Flnnle  (1966a)  provide  ero- 
sion data  for  glass,  graphite,  and  magnesia,  Hamoun  Ignores  Its  existence 


A. 2 Elastic-Plastic  Response  Regime 


After  fairly  extensive  Investigations  of  quaslstatlc  spherical 
bead  Indentations  and  particle  Impacts  by  single  spherical  beads  at 
normal  Incidence  on  ceramics,  Evans  and  co-workers  (Evans  and  Wllshaw, 
1976,  1977;  Evans,  et  al.,  1978)  formulated  two  seml-emplrlcal  erosion 
models  based  on  seml-emplrlcal  relations  between  the  magnitude  of  the 
lateral  fractures  produced  and  the  relevant  target  and  projectile  param- 
eters. The  seml-emplrlcal  relations  are  obtained  from  an  approximate 
stress  analysis  of  the  Indentation  process  and  fracture  mechanics  (quasl- 
statlc erosion  analysis)  or  from  an  approximate  analysis  of  the  Impact 
dynamics  and  fracture  mechanics  (dynamic  erosion  analysis). 


The  upper  bound  for  material  removal  per  particle  Impact  is 
determined  from  the  cylindrical  volume  encompassed  by  the  radial  extent 
of  the  lateral  cracks,  , and  their  maximum  depth,  q . Thus, 


IT  Cj^  q 


(A. 9) 


For  the  quaslstatlc  Impact  condition,  Evans  and  Wllshaw  (1976)  find  that 
U « (A. 10) 


Evans  observes  that  little  Interaction  occurs  between  lateral  fractures 
from  adjacent  Impacts,  so  a first  order  approximation  to  the  material 
removal  process  would  be  a simple  summation  of  the  volume  removed  per 
particle.  The  total  number  of  impacting  particles  Is  N,G^  Is  the 
shear  modulus  for  the  target  material,  and  p Is  the  density  of  the 
Impacting  particle.  The  dynamic  fracture  toughness  and  dynamic  hardness, 
Kq  and  Hjj  , are  introduced  without  conmient. 
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Evans  and  Ullshaw  apparently  felt  that  the  values  for  dynamic 
loadings  should  replace  the  static  values  In  their  analysis  but  do  not 
Indicate  how  they  might  be  evaluated  or  provide  justification  for  t'  > 
modification.  (Compare  Evans  and  Wllshaw,  1977.) 

The  derivation  of  Eq.(4.10)  Is  obtained  from  an  experimental  data 
base  for  spherical  Indenters  with  a liberal  Interchange  of  the  loads  and 
Indent  dimensions  for  a Vickers  Indenter  with  the  loads  and  contact  zone 
radii  for  a spherical  bead  Indentation.  The  crack  dimensions  are  related 
to  the  fracture  toughness,  however  the  depth  of  the  lateral  fractures  Is 
conjectured  to  be  proportional  to  the  radius  of  the  contact  zone  which 
Is  the  then  equated  to  the  length  of  the  hardness  Impression:  all  without 
experimental  verification.  The  load  on  the  spherical  indenter  and  the 
Vickers  Indenter  are  also  considered  equivalent  and  for  the  low-velocity 
Impact  condition  the  load  Is  expressed  In  terms  of  the  Impact  velocity 
according  to  the  Hertzian  theory  of  Impact  [Eq.(A3S)  In  Appendix  A]  as 
Is  evident  In  the  form  of  the  parameters  other  than  and  ap- 

pearing In  Eq.(4.10).  The  derivation,  although  lacking  In  rigor,  intro- 
duces an  explicit  dependence  on  the  fracture  toughness  and  hardness  Into 
the  expression  for  material  removal  which  Is  Intuitively  appealing. 


The  basis  for  the  dynamic  analysis  Is  somewhat  more  refined:  the 
relation  between  q and  the  contact  zone  Is  experimentally  defined  and 
the  characterization  procedures,  experimental  damage  measurements,  and 
analytical  procedures  have  expanded  (Evans,  et  al.,  1978).  Based  on 
the  Implicit  assumption  that  the  lateral  crack  dimensions  are  propor- 
tional to  the  radial  crack  sizes  for  which  an  expression  was  derived, 
the  Input  to  Eq.(4.9)  Is  now 


(p  p G G ) 

C.ocC  oc  ' P.  ^ P.  5/ 


1/3 


4/3 


(4.11) 
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Eq.(4.9)  becomes 


( p p G G 

oc  V P t P t) 


19/12  11/3  19/6  -4/3  -1/4 


K H 
c 


(4.13 


as  given  by  Evans  (1979). 

The  quaslstatlc  elastic-plastic  Indentation  experiments  indicate 
that  the  lateral  cracks  propagate  outward  from  the  region  of  the  plastic 
zone  (Lawn  and  Swain,  1975).  This  observation  prompted  specification  of 
an  alternative  form  of  the  material  removal  volume  in  Eq.(4.9),  as  re- 
ported by  Ruff  and  Wlederhorn  (1979).  Eq.(4.11)  is  used,  however  the 

depth  of  the  lateral  cracks  is  taken  to  be  proportional  to  the  depth  of 
the  maximum  particle  penetration.  Eq.(4.12)  is  replaced  by 


(p  V ^ h"^) 
p o 


ll/9j^ll/3^  22/9  -4/3  1/9 

^P  o c 


(4.14; 


(4.15 


The  significance  for  this  modified  form  is  not  clear. 

Note  that  Ruff  and  Wlederhorn  (1979)  comment  that  Evans,  et  al.,  (1978) 
found  a weak  experimental  dependence  of  on  Pp  and  eliminated  Pp 

from  Eq.(4.11).  However  since  pp  enters  the  theory  again  in  establish- 
ing the  erosion  rate,  the  dependence  on  Pp  is  included.  This  is  the 
reason  for  the  difference  in  the  form  of  Eq.(4.13)  as  found  in  Evans, 
et  al.  (1978),  compared  with  Ruff  and  Wlederhorn  (1979)  and  Evans  (1979) 


The  relation  between  the  volunie  removed  per  particle  and  the 
4/3  1/4 

parameters  (K^  H ) was  compared  with  Gulden's  data  (Gulden,  1979a) 
for  four  ceramics  Impacted  by  115  pm  quartz  particles  at  from  98  to 
180  ms  ^ to  show  that  the  correlation  was  nonlinear  (Evans,  et  al.,  1978). 
The  Inclusion  of  hot-pressed  silicon  nitride  may  not  be  a valid  compari- 
son, since  the  Impact  conditions  were  below  the  lateral  fracture  thres- 
hold. Thus,  while  material  was  being  removed.  It  was  not  being  removed 
In  accordance  with  the  mechanism  upon  which  Eq.(4.13)  Is  based.  The 
correlation  with  the  remaining  data  points  appears  reasonable.  On  the 
other  hand,  the  relation  U“(H^^^/K^^)^^^  In  Eq.(4.1S)  shows  a much 
stronger  nonlinear  correlation  of  the  experimental  values.  Thus,  Eq.(4.13) 
appears  to  be  a better  form  for  representing  the  material  properties. 

The  particle  velocity  and  size  dependence  Indicated  In  Eq.(4.13) 

are  In  excellent  agreement  for  magnesium  fluoride  Impacted  by  quartz 

particles  ranging  In  the  size  from  10  to  385  pm  over  a velocity  range 

from  50  to  285  ms  ^ (Gulden,  1979b) . Reaction-bonded  silicon  nitride 

4 4 

also  showed  good  correlation  with  R for  quartz  particles  when  the 

Impact  condition  was  above  the  threshold  for  lateral  crack  formation, 

however  the  volume  removed  per  particle  Impact  for  glass  bonded  alumina 

and  hot-pressed  silicon  nitride  Impacted  by  quartz  particles  showed 

3 

very  good  agreement  with  a R (the  particle  momentum)  correlating 

factor  (Gulden,  1979a).  When  hot-pressed  silicon  nitride  was  Impacted 

by  10  to  940  pm  silicon  carbide  particles,  the  data  was  consistent  with 
4 4 

an  R correlating  parameter  which  Indicates  the  significance  of 

the  particle  target  Interaction.  The  determination  of  the  range  of 
application  of  the  correlation  In  Eq.(4.13)  certainly  requires  further 
consideration. 
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Due  to  the  similarities  In  the  nature  of  the  fractures  formed  for 
quasistatic  and  dynamic  loadings  In  the  elastlc'-plastlc  response  regime, 
there  Is  some  motivation  to  use  the  quaslstatlc  fracture  measurements. 
Evans  has  maintained  an  awareness  of  the  transient  nature  of  the  particle 
Impact  event.  This  approach  appears  to  be  warranted  In  that  the  experi- 
mental data  shows  that  the  radial  and  lateral  cracks  lengths  for  zinc 
sulfide  and  magnesium  fluoride  normalized  with  respect  to  the  Impression 
radius  for  spherical  beads  are  significantly  larger  for  the  dynamic  load- 
ing conditions  compared  with  the  quaslstatlc  loadings.  The  difference 
ranges  from  a factor  of  two  to  four  for  the  larger  plastic  Indentation 
radii.  Knight,  Swain,  and  Chaudhrl's  (1977)  work  on  glass  also  indicates 
that  a quaslstatlc  analysis  of  Impact  damage  significantly  underestimates 
the  extent  of  the  damage. 
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5.0  DISCUSSION 


Both  the  experimental  and  analytical  approaches  to  solid  particle 
erosion  phenomena  as  found  In  the  publications  dating  from  the  early 
1930*8  to  the  present  time  were  reviewed.  Some  of  the  deficiencies  In 
solid  particle  erosion  testing  will  be  Indicated,  although  this  subject 
has  not  been  adequately  treated  In  this  report.  Then  some  of  the  ob- 
servations regarding  solid  particle  erosion  modeling  will  be  sunmarlzed 


5. 1 Solid  Particle  Erosion  Testlni 


The  experimental  work  has  revealed  general  trends  In  the  erosive 
response  of  materials  exposed  to  multiple  solid  particle  Impacts,  but 
an  accurate  and  reproducible  quantitative  evaluation  of  the  erosion  of 
materials  still  has  not  been  achieved.  The  lack  of  a reliable  data 
base  Is  due  to  the  need  for  a test  facility  which  will  provide  Impact 
conditions  corresponding  to  the  Impact  parameters  which  have  been  used 
In  the  past  (Impingement  angle.  Impact  velocity,  particle  mass  and 
shape,  and  so  on)  to  correlate  the  experimental  determination  of  mass 
removal.  Concern  over  the  nature  of  the  particulate  flow  and  the  actual 
Impact  conditions  at  the  surface  of  the  target  has  been  expressed  by 
several  Investigations  In  relation  to  the  results  they  obtained  In  their 
erosion  facilities.  Since  this  problem  Is  beginning  to  receive  atten- 
tion, due  to  the  need  for  more  refined  and  reliable  erosion  measurements 
the  available  experimental  results  and  procedures  were  examined  but  no 
coherent  treatment  of  this  work  evolved  which  would  be  of  value  to 
current  and  future  Investigations. 


The  survey  of  the  solid  particle  erosion  testing  facilities  and 
procedures  Indicates  that  there  Is  a definite  need  for  devising  a 
facility  to  provide  an  accurate  measure  of  the  erosive  response  of 
materials  and/or  to  provide  detailed  charaterlzatlon  of  the  actual  Im- 
pact parameters  for  existing  facilities.  This  Is  not  a simple  under- 


taking  due  to  the  number  of  parameters  which  are  associated  with  mul- 
tiple collisions  of  Irregular-shaped  particles  on  a material  surface 
whose  general  topography  and  properties  are  changing  with  the  length 
of  the  exposure.  The  characterization  process  for  a single  particle 
type  and  material  Is  In  Itself  quite  extensive.  Many  factors  have  to  be 
evaluated  to  develop  accurate  measures  of  the  mass  of  particles  actually 
contacting  the  surface,  their  Impact  velocity  and  Impingement  angle, 
and  the  nature  of  the  particles  delivered  to  the  surface  compared  fd.th 
the  Initial  properties  and  sizes  of  the  particles  (especially  for  soft 
particulates).  Concentration  effects  are  obviously  important.  Particle 
concentration  can  Influence  the  extent  of  particle  shielding  at  the  sur- 
face (the  accumulation  of  particles  on  the  surface  which  prevents  the  in- 
coming particles  from  making  direct  contact  with  the  target's  surface), 
particle  rebound  characteristics  which  change  as  a function  of  exposure 
time  also  affect  the  trajectories  of  the  Incoming  particles,  and  the  oc- 
curence of  particle  embedding  In  the  target's  surface  modifies  the  mea- 
sured erosive  response  of  the  base  material.  Current  Interests  are  di- 
rected toward  evaluating  the  response  of  materials  at  elevated  tempera- 
tures and  In  corrosive  environments  which  introduces  additional 
complications  Into  an  already  complicated  evaluation  procedure. 

In  the  past  the  above  considerations  would  generally  have  been 
dismissed  as  being  overstated,  but  the  experimental  work  underway  In 
several  laboratories  Is  demonstrating  that  these  Issues  have  to  be  re- 
solved. The  statistical  nature  of  the  multiple  particle  erosion  pro- 
cess has  to  be  recognized,  so  It  may  not  be  possible  to  provide  a de- 
finitive result:  the  erosive  response  for  restricted  impact  conditions 
may  have  to  be  expressed  In  terms  of  bounds  on  the  measure  of  erosion. 
However  even  defining  reliable  bounds  on  the  magnitude  of  the  erosive 
response  which  are  independent  of  the  test  procedure  used  would  be  an 
accomplishment  in  view  of  the  present  state  of  affairs.  Therefore 
careful  examination  of  what  Is  really  required  from  a pragmatic  view- 
point should  be  undertaken  and  then  a general  methodology  for  achieving 
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a realistic  testing  objective  established.  It  Is  therefore  suggested 
that  a significant  contribution  to  the  solid  particle  erosion  field 
would  be  to  fully  characterize  existing  test  facilities  and  test  pro- 
cedures In  order  to  establish:  (1)  a set  of  design  requirements  for 
minimizing  the  Influence  of  the  test  apparatus  on  erosion  evaluations, 

(2)  operating  procedures  to  further  assure  that  the  test  parameters  are 
within  specified  limits,  and  (3)  reporting  procedures  which  will  supply 

adequate  Information  for  reproducing  the  test  conditions  in  other  labora- 

* 

torles  in  order  to  confirm  or  extend  the  experimental  data.  It  is 
almost  Imperative  that  such  as  investigation  be  a necessary  prerequisite 
to  any  large-scale  test  and  evaluation  program. 

5.2  Solid  Particle  Erosion  Modeling 

The  various  solid  particle  erosion  analyses  which  have  been  de- 
veloped have  been  reviewed  and  compiled  for  future  reference.  This  re- 
port does  not  replace  reference  to  the  original  work,  which  is  strongly 
advised,  but  attempts  to  highlight  the  basis  for  the  approach  adopted, 
provides  some  perspective  on  the  general  Interrelationships  among  the 
various  analyses,  and  re-examines  the  validity  of  the  analytical  pro- 
cedures employed.  During  the  course  of  presenting  the  analytical  work, 
the  solid  particle  experiments  sometimes  used  as  the  basis  for  the 
modeling  effort  have  also  been  examined. 

Hutchings'  work  has  made  a significant  contribution  to  providing 
a more  realistic  view  of  the  material  removal  process  due  to  oblique 
Impacts  of  both  rounded  and  angular  rigid  particles  on  ductile  materials. 
Reasonable  correlation  with  the  mass  loss  for  spherical  body  Impacts 
is  achieved  if  numerical  Integration  of  the  equations  of  motion  for 

the  progression  of  the  particle  through  the  Mterlal  replaces  the  re- 

_ 

Round-robin  tests  are  being  conducted  by  members  of  the  ASTM  Subcomnlttee 
G2.20  on  Erosion  by  Solid  Particle  Impingement  to  begin  to  assess  what 
can  be  done  to  address  these  problems.  If  these  criteria  are  to  be  es- 
tablished within  a reasonable  time  frame,  the  magnitude  of  this  effort 
will  have  to  be  considerably  expanded. 
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strlctive  assumptions  which  have  to  be  Imposed  to  analytically  Integrate 
the  equations  of  motion  (Flnnie,  1958;  Surette,  1971).  The  model  used 
by  Hutchings  Is  still  highly  Idealized  and  requires  further  modifications 
before  it  may  adequately  reproduce  the  volume  removal  rates.  However 
Hutchings'  work  demonstrates  the  benefits  of  a well-characterized  ex- 
perimental arrangement  for  gaining  Insights  Into  the  mechanics  of  the 
material  removal  process  and  In  clarifying  by  direct  observation  the 
effectiveness  of  the  modes  of  material  removal  possible.  This  latter 
consideration  Is  usually  a matter  of  conjecture  In  the  modeling  approaches 
(Flnnie,  1958;  Bitter,  1963a, b). 

Both  Sheldon  and  Kanhere  (1972)  and  Flnnie  and  McFadden  (1978) 
attempted  to  develop  analyses  which  would  provide  a velocity  exponent 
greater  than  t«K>.  It  is  this  Investigator's  finding  that  If  Sheldon 
and  Kanhere 's  analysis  is  correctly  carried  out,  a velocity  exponent  of 
two  results  Instead  of  three  as  they  conclude.  Furthermore,  It  Is  this 
Investigator's  observation  that  Flnnie  and  McFadden  artlflcally  Intro- 
duce the  desired  velocity  dependence  into  the  analysis,  so  the  model 
they  propose  Is  not  really  representative  of  the  outcome  of  their  de- 
rivation. Bitter's  analysis  (1963b)  Is  somewhat  overshadowed  by  the 
oft-quoted  simplification  of  the  correlation  presented  by  Nellson  and 
Gilchrist  (1968),  however  Bitter's  result,  although  based  on  heuristic 
arguments  yields  a velocity  exponent  of  5/2  or  less  when  a ^ a [refer 
to  Eq.  (3.63)]. 

Velocity  exponents  not  equal  to  two  were  the  result  of  Mamoun's 
analysis  (which  Is  an  elaboration  of  Bitter's  mechanics  analysis  which 
views  the  erosion  of  ductile  metals  as  a fatigue  process).  Head's 
statistical  analysis,  and  Tabakoff's  modifications  of  Nellson  and 
Gilchrist's  purely  empirical  correlations.  ' Mamoun's  published  work  Is 
only  applicable  to  normal  colllsons  of  spherical  particles,  however, 
the  introduction  of  empirical  fatigue  correlations  results  In  material 
dependent  velocity  exponents.  Several  assumptions  are  made  with  regard 
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to  the  volume  of  material  removed  and  In  correlating  the  experimental 
data  which  weakens  Mamoun's  contention  that  he  developed  a definitive 
predictive  theory.  Head's  statistical  approach  suffered  In  part  from 
an  Inadequately  characterized  erosion  apparatus:  significant  differences 
were  found  In  the  erosion  data  when  a slightly  modified  test  arrangement 
was  used  (cp.  Head  and  Harr,  1970,  and  Head,  et  al.,  1973).  Although 
velocity  exponents  different  than  two  were  obtained,  the  inconsistency 
In  the  results  essentially  Invalidates  their  usefulness.  A much  better 
correlation  with  their  experimental  data  Is  claimed  than  can  be  supported 
by  the  comparisons  provided.  Head's  work  does  not  rule  out  the  use  of 
statistical  correlative  procedures,  as  Is  sometimes  suggested  In  the 
solid  particle  erosion  literature  since  It  has  been  attempted,  but 
simply  Illustrates  that  meaningful  results  cannot  be  obtained  from  a 
testing  program  which  does  not  adequately  address  the  prerequisites 
listed  at  the  end  of  Section  5.1.  Subsequent  analysis  along  the  same 
lines  (Jennings,  et  al.,  1976)  Indicated  the  velocity  exponent  had  a 
constant  value  of  five  which  does  not  correspond  to  the  available  test 
results  for  which  the  velocity  exponent  rarely  exceeds  three  for  ductile 
metals.  Tabakoff  investigated  the  statistics  of  the  particle  rebound 
characteristics  In  his  erosion  facility  which  he  described  empirically. 

Introduction  of  these  empirical  correlations  Inherent  In  his  erosion 

facility  Into  a modified  form  of  Nellson  and  Gilchrist's  volume  removal 

relations  also  yields  a velocity  exponent  greater  than  two  for  ductile 

metals  which  reaches  a constant  value  of  four  for  normal  xmpacts  which 

appears  to  describe  Tabakoff's  experimental  data.  The  strong  dependence 

of  the  erosion  data  on  the  size  of  the  specimen  In  Tabakoff's  wind  tunnel 

facility  has  to  be  resolved  In  order  to  differentiate  the  material's  I 

response  from  the  aerodynamics  of  particulate-laden  gases.  ; 

It  is  seen  that  the  physical  basis  for  the  correlations  (analyses) 
of  solid  particle  erosion  of  ductile  materials  Is  almost  negligible. 

It  may  be  that  the  process  Is  too  complex  to  be  described  by  other  than 
seml-empirical  correlations  of  the  test  data,  however  there  are  still 
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considerable  advances  which  can  be  made  in  this  regard.  The  best 

•> 

potential  for  progress  is  a closer  association  with  the  microscopic  de- 
tails of  the  particle/target  Interactions  responsible  for  material  re- 
moval as  observed  for  idealized,  but  well  characterized,  experimental 
conditions.  The  review  of  past  work  Indicates  that  this  approach  has 
not  yet  been  extensively  pursued  for  solid  particle  erosion  of  ductile 
metals.  Target  surface  melting  was  suggested  by  Smeltzer,  et  al.  (1970) 
futd  Jennings,  et  al.  (1976)  based  on  their  microscopic  observations  of 
eroded  metallic  surfaces.  The  significance  of  localized  melting  during 
solid  particle  impingements  is  still  a controversial  matter.  As  the 
number  of  investigators  examining  the  microscopic  details  of  eroded  sur- 
faces Increases,  the  role  of  localized  melting  should  be  clarified. 

The  approach  suggested  above  has  been  pursued  by  Evans  for  normal 
collisions  of  primarily  spherical  particles  on  polycrystalline  ceramics. 
A number  of  useful  observations  have  been  made  Including  an  estimate  of 
volume  removal  per  particle  Impact.  This  semi-emplrlcal  correlation 
indicates  an  explicit  velocity  exponent  of  3.17  for  all  ceramics  such 
that  the  Impact  conditions  are  within  the  elastic-plastic  response  re- 
gime: allowance  is  made  for  additional  velocity  dependence  in  the  pro- 
portionality coefficient,  Eq.  (4.13).  However  Evans  was  more  concerned 
with  material  properties  which  affect  the  extent  of  the  damage  due  to 
single  particle  impacts.  Fracture  toughness  and  hardness  are  identified 

as  the  most  significant  material  properties.  The  dependence  of  the 

-4/3  -1/4 

volume  removed  per  particle  on  K H , where  K is  the  fracture 

c c 

toughness  and  H is  the  Vickers  hardness  of  the  target  material,  shows 

moderate  correlation  with  the  experimental  data  (Evans,  et  al.,  1978) 

but  the  correspondence  leaves  room  for  improvement.  The  dependence  of 
-4/3  1/9 

H reported  by  Ruff  and  Wlederhom  (1979)  shows  little  cor- 
relation with  the  experimental  data. 
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Although  the  general  character  of  the  fractures  produced  In 
ceramics  by  quaslstatlc  loading  conditions  appear  similar  to  those  ob- 
served for  particle  collisions,  Evans'  quantitative  measurements  of  the 
extent  of  the  fracture  surfaces  show  the  dynamic  loading  conditions  pro- 
duce significantly  larger  cracks  than  quaslstatlc  loadings  for  comparable 
contact  zone  radii.  A transient  stress  analysis  Is  therefore  required  to 
adequately  estimate  the  extent  of  the  crack  growth  during  the  loading/ 
unloading  cycle. 

The  anlaytlcal  representations  for  volume  removal  from  brittle 
materials  are  quite  crude  In  relation  to  the  physical  process  of  removal 
[Eq.(4.1),  (4.7),  and  (4.9)],  and  as  defined,  would  overestimate  the 
volume  of  material  removed.  As  long  as  the  primary  contributors  to  the 
removal  process  are  Identified  and  are  assigned  the  proper  functional  de- 
pendence, this  volume  discrepancy  should  appear  in  the  proportionality 
coefficient.  All  of  the  brittle  erosion  analyses  assume  that  the  condi- 
tions for  the  Initial  particle  impact  prevail  during  the  erosion  process. 
That  this  was  not  the  case  was  shown  by  Adler  (1974b,  1976b)  for  a very 
Idealized  experimental  condition,  and  a counter  example  is  provided  by 
Gulden  (1979b).  The  particular  application  would  determine  If  conditions 
during  the  Incubation  period  would  be  most  significant  or  if  the  steady 
state  removal  rates  are  the  desired  erosion  measure.  In  most  applica- 
tions the  incubation  period  is  quite  short  and  not  generally  an  Impor- 
tant factor. 

There  Is  an  obvious  need  to  obtain  better  material  property  re- 
lations In  order  to  improve  the  predictive  capability  of  the  existing 
brittle  erosion  models  and  to  generalize  the  results  to  oblique  impacts 
(a  simple  generalization  was  proposed  by  Hockey,  et  al.  (1977)).  In 
contrast  to  the  erosion  models  developed  for  ductile  materials,  the 
recent  semi-emplrical  correlations  for  brittle  materials  Include  de- 
finable material  properties  for  the  target  materials.  These  correlations 
are  based  on  direct  observations  of  single  particle  Impact  damage  and 
Incorporate  concepts  from  fracture  mechanics  and  stress  wave  propagation 
In  deformable  media. 
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APPENDIX  A 


Hertzian  Theory  of  Impact 

A number  of  solid  particle  erosion  theories  employ  the  Hertzian 
theory  of  Impact,  althou^  the  limits  of  application  for  the  results 
from  this  theory  are  often  greatly  exceeded.  The  advantage  of  using  the 
Hertzian  theory  Is  that  explicit  expressions  for  the  stresses  In  an  elas- 
tic half-space  can  be  evaluated  with  very  simple  expressions  for  the 
surface  of  the  half-space  which  Includes  the  location  of  the  critical 
stress  condition  for  fractures  In  brittle  materials.  The  stresses  In  an 
elastic  half-space  Impacted  by  a spherical  body  with  different  elastic 
properties  than  the  half-space  will  be  considered  In  some  detail,  since 
there  does  not  appear  to  be  a readily  available  source  for  obtaining  ex- 
plicit expressions  for  both  the  stress  and  displacement  components. 

The  problem  originally  solved  by  Hertz  (1896)  was  the  evaluation 
of  the  stress  distribution  associated  with  two  Isotropic  elastic  bodies 
which  touch  each  other  over  a very  small  part  of  their  surface  and  ex- 
ert upon  each  other  a finite  pressure  over  the  common  area  of  contact. 
Hertz  determined  the  form  of  the  common  surface,  the  distribution  of 
pressure  over  It,  and  the  relation  between  the  distance  the  bodies  ap- 
proach each  other.  The  results  given  here  will  be  restricted  to  the 
case  of  a spherical  elastic  body  Indenting  an  elastic  half-space  which 
Is  adequate  for  the  analyses  In  the  erosion  literature.  A detailed  ac- 
count of  the  stress  distribution  for  this  case  was  given  by  Huber  (1904) 
and  the  principal  stress  trajectories  were  computed  by  Fuchs  (1913)  from 
Huber's  formulas  for  the  stress  components. 

The  radius  a of  the  circle  of  contact  and  the  distance  of 
approach  q for  an  elastic  sphere  indenting  a thick  elastic  plate  are 
given  by 


\ y.-' 


where  P is  the  normal  load  on  the  Indenter,  R is  the  Indenter  radius 
and 


where  Poisson's  ratio  and  Young's  modulus  for  the  sphere  and 

elastic  constants  for  the  half-space.  The  interfacial 
contact  pressure  between  the  sphere  and  the  half-space  is  found  to  have 
the  general  form 


The  non-vanishing  stress  and  displacement  components  in  the  half' 
space  for  the  Hertz— Huber  theory  can  be  evaluated  from  the  analysis 
found  in  a paper  by  Love  (1929).  The  notation  used  is  defined  in 
Figure  Al. 


(l+v)Ctan  (n/C)  + 2vn 


STRESS  COMPONENTS 
AT  POINT  (p,eA 
IN  THE  HALF-SPACE 


Figure  Al.  Notation  Used  In  Hertzian  Stress  Analysis 


Numerical  evaluations  of  Eq.(A5)  to  (AlO)  can  be  carried  out,  but 
the  explicit  representations  for  the  limiting  forms  of  these  equations 
at  the  surface  of  the  half-space  (z~0)  and  along  the  axis  of  symmetry 
(p-0)  will  be  displayed. 

vnien  C*0  t Eq.(All)  yields 

- (C^-1)  , 0 (A12) 


With  the  restriction  that  A Is  the  positive  root  of  Eq.(All),  two 
cases  must  be  considered 


A^  ■ when  ^ 1 

2 

A *0  when  5 < 1 


(A13) 

j 
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Then  from  the  definitions  of  the  dimensionless  parameters 


whereby 


when  C > 1 


when  C < 1 


Introducing  these  limiting  forms  into  EQ. (AS)  to  (AlO) , we  find  for 
C -0 


when  C^l 


when  ^>1 


ifhen  ^<1 


when  C>1 


when  C<1 


when  C^l 


when  C<1 


when  C>1 


when 


when 


The  limiting  forms  when  ^-0  follow.  Now  Eq.(All)  yields 


which  is  equivalent  to 


Using  this  condition  in  conjunction  with  the  limit  as  ^-*0  in  Eq«(A5} 
to  (A8), 


+ (l+v)Ctan  (1/0-2 


(l+v)Ctan  (1/0+2V 


Using  either  Eq.(A16)  to  (A18)  and  letting  • or  using 

Eq.(A24)  to  (A26)  and  letting  (-*0  , we  have  the  following  limiting 
cases  at  the  initial  point  of  contact  . 


(l+2v) 


(l+2v) 


The  corresponding  displacements  are 


The  generally  quoted  results  follow  from  the  limiting  cases  for 
stress  components.  From  Eq.(A16)  and  (A17)  it  is  seen  that  the  circum- 
ferential stress  component,  Oqq  , is  compressive  over  the  entire  sur- 
face and  of  the  same  magnitude  as  for  p>a.  The  a stress  coai- 

ponent  is  compressive  and  corresponds  to  the  loading  condition  specified 

in  Eq.(A4).  According  to  Eq.(A24)  to  (A26),  and  O are  all 

pp’  00’  zz 

compressive  along  the  z-axls  (axis  of  symmetry)  between  the  surface  and 

a critical  depth  which  depends  on  V . At  this  depth  both  a and 

PP 

^96  vnnish  and  are  tensile  for  all  greater  depths.  The  o stress 
component  is  always  compressive  along  the  axis  of  symmetry.  The  maxi- 
mum radial  tensile  stress  is  located  at  the  circular  boundary  of  the 
contact  surface  which  according  to  Eq.(A16)  has  the  magnitude 


The  general  stress  and  dlsplacesient  components  in  Eq.(A16)  to  (A21)  can 
be  readily  programed  for  digital  computer  evaluation  at  any  point  in  the 
half-space. 


Hertz  assumed  that  the  collision  of  two  bodies  could  be  regarded 
•s  a statical  problem;  the  compression  at  the  place  of  contact  is  re- 
garded as  gradually  applied  and  as  subsiding  coaq>letely  by  reversal  of 
the  process  by  which  it  was  produced.  To  study  the  impact  of  a spheri- 
cal body  on  an  elastic  surface,  the  previous  relations  from  the  Hertzian 


theory  are  restated  in  terns  of  the  Impact  velocity  Instead  of  the 

load  P which  is  now  an  unknown.  In  order  to  establish  a relation  be- 
tween P and  , the  kinetic  energy  of  the  impacting  sphere  is  equated 
to  the  strain  energy  at  the  time  of  maximum  impression. 

qmax 

1(3  *^1^  '^o  " / (A34) 

o 


where  is  the  density  of  the  impacting  sphere. 


P 

max 


-2/5  6/5 

K V 

o 


2 

R 


Using  Eq.(A2) 

(A35) 


The  duration  of  contact  T between  the  two  impacting  bodies  is  given 
by 


T 


2.943 


(!"v) 


2/5 


(A36) 


Thus  knowing  the  elastic  properties  of  the  colliding  bodies,  the  load 
P in  Eq.(A16)  to  (A21)  can  be  replaced  by  Eq.(A35)  in  terms  of  the  im- 
pact velocity  so  the  stress  and  displacement  components  can  be 

evaluated  for  the  Hertzian  Impact  conditions. 
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APPENDIX  B 


Geometrical  Relations  for  a Rigid  Sphere 
Indenting  a Compressible  Half-Space 


The  volume  of  material  removed  from  a material  surface  due  to 
particle  Impingement  Is  often  related  to  the  volume  Included  within  a 
spherical  segment.  The  geometrical  relations  for  a spherical  segment 
are  summarized  for  general  reference.  The  geometry  of  the  Indent  Is 
shown  In  Figure  Bl. 


The  approximation  to  d can  be  used  when  the  depth  of  the  Indent  Is 
small. 


The  surface  area  of  the  Indent  Is 


where  the  geometric  relation  In  Eq.(Bl)  was  used 


The  volume  of  the  spherical  Indent  is  determined  from 


using  Eq.(B2) 
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